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FOREWORD 


This  report  was  prepared  by  Wyle  Laboratories,  Scientific  Services  & Systems  Group, 
for  the  George  C.  Marshall  Space  Flight  Center,  National  Aeronautics  and  Space 
Administration.  The  work  was  performed  under  contract  NAS8-33508,  entitled  "Dynamic 
Analysis  of  SSME  Vibration  and  Pressure  Data." 

Mr.  T.  Coffin  and  Mr.  W.  L.  Swanson,  of  the  Wyle/Huntsville  Research  Department, 
served  as  Program  Manager  and  Project  Engineer,  respectively,  on  this  study. 
Messrs.  T.  Gardner,  R.  Dandridge,  and  Dr.  J.  Jong  provided  valuable  contributions 
to  SSME  data  analysis  and  component  modeling  tasks.  Mr.  W.  C.  Smith,  MSFC/ED24 
served  as  Contracting  Officer’s  Technical  Representative  for  the  study.  In  this 
capacity  he  provided  valuable  assistance  in  the  coordination  of  test  evaluation 
activities  and  served  as  a focal  point  for  defining  task  requirements  and  priori- 
ties. 


ABSTRACT 

This  report  summarizes  an  engineering  evaluation  of  dynamic  data  from  SSME  hot 
firing  tests  and  SSV  flights.  The  basic  objective  of  the  study  was  to  provide 
analyses  of  vibration,  strain  and  dynamic  pressure  measurements  in  support  of  NKFC 
performance  and  reliability  improvement  programs.  A brief  description  of  the  SSME 
test  program  is  given  and  a typical  test  evaluation  cycle  reviewed.  Data  banks 
generated  to  characterize  SSME  component  dynamic  characteristics  are  described  and 
statistical  analyses  performed  on  these  data  base  measurements  are  discussed. 
Analytical  models  applied  to  define  the  dynamic  behavior  of  SSME  components  (such 
as  turbopump  bearing  elements  and  the  flight  accelerometer  safety  cut-off  system) 
are  also  summarized.  Appendices  are  included  to  illustrate  some  typical  tasks 
performed  under  this  study. 
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SECTION  I 

EXECUTIVE  SUMMARY 

The  Space  Shuttle  Main  Engine  (SSME)  system  and  components  have  been  and  are 
presently  undergoing  extensive  development  and  certification  tests,  at  which  time 
some  250  measurements  are  acquired,  including  vibration,  strain,  and  dynamic 
pressure  at  critical  engine  locations.  Limited  engine  data  is  also  recorded  during 
operational  flights.  Under  the  severe  temperature,  pressure,  and  dynamic  environ- 
ments sustained  during  operation,  engine  systems  and  components  have  been  subject  to 
malfunction  and  failure.  Over  the  past  13  years  of  SSME  development,  28  major 
component  failures  have  occurred,  causing  extensive  damage  to  engine  hardware  and 
test,  facilities,  at  considerable  expense  in  cost  and  schedules.  In  addition, 
numerous  off  normal  operations  of  a less  severe  nature  have  been  observed. 

Detailed  analysis  and  evaluation  of  the  dynamic  measurements  obtained  is  mandatory 
to  permit  quick  assessment  of  engine  component  condition  and  the  initiation  of 
corrective  measures  when  necessary.  A wealth  of  dynamic  data  is  available  to 
support  this  effort.  Efficient  performance  of  this  task  is  especially  critical  when 
considering  the  high  test  and  launch  rate  in  progress  and  the  significant  impact  of 
dynamic  evaluation  results  on  test  and  launch  turn-around  time.  Dynamic  analysis, 
modeling  and  evaluation  of  SSME  measurements  and  components  have  been  performed  by 
Wyle,  under  NASA  contract  NAS8-33508,  in  support  of  these  efforts.  This  report 
documents  the  results  of  this  investigation. 

The  basic  objective  of  the  effort  was  to  provide  detailed  analyses  and  evaluation  of 
vibration,  dynamic  pressure,  and  strain  data  available  and  being  acquired  during 
SSME  tests.  Additionally,  analyses  were  performed  on  data  obtained  from  Space 
Shuttle  Vehicle  flights.  Statistical  analyses  were  performed  to  characterize 
nominal  and  abnormal  SSME  component  dynamic  behavior  under  available  operating 
conditions.  Data  banks  representing  component  characteristics  were  generated  and 
updated,  and  analyses  performed  to  assess  system  response  under  actual  and  hypothe- 
sized operating  conditions.  This  included  update  of  the  MSEC  Diagnostic  Data  Base 
and  application  of  the  SSME  Isospectral  Automated  Data  Base  System.  A data  base 
was  also  generated  and  updated  in  a format  compatible  with  the  Flight  Accelerometer 
Safety  Cut-Off  System  (FASGOS)  filter  characteristics  to  provide  a basis  for 
comparison  of  static  test  and  flight  results.  FASCOS  operational  characteristics 
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have  recently  been  defined  statistically  through  system  analysis  and  simulation,  to 
permit  assessment  of  the  system  logic. 

The  work  was  performed  under  three  broad  tasks,  which  are  summarized  as  follows. 

TASK  I:  Analysis,  Evaluation,  and  Documentation  of  SSME  Dynamic  Test  Results 
Under  this  task  Wyle  performed  analysis,  evaluation  and  documentation  of  SSME 
dynamic  test  results.  This  task  represented  the  mainstream  of  the  contract  effort 
and  included  data  verification,  analysis,  evaluation,  and  documentation  for  each 
SSME  ground  test  and  additionally  for  SSV  flight  measurements.  Results  included 
definition  of  temporal  and  spectral  characteristics  observed.  Subsynchronous, 
synchronous,  and  higher  order  spectral  characteristics  were  summarized  for  engine 
components  under  all  available  operating  conditions  to  characterize  SSME  component 
behavior.  Full  utilization  of  the  SSME  Isospectral  Automated  Data  Base  System  was 
employed  to  provide  informative  data  summaries.  Strain  measurements  were  analyzed 
as  above,  and  improved  strain  data  reduction  procedures  evaluated  for  subsequent 
strength  and  fatigue  analysis.  Oral  presentation  and  written  summaries  were 
provided  to  document  the  results  of  each  test  and  flight. 

TASK  III:  Development  and  Documentation  of  Statistical  Models  of  SSME  Component 
Dynamic  Behavior 

Under  this  task  Wyle  developed,  maintained,  and  updated  the  data  base  and  statis- 
tical models  of  SSME  component  dynamic  response  measurements  to  provide  characteriz- 
ing profiles  of  observed  parameter  ranges,  distributions,  etc.,  under  normal  and 
abnormal  operating  conditions  at  available  power  levels.  These  models  included 
wide-band  root-mean-square  values,  narrow-band  spectra,  and  band-pass  results  in 
the  flight  data  format.  Analytical  models  and  computer  schemes  were  applied  to 
define  the  dynamic  and  statistical  behavior  of  SSME  components,  such  as  bearing 
elements  and  the  Flight  Accelerometer  Safety  Cut-off  Systems. 

The  above  tasks  are  seen  to  be  intimately  related,  since  promising  statistical/ 
analytical  models  may  be  immediately  integrated  into  the  SSME  evaluation  process. 
Also,  data  base  statistics  were  immediately  updated  as  test  and  flight  measurements 
became  available.  It  should  be  noted  that  the  above  evaluations  were  performed 
under  extremely  limited  time  constraints  consistent  with  test  and  flight  schedules. 
The  extent  of  a given  investigation  varied  widely  depending  on  the  specific  measure- 
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ments  acquired,  whether  or  not  observed  engine  operation  was  nominal,  and  the 
severity  of  any  observed  component  malfunction  of  failure. 

The  following  section  presents  an  overview  of  the  study.  The  SSME  test/measurement 
program  is  briefly  described.  Dynamic  evaluation  considerations  are  discussed. 
SSME  data  base  development  and  application  tasks  performed  under  this  contract  are 
summarized  along  with  analytical/statistical  modeling  efforts.  A list  of  engineer- 
ing investigations  recently  performed  in  support  of  SSME  dynamic  evaluations  is 
presented.  This  table  indicates  the  wide  diversity  of  engineering  effort  required 
in  accomplishing  contract  objectives. 

As  with  most  test/evaluation  programs,  quick  turnaround  of  investigative  results 
was  imperative  to  successful  task  accomplishment.  These  results  were  therefore 
provided  directly  to  the  MSEC  COTR  and  cognizant  Program  personnel  through  informal 
data  packages  and  presentations.  A detailed  discussion  of  each  evaluation  performed 
is  given  in  the  technical  progress  and  interim  reports  provided  under  this  contract. 
Several  of  these  reports,  addressing  specific  task  evaluations,  are  included  as 
Appendices  for  reference. 
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SECTION  II 

PROGRAM  REVIEW  AND  SURVEY  OF  TASKS 
2.1  Background 

The  Space  Shuttle  Main  Engines  (SSME)  are  required  to  operate  under  extreme  tempera- 
tures with  high  fluid  pressures  and  rotational  pump  speeds.  Developmental  work  is 
presently  in  progress  to  uprate  SSME  performance,  including  engine  certification  at 
FPL  (109%).  The  SSME  and  components  have  been  subjected  to  extensive  hot  firing 
tests.  Acceleration,  dynamic  pressure  and  strain  data  have  been  acquired  during 
these  tests  and  additionally  from  Space  Shuttle  Vehicle  (SSV)  flights.  Analysis 
and  evaluation  of  results  obtained  from  SSME  operation,  to  aid  in  the  identification 
and  resolution  of  sources  of  malfunction,  were  performed  by  Wyle  Laboratories  under 
NASA  Contract  NAS8-33508.  Dynamic  and  statistical  modeling  to  assess  system 
behavior  and  component  condition  has  also  been  accomplished. 

This  section  presents  an  overview  of  the  Space  Shuttle  Vehicle  (SSV)  system  and 
SSME  in  particular.  The  SSME  test  program  is  briefly  described,  and  a typical  SSME 
test  evaluation  cycle  is  outlined.  Some  methods  applied  by  Wyle  to  the  assessment 
of  acceleration,  strain  and  dynamic  pressure  measurements  are  discussed,  A recent 
investigation  of  the  FASCOS  system  is  also  described. 

2.1.1  The  Space  Shuttle  Vehicle  System 

The  SSV  is  composed  of  the  Orbiter,  an  External  Tank  (ET),  which  contains  the  ascent 
propellant  to  be  used  by  the  Orbiter’s  three  main  engines,  and  two  Solid  Rocket 
Boosters  (SRB).  The  Orbiter  and  SRBs  are  reusable;  the  ET  is  expended  on  each 
launch. 

A Space  Shuttle  mission  begins  with  installation  of  the  mission  payload  into  the 
Orbiter  cargo  bay.  The  SRBs  and  the  SSMEs  fire  together  at  liftoff.  The  two  SRBs 
are  jettisoned  after  burnout  --  about  45  kilometers  (28  miles)  high  — and  recovered 
for  reuse  by  means  of  a parachute  recovery  system.  The  SSMEs  continue  to  burn  until 
the  Orbiter  is  just  short  of  orbital  velocity,  at  which  time  the  engines  are  shut 
down  and  the  ET  jettisoned.  During  its  return  through  the  atmosphere,  the  tank  will 
tumble,  break  up  and  be  destroyed. 

The  orbital  maneuvering  system  is  used  to  attain  the  desired  orbit  and  to  make  any 
subsequent  maneuvers  that  may  be  needed  during  a mission.  After  orbital  opera- 
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tions  are  completed,  normally  about  seven  days,  deorbiting  maneuvers  are  initiated. 
The  Orbiter  reenters  the  Earth’s  atmosphere  at  a high  angle  of  attack.  It  then 
levels  into  horizontal  flight  at  low  altitude  for  an  unpowered  aircraft-type 
approach,  landing  at  a speed  of  about  335  kilometers  per  hour  (208  miles  per  hour). 

2.1.2  The  Space  Shuttle  Main  Engines 

The  Orbiter  vehicle  main  propulsion  system  consists  of  three  SSMEs.  The  SSMEs  are 
reusable,  high-performance,  liquid-propellant  rocket  engines  with  variable  thrust. 
They  are  ignited  on  the  ground  at  launch  and  operate  in  parallel,  with  approximately 
500  seconds  total  firing  duration.  Each  of  the  rocket  engines  operates  at  a mixture 
ratio  (liquid  oxygen/liquid  hydrogen)  of  6:1  and  a chamber  pressure  of  approximately 
3000  psia  to  produce  a sea-level  thrust  of  375,000  pounds  and  a vacuum  thrust  of 
470,000  pounds.  The  engines  are  presently  throttleable  over  a thrust  range  of  60  to 
109  percent  of  the  design  thrust  level.  This  provides  a higher  thrust  level  during 
liftoff  and  the  initial  ascent  phase,  and  allows  Orbiter  acceleration  to  be  limited 
to  3 g’s  during  the  final  ascent  phase.  The  engines  are  gimbaled  (+10.5  degrees  for 
pitch  and  ±8.5  degrees  yaw)  to  provide  pitch,  yaw,  and  roll  control  during  the 
Orbiter  boost  phase. 

Significant  to  meeting  performance  requirements  is  the  use  of  a staged  combustion 
power  cycle  coupled  with  high  combustion  chamber  pressures.  In  the  SSME-staged 
combustion  cycle,  the  propellants  are  partially  burned  at  high  pressure  and  rela- 
tively low  temperature  in  the  preburners,  then  completely  combusted  at  high  tempera- 
ture and  pressure  in  the  main  chamber  before  expanding  through  the  high-area-ratio 
nozzle.  Hydrogen  fuel  is  used  to  cool  all  combustion  devices  in  contact  with  high- 
temperature  combustion  products.  An  electronic  engine  controller  automatically 
performs  checkout,  start,  mainstage,  and  engine  shutdown  functions.  Major  compo- 
nents of  the  SSME  are  illustrated  in  Figure  1.  A more  detailed  view  of  the  SSME 
power  head  is  shown  in  Figure  2.  This  figure  provides  an  indication  of  the  complex- 
ity of  the  SSME  turbomachinery. 
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FIGURE  1.  Space  Shuttle  Main  Engine 


2.1.3  SSME  Development  and  Acceptance  Testing 
To  validate  system  performance  and  ensure  equipment  reliability,  the  SSME  and 
components  have  been  and  are  presently  undergoing  extensive  development  and  qualifi- 
cation tests.  Testing  of  the  engine  and  components  is  conducted  at  several  NASA  and 
contractor  locations.  Full  scale  engine  test  firings  for  development  and  flight 
acceptance  are  performed  on  two  single-engine  test  stands  at  the  National  Space 
Technology  Laboratories  (NSTL),  Bay  St.  Louis,  Mississippi,  and  at  one  stand 
operated  by  Rockwell  International  near  Santa  Susana,  California,  with  plans  to 
include  a development  test  stand  at  MSFC.  In  addition,  main  propulsion  testing 
(MPT)  is  performed  at  NSTL  on  a stand  designed  to  accommodate  the  Shuttle  main 
propulsion  system  elements~the  three-engine  cluster,  the  ET,  and  the  Orbiter 
systems. 
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Testing  is  being  performed  on  a continuing  basis.  The  length  of  a given  test  is 
dependent  on  specific  test  objectives  and  may  run  from  several  seconds  to  over  800 
seconds.  Tests  are  generally  designed  to  satisfy  multiple  specific  objectives, 
which  fall  into  two  broad  categories;  acceptance /certification  firing  of  flight 
hardware  and  development  testing  directed  towards  design  verification,  performance 
and  reliability  improvement.  Test  operations  are  controlled  by  a computer  called 
the  Command  and  Data  Simulator  (CADS)  which  communicates  with  the  engine,  displays 
vital  measurements  for  on-line  observation/control  and  initiates  pre-  and  post-test 
procedures. 

Approximately  250  measurements  are  recorded  on  a given  test  including  wide  band 
vibration,  dynamic  pressure  and  strain  at  critical  engine  locations.  Some  of  these 
measurements  are  utilized  on-line  as  emergency  cut-off  indicators  and  all  are 
recorded  on  magnetic  tape  for  subsequent  analysis  and  evaluation.  Limited  SSM£ 
vibration  measurements  are  recorded  on  magnetic  tape  during  SSV  flights  for  evalua- 
tion with  orbiter  return.  Typical  dynamic  measurements  obtained  during  SSME 
operation  are  illustrated  in  Figure  3. 

2.2  Engineerina  Evaluation  of  SSME  Dynamic  Data 

2.2.1  Dynamic  Analysis  and  Evaluation  Considerations 

Acceleration  measurements  are  obtained  at  fuel  and  oxidizer  turbopump  locations 
during  all  test  firings,  providing  an  extensive  vibration  data  base  representing 
various  turbopump  builds  under  widely  differing  operating  conditions.  Additional 
measurements  are  obtained  on  a test-specific  basis,  depending  on  performance, 
structural  integrity,  or  rotor  dynamic  characteristics  under  evaluation.  For 
example,  test  series  have  been  performed  with  some  80  strain  measurements  to  support 
engine  nozzle  and  injector  dynamic  stress  evaluations.  Recent  firings  have  also 
been  conducted  with  internally  instrumented  turbopumps  to  define  component  dynamic 
load  and  signature  characteristics.  Thus  it  is  seen  that  the  extent  of  the  evalua- 
tion process  varies  widely  from  test  to  test,  even  though  engine  performance  is 
nominal.  In  the  event  of  anomalous  performance  or  component  malfunction,  the  extent 
of  this  process  is  increased  significantly.  Limited  turbopump  measurements  are  also 
obtained  from  the  three  SSMEs  on  each  SSV  flight.  Data  bandwidths  available  from 
SSV  flight  instrumentation  differ  from  the  wide-band  capability  used  during  ground 
testing,  thus  the  need  for  generating  a data  base  of  filtered  ground  test  measure- 
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ments  to  permit  direct  comparison  with  flight  results.  Typical  engineering  activ- 
ities supported  by  Wyle  in  an  SSME  test  evaluation  cycle  are  summarized  as  follows. 

Data  evaluation  and  documentation: 

Data  verification  and  validation 
Events  analysis 

Temporal,  and  spectral  correlation  with  operating  profile  and  machine 
dynamic  characteristics 
Test/flight  data  summary 
Analytical/statistical  modeling  and  classification: 

Develop  statistical  models  characterizing  normal  and  abnormal  behavior 
Update  SSME  diagnostic  data  base 
Establish  and  update  engine  cut-off  redlines 

Develop  and  apply  computer  programs  to  define  SSME  component  dynamic 
behavior 

Failure  investigation: 

Time/event  correlations  with  other  test  parameters/observations 

- Temporal  and  spectral  comparison  with  structural  dynamic  and  statisti- 
cal models  and  associated  failure  modes 

- Evaluation  and  recommendations  of  probable  failure/effect  scenar- 
ios, and  means  of  resolution. 

It  should  be  noted  that  the  above  evaluation  must  be  performed  under  extremely 
limited  time  constraints  consistent  with  test  and  flight  schedules.  Also,  the 
extent  of  a given  evaluation  depended  on  the  specific  measurements  acquired  and 
whether  or  not  observed  engine  operation  was  nominal. 

2.2.2  SSME  Data  Base  Development  and  Application 

Wyle  personnel  have  been  instrumental  in  the  development,  modification  and  applica- 
tion of  the  MSFG  Diagnostic  and  Isospectral  data  base  programs.  These  routines  are 
used  extensively  in  routine  test  evaluation  and  also  in  diagnostic  investigations. 
The  SSME  diagnostic  data  base  and  software  implemented  in  the  last  year  has  greatly 
facilitated  the  generation  of  quick-look  performance  summaries  and  comparisons  for 
input  to  the  SSME  data  reviews  conducted  immediately  after  each  test,  as  well  as  the 
maintenance  of  historic  statistical  profiles.  Figure  4,  from  Wyle  TM  64058-01, 
illustrates  a statistical  summary  of  measurements  obtained  from  the  oxidizer 
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preburner  pump  during  141  test  firings  and  the  fuel  turbopump  flight  data  as  of 
August  1985.  These  data  represent  measurements  obtained  under  nominal  engine 
operating  conditions.  Separate  statistical  definition  is  given  for  wide-band 
(composite)  levels  and  the  spectral  levels  at  pump  shaft  speed  (synchronous)  and 
selected  harmonics.  Figure  5 illustrates  cumulative  and  probability  distributions 
characterizing  the  measurements.  The  cumulative  distribution  provides  a useful  tool 
for  the  quick-look  assessment  of  test  results  with  the  historical  distribution  of 
measurements  observed  over  a number  of  test  firings,  and  also  indicates  the  statis- 
tical risks  implied  by  vibration  level  redlines. 


A plot  of  the  classical  gamma  distribution  is  included  in  these  figures  for  refer 
ence,  with  parameters  defined  by 


r*(G) 


where 


and 


k 

a 

in 


G 

ik°^  /riot) ) J e-  dx 

0 


sample  mean 
sample  variance 


The  application  of  classical  statistical  models  is  desirable  for  data  characteriza- 
tion since  this  permits  continuous  statistical  definition  and  manipulation  from 
discrete  measurement  observations.  Automated  optimal  function  definition  has  also 
been  developed  for  implementation  in  the  diagnostic  data  base  program.  Software  has 
been  written  to  assess  a wide  class  of  competing  functions  including  the  Weibull, 
Gumbel  and  generalized  lambda  distributions. 


II-9 


STATISTICAL  SUHHARY  OF  SSHE  OIBRATIOH  DATA  16-MAR-84 


Test 

LOX  PBP  RAO 

i e 

45 

Max 

G 

LOX  PBP  RAO 

i 5 

135 

-1 

Max 

G 

LOX  PBP  RAD 

i 1 

135 

-2 

Max 

G 

Stand 

Tests 

rns 

Sis 

rns 

Tests 

rns 

SI0 

rns 

Tests 

rns 

Sis 

rns 

— 

— — 

eeeeee 

eweeee 

— 

A1 

41 

6.4 

2.8 

15.5 

47 

5.2 

2.1 

18.5 

38 

6.3 

2.2 

13.3 

A2 

44 

5.3 

2.3 

11.0 

52 

5.3 

1.6 

18.0 

24 

9.5 

1.9 

9.8 

A3 

41 

7.2 

2.4 

15.0 

42 

7.5 

2.8 

12.0 

26 

5.0 

1.8 

10.8 

Conbined  126 

6.3 

2.6 

15.5 

141 

5.9 

2.2 

12.8 

88 

5.6 

2.0 

13.3 

Syikchrotkous  81095C  Pou«r  Level  . 

LOX  PBP  RAD  45  LOX  PSP  RAD  135-1  LOX  PBP  RAD  135-2 

Max  _ Max  - Max 

Test  #6.  G # G ^ ® ^ ^ 

Stand  Tests  rfis  Si«  r«s  Tests  r«s  Sis  »'«*  Tests  r«s  Sis  rns 


40  2.9  1.8  7.S  44  3.2  2.0  8.6  31  3.5  1.9  8.5 

39  2.5  1.6  6.4  45  3.8  1.4  6.8  24  3.3  1.9  6.8 


37  2.5  1.6  7.6 


48  2.8  1.5  8.9 


24  2.5  1.8  9.1 


Conbined  116  2.6  1.7  7.6  129  3.0  1.7  8.9  79  3.1  1.9  9.1 

TEST  •'!  ^244-390, 81353-355*01357-362, A1364-388»A1398-392»A1394-398»  , 

A1 40 1 %08  * A 1 4 1 0-436 , A2 1 95-305 , A2307-332 * A3847- 1 12 » A3 1 1 4- 1 52  * A3194- 160*03 
169-172*03174-181*03184-223*03225-232, 

CoNooaite  0l04X  Po»»r  Level 

FUEL  PUnP  RAO  0 FUEL  PUMP  RAO  80  FUEL  PU1P  RAO  174 


m 

c 

rmm 

!!? 

riox 

s 

r»« 

■ 5 

Tea  t a me 

!!? 

Max 

G 

rma 

TM?a 

e 

raa 

S19 

flax 

C 

rmm 

STS 

24 

3*0 

1 .0 

5.0 

3 t.e 

t .0 

3.0 

!8 

3.4 

1 .5 

7.3 

STATISTICAL  SUnriARY  OF  SSflE  VIBRATION  BATA  2S-AUC-85 
Coepeeite  t104X  Power  Level 


FUEL  TURB  RAB  80 


FUEL  TURB  AXIAL 


Teat 

Stand 

• 

Teste 

e 

me 

Sis 

Max 

G 

rmm 

« c 

Immtm  rmm 

Sio 

flax 

C 

rmm 

m & 

Taafs  raa 

SlQ 

STS 

27 

S.8 

1.0 

5.2 

8 0.0 

0.0 

0.0 

0 0.0 

8.0 

Syrettronoue  •IB4X  Power  Level 

FUEL  PUflP  RAB  0 Fuel  pump  RAB  OB  fuel  PUNP  rad  174 

_ No*  _ Hex  Hox 

Tcei  C *C  C «C  C 

Stand  Teeie  rme  S19  me  Teaie  rme  Sig  me  Ta«i«  rme  Sig  me 


STS 

24 

2.7 

1.0  4.7 

S IS  1.0  2.3 

V8 

2.8 

0.8 

3.0 

STATISTICAL  SUnnARY  OF  SStlE  VIBRATION 

SATA 

23- 

-Auc-es 

Syachronoua  0104X  Poaar  Laval 

ruEL 

PUflP 

PAS  188 

fuel  TURB  RAB  80 

FUEL 

TUPS 

AXIAL 

Max 

_ Max 

flax 

Taat 

m 

C 

C 

e C C 

8 

e 

C 

Siond 

Taaia 

rmm 

SI9  raa 

Taava  raa  Sig  raa 

Taa?a 

raa 

Sio 

raa 

STS 

27 

2.5 

1.0  5.0 

0 0.0  0.0  0.0 

0 

0.0 

0.0 

0.0 

Figure  4.  Statistical  Summary  of  High  Pressure  Fuel  and 
Oxidizer  Pump  Measurements  for  Ground  and  Flight  Tests 

II-IO 


Of  POOR 


flLltl  CK  U lAi  X « O Ul 


CoNposite  9109^:  PMR  LUL  l6-nflR-84 


— — Sjmei^anou*  IMX  Wl*  LVU  87^AUC-i8 
TESif  •'%  A1183-489. >2221^374. *3864-280. 

Nu!M»*r  ef  tpsta  atlOS 


Figure  5.  Cumulative  and  Probability  Distribution  Plots 
for  a Selected  Vibration  Data  Sample 

ii-n 


In  addition  to  the  SSME  diagnostic  data  base,  the  SSME  Isospectral  Data  Base  System 
was  applied  extensively.  Historically,  analog  tapes  were  flown  to  MSFC,  dubbed, 
reduced,  and  evaluated.  Under  the  new  system,  spectra  are  extracted  from  the  tapes 
every,  0.4  second  for  each  measurement  throughout  each  test  and  stored  at 
NSTL/Slidell.  On  command,  the  spectral  data  is  telemetered  to  MSFC  via  satellite. 
Based  on  MSFC-developed  software,  these  preprocessed  data  are  then  manipulated  and 
printed  to  display  isoplots,  bandpass  trends,  engine  speed,  etc.,  on  user  command. 
As  the  diagnostic  data  base  system  has  made  test-to-test  comparisons  more  efficient, 
the  SSME  isospectral  analysis  system  has  been  a significant  aid  in  the  evaluation  of 
time/amplitude/frequency  trends  within  a test.  In  addition  to  the  above  preprocess- 
ed data,  analog  tapes  are  also  obtained  by  NCFC  for  special-purpose  time  and 
frequency  domain  analysis.  Modification  of  the  Isospectral  data  base  is  under 
evaluation  to  include  transmission  of  the  Fourier  coefficients  (amplitude  and  phase) 
as  well  as  PSD  samples.  This  will  permit  immediate  access  to  the  vibration  time 
history  for  detailed  analysis,  by  straightforward  inverse  transformation. 

Figure  6 is  an  example  of  an  isoplot  representing  a high  pressure  fuel  turbine 
measurement  throughout  test  902-378.  A frequency  band  of  60  Hz  to  1310  Hz  and  1310 
to  2560  Hz  was  selected  and  a spectrum  plotted  every  8 seconds.  The  displayed 
amplitude  range  is  selectable,  permitting  clear  representation  of  major  spectral 
peaks  or  identification  of  low  level  spectral  components.  Interpretation  obviously 
requires  correlation  with  engine  speed  and  other  parameters.  Figure  7 illustrates 
the  root-mean-square  acceleration  time  history  composite,  synchronous  and  selected 
harmonics  for  the  same  measurement.  This  time  history  is  synthesized  from  the 
stored  isospectral  data  by  integrating  over  the  PSDs  obtained  each  0.4  second  during 
the  test. 
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A sample  PSD  for  the  same  measurement  is  illustrated  in  Figure  8.  An  "Anomalous 
Frequency"  routine  was  developed  to  track  spectral  components  within  a selected 
frequency  band  and  relate  these  to  the  synchronous  behavior.  Figures  9 and  10 
illustrate  the  output  of  the  procedure  in  the  selectable  frequency  bands  of  200-400 
Hz  and  1700-2000  Hz.  Another  feature  of  the  Isospectral  Data  Base  Program  is  the 
cumulative  speed  analysis  routine.  Certain  components  on  the  SSME  respond  more 
strongly  when  exposed  to  a narrow  speed  range.  Software  was  developed  to  track  and 
evaluate  the  total  time  a particular  serial  number  pump  has  accumulated  in  each 
speed  range  for  a given  test  or  series  of  tests.  An  example  plot  and  tabulated 
values  for  two  recent  tests  (Al-377  and  Al-378)  are  shown  in  Figure  11. 

2.2.3  Special  Considerations  in  Dynamic  Strain  Assessment 

Wide  band  strain  measurements  have  generally  been  characterized  in  terms  of  power 
spectra  and  rms  time  histories.  In  Special  cases  (such  as  a LOX  injector  post 
fatigue  evaluation  supported  under  this  contract)  peak  strain  distributions  were 
estimated.  More  refined  cycle  counting  techniques  appeared  desirable  to  character- 
ize dynamic  strain  measurements.  A cursory  review  of  the  experimental  literature 
on  high  cycle  fatigue  was  therefore  performed. 

Miles’  classical  relation  for  fatigue  under  random  vibration  is  based  on  an  assumed 
distribution  of  stress  or  strain  peaks.  In  a complex  time  history,  consideration 
of  adjacent  peaks  (and  valleys)  only  can  mask  wide  variations  (the  range)  of  the 
dynamic  stress  or  strain  incurred.  A number  of  empirical  cycle  counting  techniques 
have  therefore  been  developed  by  researchers  to  better  account  for  the  ranges  of 
stresses  in  a complex  wave.  When  considering  that  damage  has  been  shown  to  be 
roughly  proportional  to  the  sixth  power  of  the  strain  range,  importance  of  cycle 
counting  and  range  definition  techniques  is  underlined. 

The  objective  of  all  cycle  counting  methods  is  to  relate  the  effect  of  a complex 
load,  stress  or  strain  time  history  to  material  fatigue  properties  obtained  under 
simple  cyclic  loading  conditions.  "Overall  range"  counting  methods  have  been  shown 
to  correlate  more  highly  with  observed  fatigue  life  than  simple  peak  value  distribu- 
tions. 
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Figure  8.  Power  Spectral  Density  of  a High  Pressure  Fuel  Turbopump  Measurement 


Three  cycle  counting  techniques  and  associated  computational  algorithms  appear  most 
prominent  in  the  literature  for  fatigue  evaluation  of  complex  strain  (stress)  time 
histories.  These  are  known  as  the  range-pair,  Rainflow  (or  hysteresis  loop  closure) 
and  ordered  overall  range  (or  racetrack)  methods.  These  three  procedures  yield 
virtually  identical  results.  Compu^t|^t|on^  schemes  for  implementation  of  these 
methods  were  also  obtained  from  the  references.  These  techniques  permit  significant 
data  compression,  in  that  only  major  strain  excursions  need  be  retained  to  charac- 
terize fatigue  growth.  The  strain  ranges  generated  as  above  are  generally  used 
concurrently  with  a fatigue  analysis.  However,  the  techniques  may  also  be  used 
directly  in  the  data  reduction  procedure  to  generate  strain  range  spectra  depicting 
the  number  of  occurrences  as  a function  of  strain  range  (a  probability  density 
approximation).  Statistical  definition  of  strain  range  spectra  per  event  (such  as 
engine  startup)  or  per  unit  time  (during  constant  power  operation,  for  example) 
should  constitute  a valuable  component  fatigue  life  estimation  tool. 

2.3  Analvtical/Statistical  Modeling 

The  above  discussion  illustrates  some  capabilities  of  the  SSME  data  base  systems 
and  software.  Wyle  personnel  have  made  extensive  use  of  this  capability  to  define 
component  test  heritage  and  response  characteristics.  Numerous  dynamic  modeling 
investigations  have  also  been  conducted  ranging  from  quick-look  failure  evaluations 
to  extensive  computer  simulations  of  system  behavior.  As  an  example  of  these 
efforts,  we  summarize  a recent  evaluation  of  the  SSME  Flight  Accelerometer  Safety 
Cut-Off  System  (FASCOS)  logic.  Some  additional  investigations  are  described  in  the 
Appendix. 


Certification  of  the  SSME  at  full  power  level  (109%)  is  in  progress.  Implementation 
of  the  FASCOS  system  is  an  essential  element  of  the  certification  program.  The 
FASCOS  is  a vibration  based  engine  condition  monitor.  The  system  continuously 
monitors  accelerations  measured  on  the  SSME  turbopumps  during  engine  operation. 
Based  on  these  measurements  and  assigned  vibration  "redlines,"  the  FASCOS  and 
engine  control  computer  determine  the  "operational  health"  of  each  engine  through  a 
series  of  nonlinear  operations  on  the  vibration  signals  and  a voting  protocol.  The 
outcome  of  the  decision  process  (run/cut-off)  is  highly  dependent  on  the  redline 
values  assigned  to  the  FASCOS  comparator  and  the  (statistical)  behavior  of  the 
monitored  signals.  For  this  reason  an  intensive  analysis/simulation  effort  was 
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recently  initiated  by  NBFC  to  define  the  FASCOS  decision  process  as  a guide  to  the 
assessment  of  redlines  with  respect  to  mission  success  and  safety.  Wyle  performed 
analytical  modeling  and  computer  simulations  in  support  of  this  challenging  and 
significant  task. 

A functional  schematic  of  the  FASCOS  logic  is  illustrated  in  Figure  12.  Three 
accelerometer  signals  are  processed  simultaneously  to  perform  the  required  turbopump 
operational  condition  assessment.  Major  elements  of  the  system  are 

• A filter/detection  stage 

• A comparator  with  time  delay 

• A voting  protocol  based  on  the  simultaneous  state  of  the  three  compara- 
tor outputs. 

The  filter /detection  stage  is  shown  in  more  detail  in  Figure  13.  This  circuit 
consists  of: 

• Band-pass  filter  (approximately  50  Hz-1  KHz) 

• Square-law  detector 

• Exponential  low-pass  filter 

• Square-root  detector 

• Second  low-pass  filter. 

It  is  seen  that  this  nonlinear  sequence  of  operations  transforms  the  (possibly) 
wide-band  vibration  input  into  a slowly  varying  "weighted  RMS"  signal  output  which 
is  fed  to  the  comparator. 

The  comparator  is  essentially  a voltage  gate  with  an  adjustable  threshold  level  and 
time  delay  required  to  activate  the  "switch."  The  threshold  level  defines  the 
vibration  redline  which  must  be  exceeded  to  initiate  an  output.  The  time  delay 
precludes  a threshold  exceedance  indication  due  to  spurious  noise  sources.  Transfer 
characteristics  of  the  FASCOS  comparator  are  illustrated  in  Figure  14.  The  output 
is  seen  typical  of  a random  telegraph  signal.  The  length  of  each  pulse  represents 
the  (random)  persistence  time  (minus  delay)  of  the  input  signal  above  the  redline 
threshold  level. 
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The  controller  periodically  samples  the  state  of  the  three  comparator  outputs  and 
determines  whether  engine  shut-down  should  be  initiated  based  on  a counting  proce- 
dure. The  voting  protocol  is  quite  simple.  On  a given  sampling  cycle  a cut-off 
vote  is  registered  only  if  all  three  controller  inputs  simultaneously  indicate  a 
redline  exceedance.  Shut-down  is  initiated  only  after  three  such  consecutive 
cycles  have  been  observed;  otherwise  the  voting  process  is  reinitialized.  An 
effort  to  demonstrate  this  voting  process  is  illustrated  in  Figure  15,  which 
indicates  a sequence  of  events  leading  to  a cut-off  decision. 

The  above  qualitative  description  of  the  FASCOS  system  indicates  a number  of  key 
issues  to  be  addressed  in  system  modeling,  simulation  and  performance  assessment, 
including 

• Statistical  characterization  of  the  class  of  (random)  functions  repre- 
sentative of  operational  vibration  signals. 

• Modeling  of  the  nonlinear  detector  output  to  representative  stochastic 
inputs. 

• Derivation  of  threshold  exceedance  persistence  time  probabilities  repre- 
sentative of  the  comparator  operation. 

• Estimation  of  joint  statistical  properties  required  to  model  the  probabil- 
ity of  combined  level  exceedance  as  utilized  in  the  controller  voting 
process. 

The  above  list  is  not  exhaustive,  but  it  does  underline  some  significant  technical 
considerations.  At  the  outset,  the  (band-passed)  vibration  signals  were  considered 
representative  of  a random  sine  wave  plus  independent  Gaussian  noise  process.  This 
model  admits  a wide  range  of  spectrum  shapes  and  signal-to-noise  level,  and  is  yet 
mathematically  tractable.  Under  this  assumption  a mathematical  description  of  the 
filter/detection  stage  output  (power  spectrum  and  covariance  function)  was  derived 
as  a function  of  input,  filter,  and  averaging  characteristics.  Limited  computer 
results  have  demonstrated  excellent  agreement  with  empirical  simulations  performed 
by  MSFC.  The  definition  of  threshold  exceedance  time  persistance  probabilities,  to 
model  the  comparator  output,  represents  a formidable  theoretical  problem.  Several 
approximate  formulations,  however,  have  been  implemented  for  evaluation.  These  are 
called  the  Rice  method,  the  method  of  non-correlated  pulses  and  the  quadratic 
approximation  method.  (See  "Non-Linear  Transformations  of  Stochastic  Processes," 
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by  Tikhonov  and  Amiartov  [RSIC  No,  QA273],  for  a comparison  of  these  techniques  and 
enabling  assumptions.)  Persistence  time  probability  density  estimates  obtained  by 
these  approximations  are  illustrated  in  Figure  16,  for  an  idealized  white  noise  plus 
sine  wave  input.  The  curves  represent  the  persistence  time  distribution  of  the 

detected  mean-square  above  a three  sigma  threshold  from  the  average  output.  Basic 
assumptions  employed  in  all  three  approximations  rentier  “*^t!iem  most  valid  for  high 
threshold-to-mean  values  and  short  time  periods.  A theoretical  extension  of  the 
Rice  method  was  developed  to  relax  these  constraints.  With  regards  to  the  final 
issue  above,  the  estimation  of  joint  properties  between  vibration  signals  appears 
best  approached  through  empirical  simulation.  In  the  absence  of  such  information, 
however,  bounds  on  system  behavior  may  be  obtained  by  assuming  joint  properties 
(total  statistical  dependence/independence  between  monitored  signals,  for  example). 

By  the  above  approach,  a computer  program  was  developed  to  define  the  probability 
of  engine  cut-off  as  a function  of  assigned  redline  values,  turbopump  vibration 
characteristics  and  FASCOS  logic  operation.  Computer  analyses  were  performed  with 
the  above  model,  based  on  vibration  characteristics  provided  by  MSFC,  to  estimate 
engine  cut-off  probabilities  as  a function  of  redline  values.  Finally,  the  program 
was  installed  on  the  MASSCOMP  processor  at  the  MSFC  Systems  Dynamics  Laboratory. 
SDL  personnel  were  provided  computer  program  orientation  at  the  conclusion  of  this 
task  to  permit  FASCOS  parameter  evaluations  by  MSFC.  A program  listing  of  this 
analytical  model  is  included  at  the  end  of  this  section. 
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Figure  16.  Three  Approximations  of  the  Persistence  Time  Distribution 


2.4  Concluding  Remarks 

This  report  provides  a cursory  overview  of  study  objectives  and  approaches  applied 
by  Wyle  in  the  performance  of  Contract  NAS8-33508.  As  a test/evaluation  program, 
rigid,  long  term  task  planning  was  neither  feasible  nor  desirable.  On  the  contrary, 
most  tasks  performed  under  this  contract  were  initiated  on  an  ad  hoc  basis,  moti' 
vated  by  observed  or  suspected  SSME  component  failure  modes.  Continued  coordination 
with  the  hGFG  COTR  was  therefore  maintained  to  revise  task  priorities  based  on  SSME 
test  results.  Consistent  with  stringent  SSME  test  schedules,  evaluation  results 
were  immediately  provided  the  COTR  in  the  form  of  presentations  and  informal  data 
packages.  A detailed  chronology  of  these  evaluations  is  given  in  the  technical 
progress  reports  provided  under  this  contract.  To  illustrate  the  diversity  of  tasks 
accomplished  under  this  contract,  some  topics  recently  investigated  by  Wyle  in 
support  of  SSME  dynamic  evaluations  are  summarized  in  Table  I at  the  end  of  this 
section.  Several  representative  evaluations  performed  under  this  contract  are 
described  in  the  Appendices. 
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TABLE  I.  SOME  TOPICS  RECENTLY  INVESTIGATED 
IN  SUPPORT  OF  SSME  DYNAMIC  EVALUATIONS 

• Pseudoharmonic  Study.  Performed  nonlinear  coherence  analysis  to  indicate  . 
whether  "apparent"  harmonics  were  correlated  with  synchronous  response. 

• Pump  configuration  analysis.  Applied  linear  discrimination  analysis  to  pump 
measurements  to  assess  statistical  difference  between  design  configurations. 

• Trimodal  density  investigation.  Performed  analysis  and  simulations  to  deter- 
mine time  series  characteristics  yielding  three  mode  probability  densities 

observed  in  turbopump  measurements. 

• Source  location  investigation.  Developed  and  applied  analysis  to  estimate 

apparent  source  location  based  on  geometry,  wave  speed  and  crosscorrelation 
measurements. 

t FPB  oxidizer,  ASl  line  accelerometer  and  strain  gauge  evaluation. 

• MMC  inlet  manifold  elbow  strain  gauge  analysis  for  strength  calculations. 

• Developed  algorithm  for  computer  data  band  program: 

Computer  selection  of  PSD  time  slices 
Cumulative  distribution 
Amplitude  indian  belt  format 

• Analysis  of  internally  instrumented  pump: 

• Transfer  functions 

Evaluation  of  ballpass  frequencies 
Calculation  of  contact  angle  from  measured  data 

• Whirl  study.  This  required  compilation  of  a complete  history  of  the  pumps, 
by  serial  number,  response  characteristics,  and  time  of  occurrence. 

• Rubbing  model  investigation.  Statistical  moments  were  derived  for  random 
sine  wave  plus  noise  processes  with  asymmetrie  clipping  for  comparison  with 
observed  turbopump  measurements. 

• STS-8.  Investigation  of  ASI  line  malfunction  required  analysis  of  previous 
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test  data  to  determine  whether  the  phenomenon  could  be  detected  on  the  high 
pressure  pumps  during  ground  test  firings. 

• Calculation  of  ball  train  frequency  for  high-speed  ball  bearing.  An  attempt 
was  made  to  modify  frequency  equations  to  include  the  effect  of  varying  ball 
size  within  a single  bearing. 

• High  synchronous  study.  Evaluated  maximum/minimum  amplitudes  by  pump  and 
component  serial  number  to  define  statistical  changes  in  response  levels. 

• Developed  relationships  for  data  extrapolation  as  a function  of  engine  power 
level. 

• Cumulative  distribution  of  HPOT.  Fit  data  to  normal,  Rayleigh  and  gamma 
distributions.  Gamma  was  found  to  be  best  fit  to  empirical  observations. 

• Modulation  investigation.  The.  thrust  was  to  determine  the  time  functions 
that  produce  certain  side-band  frequencies  which  often  appear  in  the  data 
and  are  usually  considered  noise.  A next  step  would  be  to  evaluate  the 
physical  functions  that  produce  such  time  histories. 

Results  of  the  above  investigations  have  been  provided  MSFC  through  technical 

progress  reports,  informal  data  packages  and  presentations. 
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FASCOS  PROGRAM  LISTING 


7 Microsoft  F0RTRftlM77  V3.  £0  0£ 

LINK  FASC0S4'FFTSUS+FFTa4£ 

OUTPUT  PSD  IS  ODE-SIDE  PSD  OF  XS 
DIMENSION  BBWW(£) 

CHARACTER*! 6 FPSD5, FC0R5,  F1T05 

COMMON  WTR (4096  ) , WTI ( 4096) , X R (4096) , XI (4096) , S (4100) 

INTESER  TTI,TTO 

DOUBLE  PRECISION  NFS, NF3, NF4 

WRITE  (*,  MA>  ’ ) ’ ENTER  VARIANCE  OF  XI  (T)  AFTER  BANDPASS  FILTER’ 
READ(*,*)VVV 

WRITE(*, ’ (A) ’ ) ’ ENTER  1 FOR  RC  AVERASE  FILTER’ 

READ(*,*) IFRC 

WRITE (*,’ (A)’)’  ENTER  1 FOR  ANAL08  RC  FILTER?  £ FOR  DIGITAL  RC’ 

READ(*,  *) IFAD 

IFAD-£ 

IF  (IFAD.  NE.  DTHEN 

WRIT£(*, ’ (A) ’ ) ’ ENTER  MAX  FREQ  (HZ)  ( ENTER  SAMPLING  FREQ  )’ 
READ(*, *)FS 
END  I F 

WRITE(*, ’ (A)’)’  ENTER  # OF  DATA  IN  OUTPUT  PSD  NN  ( LE  4096)’ 
WRITE<*, ’ (A) ’ ) ’ NOTE?  THIS  IS  UPTO  FS  (SAMPLING  FREQ)’ 

WRITE (*,’ (A)’)’  OUTPUT  PSD  OF  X5  HAS  NN/£  DATA  UPTO  FS/£’ 

READ(*, *)NN 

WRIT£(*, ’ (A) ’ ) ’ ENTER  FREQUENCY  OF  SIN  WAVE’ 

READ(*,  *)FSIN 

WRITE(*,  ’ (A)’  )’  ENTER  FREQ  FI  & F£  OF  FASCOS  FILTER’ 

READ(*,  *)F1,  FS 

WRITEC*,  ’ (A)  ’ ) ’ ENTER  AVERAGE  TIME  IN  SEC  ( ==£*RC)  ’ 

READ(*,  *)T 

WRITE (*,’  (A)’)’  ENTER  SIGNAL  TO  NOISE  RATIO  (NASA  DEFINITION)  ’ 
READ(*, *)ET 
ET=^£T*ET/  ( 1-ET*ET) 

WRXTE(*, ’ (A)’ ) ’ ENTER  CUTOFF  FREQ  IN  POST-RIPPLE  FILTER’ 

READ (*^  *) FC 

WRITE (*,’ (A)’)’  ENTER  ORDER  OF  BUTTERWORTH  FILTER’ 

READ(*, *)NQRDER 

WRITE(*, ’ (A) ’ ) ’ ENTER  1 FO  ANALOG?  £ FOR  DIGITAL  BUTTERWORTH’ 
READ(*, *) IFBW 

WRIT£(*, ’ (A) ’ ) ’ ENTER  MAX  TIME  LAG  OF  OUTPUT  R(T)  IN  SEC’ 

READ (*, *) TMAX 

WRITE (*,’ (A)’)’  ENTER  # OF  DATA  OF  OUTPUT  R(T)  NOUT’ 

WRIT£(*,  ’ (A)  ’ ) ’ NOTE:  NN/£-hNOUT  MUST  LE.  4096’ 

READ(*,  *)NOUT 

WRITE(*, ’ (A)’)’  ENTER  OUTPUT  FILNAM  OF  PSD  OF  X5’ 


ENTER  1 FOR  ANALOG  RC  FILTER?  £ FOR  DIGITAL  RC’ 


( ENTER  SAMPLING  FREQ  )’ 


# OF  DATA  IN  OUTPUT  PSD  NN  ( LE  4096)’ 
THIS  IS  UPTO  FS  (SAMPLING  FREQ)’ 

OUTPUT  PSD  OF  X5  HAS  NN/£  DATA  UPTO  FS/S’ 


ENTER  FREQUENCY  OF  SIN  WAVE’ 


ENTER  FREQ  FI  & F£  OF  FASCOS  FILTER’ 


ENTER  AVERAGE  TIME  IN  SEC  ( ==£*RC)  ’ 


ENTER  SIGNAL 


NOISE  RATIO  (NASA  DEFINITION) ’ 


CUTOFF  FREQ  IN  POST-RIPPLE  FILTER’ 


ENTER  ORDER  OF  BUTTERWORTH  FILTER’ 


FOR  DIGITAL  BUTTERWORTH’ 


ENTER  MAX  TIME  LAG  OF  OUTPUT  R(T)  IN  SEC’ 


ENTER  # OF  DATA  OF  OUTPUT  R(T)NOUT’ 


NOTE:  NN/£-hNOUT  MUS*! 


4096’ 


ENTER  OUTPUT  FILNAM  OF  PSD  OF  X5’ 


READ  C*,  ’ (A)  ’ ) FPSD5 


WRITE(*, ’ (A) ’ ) ’ ENTER  OUTPUT  FILNA 
READ(*, ’ (A) ’ >FC0R5 

WRITE (*, ’ (A) ’ ) ’ ENTER  OUTPUT  FILNA 
R£AD(*, ’ (A) ’ )F1T05 

WRITE(*, ’ (A) ’ ) ’ ENTER  1 FOR  SETTING  S(0)= 
R£AD(*, *) IFSp 

OPEN  (£,  FILE=*FPSD5,  STATUS==:’  NEW’  ) 

OPEN  {&,  FIL£==:F1TQ5,  STATUS==’  NEW’  ) 
********************************** 

PI  ==3.  14159H653 
NIN=NN/£ 

IFRC«1 
Fjv|AX=«FS/£. 

DT=:l/FS 
RC=T /£. 


FILNAM  OF  R(T),  R"  (T)  8c  R'  (T) 


NAM  AVE  8c  DEV  OF 


TO  X5’ 
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PRC=:£XP<-DT/RC) 

AlRC«(i-lPRC)**Jt2 
P2RC=:l+ARC*ftRG 
C2=:2*PI*DT 
C1=2*PI*RC 
BW=F2-F1 
DF=:FS/  (NN-1) 

NNMAX=2*F2/DF 
FF1=2*F1 
FF2=='BW 
FF3==BW4-FF1 
FF4=FF34.FF2 
NF1=(FF1/DF>  •+•! 

NF2«FF2/DF4.1 
NF3=FF3/DF+1 
MF4==FF4/DF+i 
ft=:VVV*0.  5/  ( BW*  < i+ET  ) ) 

VAR1=2*BW-«'A 
A2=A*A 
T2A2==^2*A2 
T4A2=4*A2 
Si==-i*T4A2 
S2=Sl/2. 

S3=-S2 
S4=S2 

B1=:sT4A2*BW 
B2«T4A2*  < BW-2*F1 ) 

B3==^B2/2. 

B4=Bl/2. 

NL2=NF2 
NL3=NF3 
NL4~NF4 

■«•*****************•«'*************** 

IF <NF 1. ST. NN) NF1=NN 
DO  11  1=1, NFl 

iFa.eE.NN+DeoTa  20 

S ( I > =B14-S1*  < I-l ) *DF 
CONTINUE 

IF<NF2. ST. 10000)  SOTO  20 
NL2=NF2 

DO  12  I=NF1-J-1,  NL2 
IFd.SE.  NN4*l)GGT0  20 
S< I ) =B2+S2* ( ( I-l ) *DF-2*F1 ) 

CONTINUE 

IF<NF3. ST. 10000)  GOTO  20 
NL3=NF3 

DO  13  I=NL2+1,NL3 
IF (1. GE. NN+l ) GOTO  £0 
3(1) =B3+S3*  < (I - 1 ) *DF"BW ) 

CONTINUE 

IF (NF4. ST. 10000)  SOTO  20 
NL4=NF4 

DO  14  I=NL34.1,NL4 
IF(I.GE.NN*4-i)GGT0  20 
S ( I ) ==B4‘+*34*  < ( I — 1 ) *DF—EiW— 2*F  1 ) 

CCNTINUE 

DO  15  I=NL4“f'l,NN 
I F ( I . GE . NN*^  1 ) GOT  □ £0 
3(1) =0. 


Woe  is 
**  qualitv 
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D L i ne#  1 
1 119  15 

120  20 
121 

1 122 
1 123 

1 124 

1 125 

1 126 
1 127 

1 128 
1 129 

1 130 

1 131 

1 132 

1 133  21 

134 

135 

136  C 

137 

138 

139 
1 140 

1 141  478 

142 

143 

144 

145 

146 

147 
1 148 

1 149 

1 150 

1 151 

1 152 

1 153 

1 154  222 

155 

156 
1 157 

1 158 

1 159 

1 160 

1 161  £2 
162 

163 

164 

165 

166 

167 

168 
169 
17-0 

171 

172 

173 

174 

175 

176 

177 
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CONTINUE 
CONTINUE 
DO  21  1=1, NN 
IF<I. BT. NL4)B0T0  21 
W= < 1-1.  )#DF 
IF  <W.  LT.  (F£rrFSINr)lT-14iN 
S (I ) =S  it  > +8*BW*ET*FI2'«C 
ELSE 

IF(W.LT.  <FSIN-F1).0R.  <W.  BT.  <F1+FSIN) , fiND.  W.  LT.  <F£+FSIN)  ) ) THEN 
S < I ) =S ( I ) +4*BW*ET*ft2 
ELSE 

S(I)»S(I> 

END  IF 
END  IF 
CONTINUE 
NC=£*FSIN/DF+1 
S (NO  = < A*BW*ET)  **2+S  (NO 
WRITEd,  *)  (S(I) , 1 = 1,  NN) 

VRV=S(1) 

I F ( I FSO.  EQ.  1 ) VRV=0. 

DO  478  1=2,  NN 
VRV=VRV+S(I) 

CONTINUE 
VRV=0RV*2*DF 
SMO=S ( 1 ) 

IF  < IFSO. EQ.  1 ) SM0«0. 

F=0 

IFdFRC.  EQ.  1)  THEN 
DO  222  1=2, NN 
F=F+DF 

IF (IFAD.  NE.  1) THEN 

S ( I ) =S ( I ) *A1 RC/ (A2RC-2»ARC*C0S (C2*F) ) 

ELSE 

S ( I ) »S  ( I ) / ( 1 + ( C 1 *F  ) **£ ) 

END  IF 

SMO=SMO+S(I) 

ELSE 

DO  22  1=2, NN 

F=F+DF 

a=PI*T*F 

S ( I ) =S ( I ) * ( ( S I N ( Q > /Q ) **£ ) 

SMO=SMO+S(I) 

CONTINUE 
END  IF 

SM0=SM0*2*DF 
AAA=2*A*BW* ( 1 +ET) 

AA4=SaRT(0. 5*SQRT (4*AAA*AAA-£*SMO) ) 

SM4=AAA-AA4*AA4 

DEV1=SQRT(WV) 

WRITE (*, 121 ) ZQ, VVV, DEVI 
DEV2=SQRT(VRV) 

WRITE (»,  1£2)WV,  VRV,  DEV2 
DEV3=SQRT(SM0) 

WRITE (*, 123) AAA,  SMO,  DEV3 
DEV4=SQRT(SM4) 

WRITE (*, 124) AA4, SM4, DEV4 
WRITE(6, *)Z0,  DEVI 
W R I TE ( 6,  * ) V VV , DE V£ 

WRITE  <6, *) AAA, DEV3 
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Line# 

1 

7 

178 

WRITE  <6, *) ftA4,  DEV4 

179 

121 

FORMAT ( » AVEl*  ’ , E16. 8, » 

VAR1  = 

180 

122 

FORMAT  < ’ AVE2*  ’ , El 6, 8, » 

VAR2=» 

181 

123 

FORMATS’  AVE3*  ’,E1S. 8, ’ 

VAR3= 

182 

124 

FORMAT <»  AVE4*  ’ , E16. S, ’ 

VAR4=» 

183 

126 

FORMAT (»  AVES*  ’ , El 6. 8, ’ 

VARS* 

184 

NFFT*NN 

185 

A4*AA4 

186 

NNP2*NFFT+e 

V 

187 

NMP22*NNP2/2 

188 

ft4S*A4*A4 

189 

A4F*A4S*A4S 

190 

PI2*2*PI 

191 

TIME*NFFT/FS 

192 

DF*FS/NFFT 

193 

DT*1/FS 

194 

DT2*DT*2 

195 

DF2*DF*2 

196 

IF<IFSO.EQ.  1)S<1)*0. 

197 

CALL  FFT (S, NNP2,  1 ) 

198 

DO  601  I*1,NNP22 

1 

199 

601 

S ( I ) *S  < 2* I - 1 ) #NFFT 

200 

DO  602  I*l,NNP2£-2 

1 

201 

602 

S < NNP22+I ) *S ( NNP22-I > 

202 

S<NFFT+1)*0. 

203 

S<NFFT+£)*0. 

204 

DO  3 I*1,NFFT 

1 

205 

S(I)*DF2*Sa> 

1 

206 

3 

CONTINUE 

207 

C 

********  CALCULATE  R4(T) 

FROM 

208 

R30=S<1) 

209 

DO  66  1=1,NFFT 

1 

210 

S < I ) *-A4S+0. 5»SQRT  <4*A4F*2*S  < I 

1. 

211 

66 

CONTINUE 

212 

R40=S(1> 

213 

125 

FORMAT  <’  R3  <0)=  ’ E16.  8,  ’ 

R4(0) 

214 

CALL  FFT  <S,  NNP2,  1 ) 

215 

DO  606  I*1,NNP22 

1 

216 

606 

S ( I ) *S  <2* I - 1 ) *NFFT 

217 

DO  607  I*l,NNP£2-£ 

1 

218 

607 

S ( NNP22+ 1 ) =S  < NNP22- 1 > 

219 

DO  4 I=1,NNMAX 

i 

220 

4 

S(I>=S(I)*DT2 

221 

S(l> =S(1) /£. 

222 

DO  224  I=NNMAX+1, NFFT 

1 

223 

224 

S < I > =0. 

£24 

C 

W R I TE  < 1 1 , * ) <S  < I ) , I « 1 , NFFT/2 ) 
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,E16.  a,  ’ DEVI*  ’ ,E16.  a) 

,E16.  a,  V DEVE*  ’ , E16.  8) 

, E16.  a,  ’ DEV3*  » , E16,  8) 

, E16.  8,  ’ DEV4*  ’ , E16.  8> 

, E16. 8, ’ DEVS*  ’ , E16. 8) 


OWGiMAL  PAGE  !S 
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R3 ( T ) ********* 


= %E16,8) 


CALL  POSTR  <S,  NFFT/a, FC,  NORDER, 
DEV5«SQRT<VflR5) 

WRITE  <*, 126) PA4, VPR5, DEVS 
WRITE <6, *) AP4, DEV5 
WRITE<£,*)  <S(I)  , I«i,NIN) 
C*********************** 

FSEQ«i.O 

SISO-1 

S 1 60«=FSEQ*  ALQG  < S I GO ) 

DLTSIG=1 

DLTSIG=FSEa*ftLOG  < DLTSI 8) 
0ME0==0- 


******  PSD  OF  X4  ********** 
DF, VAR5, FS, I FEW) 
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DLTOMG=FMPX*TMfiX/ <NIN*N0UT) 

DF=«FMAX  / <NIN-i  ) 

dt=tmax/nout 

DF£*=DF*£ 

PI =4. 0*ATAN<1. O) 

PI£=PI*£. 

DFPI£=DF*PI£ 

C******************************«**** 

C WRITE<*, ’ (A) ’ ) ’ ENTER  STARTING  FREQ  IN  HZ’ 

C R£AD(^^,  *)  0ME0 

C WRITE(*, ’ (A)’)’  ENTER  ENDING  FREQ  IN  HZ’ 

C'  READ <^,*> FEND  . * 

C DLTOMG=  < FEND-OMEO  > / < NOUT- 1 , ) 

OPEN  (12^  FILE«Ft:bR5,  STATUS^’  NEW’  ) 

NFFT=NIN-HNOyT. 

DO  30  1 = 1,  12 
NTEST=£*^«'I 

IF  <NT£ST. GE*  NFFT) GOTO  40 
30  CONTINUE 

WRITE<*,’  <A)’ )’  N TOO  BIG  FOR  FFT’ 

STOP 

40  NFFT=NTEST 

C************************************** 

IFdFSO.  EQ.  1)S(1)=0. 

DO  82  IL=1,NIN 
S ( ID  =S  (IL)/£. 

XR(ID=S<IL> 

8£  CDHTmUE 

C WRITEC*,’  <A>’ )’  BEFORE  R(T)’ 

C*******^«‘**********’t'*******  R ( T ) **»*#**************^*J«-***^****3f4«5* 

CALL  ZERO <XI, NFFT) 

CALL  CZT (XR, XI, NIN, NOUT, DLTSIG, DLTOMG, WTR, WTI, SI GO, OMEO, 

1 0,NFFT,FS£a) 

C WRITE<*, ’ (A) ’ ) ’ FINISH  R(T)’ 

DO  51  1=1, NOUT 
51  XR<I)=XR(I)*DF2 

EEE=XR(NaUT  ) *1.  0£ 

DO  301  1=1, NOUT 

IF (XR (I ) . LE. SEE) GOTO  302 

301  CONTINUE 

302  K=0 

BBB=10. / (NOUT-I ) 

DO  303  J=I,NOUT 
K=K-f-l 

XR  < J) =XR  < J) *EXP  <-BBB*K) 

303  CONTINUE 

WRITE (12, *) (DBLE(XR(I) ) , 1=1, NOUT) 

C***********-J«‘**-»'^*^**‘i«'****^  R**  (T>  ******->f**************Jt********-x*^ 
FF=-DFPI£ 

DO  502  1 = 1,  NIN 
FF=FF+DFPI£ 

502  X R (I ) =S ( I ) *FF^FF 

DO  50££  I=NIN4.1,NFFT 
50££  XR ( I ) =0, 

CALL  ZEROXXI, NFFT) 

CALL  CZTXXR,  XI, NIN,  NOUT,  DLTSIG,  DLTOMG,  WTR,  WTI, SIGO, DMEO, 

1 0,NFFT, FSEQ) 

DO  5£  I=i,NOUT 
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5£  XRa)=-DF£*XR<I) 

WRITE<i£, *) <DBLE(XR<I) ) , I=1,N0UT) 

RTT0«--XR<1) 

WRITE (*,  1£7) RTTO 
WRITECB,*) RTTO 

1£7  FORMAK’  VARIANCE  OF  SLOPE  OF  RNS  TIME  eERIES-*-R”  <0)  ^ , E16. 8) 

WRITE  (A)M’  N*  « EXPECT  # OF  UPCROSSING  OF  R IN  TIME  T’ 

WRITE (A)»)’  TBAR«  EXPECT  PERSISTENCE  TIME  OF  REDLINE  R’ 

U9^SQRT<RTJ0) /iPia^DEVm 

UPI=1. /UP 

WRITE <*, 1£8)UP 

WRITE!*, 1£9>UPI 

1£8  FORMATC’  N*  , F15.8,^»  * T * EXP  <-0.  5*R**£)  ’ ) 

1£9  FORMAT!’  TBAR  =’,F15. 8,’  * EXP <0. 5*R**£)  *ERR(R)’) 

WRITE!*,  ’ <A)  ’ )’  £RR(R>==  Pr  ( X SREATER  THEN  R ) X IS  UNIT  NORMAL’ 

C WRITE  !*,  ’ !A) ’)’  FINISH  R‘*  (T)  ’ 

C**************************  R ’ ( T > ****************************** 
FF^-DFPra 
DO  503  I«1,NIN 
FF=FF4-DFPI£ 

503  XR!I)-S(I>*FF  v 

DO  5033  I=:=NIN4-1VNFFT 

5033  XRa)=0. 

CALL  ZERO(XI,NFFT> 

CALL  CZT! XR, XI , NIN, NOUT, DLTSIG,  OLTOMS,  WTR,  WTI , SIEO,  OMEO, 

1 0,NFFT,FSEQ) 

DO  53  I«1,N0UT 

53  XR(I)=DF£*XI (I) 

WR I TE ( 1 £ , * ) < DBLE  < X R < I ) ) , I = 1 , NOUT  > 

C WRITE!*,  ’ <A) ’)’  FINISH  RV(T)  ’ 

C************************  R **  ” ( 0 ) ******************************* 

R40==0. 

FF«-DFPI£ 

DO  504  I==i,NIN 
FF=FF-4-DFPI£ 

XR<I)«=S(I)*FF**4 

504  R40«:R40+XR(I) 


333  R40«R40*DF£ 

334  WR I TE !!£,*>  DBLE  < R40 ) 

335  C***************************************************** 

336  CLOSE !£) 

337  STOP 

338  END 


Name 

Type 

Offset 

P Class 

A 

REAL 

£18 

A IRC 

REAL 

146 

A£ 

REAL 

££6 

A£RC 

REAL 

150 

A4 

REAL 

636 

A4F 

REAL 

65£ 

A4S 

REAL 

648 

AA4 

REAL 

350 

AAA 

hLQG 

REAL 

346 

INTRINSIC 

ARC 

ATAN 

REAL 

14£ 

INTRINSIC 

B1 

REAL 

£54 

11-35 


n Line#  1 
1 C 

a 

*3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

/13 

14 

15 

16 
17 
IS 
IS 
20 
21 


24 

25 

26 

27 

28 

29 

30 

31 


34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


53  501 

54  502 

55  503 

56  504 

57  505 

58  506 

59  507 
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PROGRAM  PSTAR-FOR^ — 3 STEP  MEMORY  MARKOV  MODELLING  OF  P.  T« 


IMPLICIT  REAL*8  <A-H, 0-Z) 

REAL*8  AO (50) , A1 (50) , VO (300) , PVO (300) , PSTAR (20) , QSTAR (20) 

REALMS  SM  <3, 3) , TM (3) , RM (3) ^ SMI <3, 3) , TP (3) , REDV ( 10) 

DIMENSION  L3(3),M3(3) 

REALMS  MU (20)  , DEL (20,  20) , ERR ( 1005) 

CHARACTER*16  FI , F2, FRED, FOOT 

COMMON  PI 2, ERR 

KK=3 

WRITEC*, ’ (A)’ )’  ENTER  # OF  REDLIN’ 

READ(*, *)NRED 
IF(NRED.£Q.  DTHEN 

WRITE(^«^,  MA)  ’ )»  ENTER  REDLINE  ( ? DEV  )» 

READ(*, *)REDV(1) 

ELSE 

WRITE(*,MA)’  )’  ENTER  FILNAM  OF  REDLIN’ 

READ(*, ’ (A)’ )FRED 

OPEN  <1 1 , FIL£=:FR£D,  STATUS^:’  OLD’  ) 

READdl,  *)  (REDV(I)  , i«l,  NRED) 

CLOSE (11) 

END  IF 

WRITE  (*,’  (A)  D’  ENTER  INPUT  R(T>,RMT)  FILNAM’ 

READ  (*, ’ (A)  ’ )F1 

WRIT£(*, ’ (A)’)’  ENTER  OUTPUT  INFOMATION  FILNAM  ’ 

READ(*, ’ (A) ’ )F2 

WRITEC*, ’ (A)’)’  ENTER  OUTPUT  FILNAM  OF  PSTAR’ 

R£AD(4^,  ’ (A)  ’ ) FOUT 

WRITE(#, ’ (A)’)*  ENTER  ORDER  OF  RICE  THEORY^NN  (NOT  INCLUDE  XO) ’ 
READ(*,*)NN 

WRITE (*,’ (A)’)’  ENTER  # OF  DATA  IN  ERROR  FUNCTION’ 

READ(*«^,  *)  NERR 

WRITE(*, ’ (A)’)’  ENTER  # OF  VO’ 

READ  ( *,  *)NVO 

WRITE(*, ’ (A)’ )’  ENTER  MAX  VALUE  OF  VO  ( ? DEV (VO)  )’ 

READ(*, *) VOMAX 

WRITE!*,  ’ (A)’ )’  ENTER  # OF  X’ 

READ (*, *) NX 
NX2«NX*2 


WRITE(*, ’ (A) ’ ) ’ ENTER  MAX  VALUE  OF  X ( ? DEV(X)  )’ 

READ(*, *>XMAX 

WRITE!*,’ (A)’)’  ENTER  1 FOR  PRINTING  INNER  LOOP  PARAMETERS’ 
READ (*,*) INNER 
CALL  SERRXERR,  NERR) 

OPEN  (6,  FIL£=F2,  STATUS==’  NEW’  ) 

WRITE(6, ’ (A)’)’  PROGRAM  JEN3’ 

WRITE (6,  50i)  KK 

WRITE(6, 502) NERR 

WRITE (6,  503) NN 

WRITE <6,  504) NVO 

WRITECS, 505) VOMAX 

WRITE (6,  506) NX 

WRITE  (6,  507) XMAX 

FORMAT  (’  MEMORY  STEP  =:  ’ ,12) 

FORMAT (’  # OF  DATA  IN  ERROR  FUNCTION  = ’.IS) 

FORMAT(’  ORDER  OF  RICE  THEORY  = ’,13) 

FORMAK’  # OF  INITIAL  SLOPE  VO  = ’,15) 

FORMAT(’  MAX  VALUE  OF  VO  = ’ , FS. 4, ’ DEV < VO)’) 

FORMAT  (’  # OF  X , 14) 

FORMAT<  ’ MAX  VALUE  OF  X ==  ’ , FS.  4,  ’ DEV  ( X ) ’ ) 


1 


ro  ro  u 
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D 

Line# 

1 

7 microBoft  F0RTRRN77  03. £0  0£ 

60 

OPEN  < 1 , FILE^Fl , STATUS-’  OLD’  ) 

61 

□PEN (£, FILE«F0UT, STATUS®’  NEW’  ) 

6£ 

READ<1,*) <AG<I), I®0,NN) 

! 

63 

READd,  *)  <A1  (I),  I®0,  NN) 

i 

1 

64 

65 

PI*4.D0*DftTflN  <l  . DO)  PAGE  IS 

t 

1 

66 

Pi£=Pi*£  POOR  QUALITY 

i 

67 

DVO«VOMAX*DSQRT(VftRV) /NVG 

1, 

68 

ftREfl*0. 

63 

DO  1 I = i,IMV0 

1 

70 

V0(I>==I*DV0 

, i 

71 

PVO ( I ) =DEX  P t -0.  5*V0  < I ) »V0  < I > / VAR V > *V0  < I ) 

1 

7£ 

AREA=AREA+PVO(I) 

^ 1 

73 

1 

CONTINUE 

74 

DO  £ 1=1, NVG 

i 

75 

a 

PVO ( I ) =PV0 (I ) /AREA 

76 

WRITE(*,*> (PV0<I), 1=1, NVO) 

77 

DO  3 1=1, NN 

1 

78 

DO  3 J=l  ,NN 

a 

73 

DEL  <1, J) =A0 (lABS ( I-J) ) -AO  < I ) *A0  < J) -A1 < I > *A1 (J) /VARV 

a 

SO 

3 

CONTINUE 

' 

81 

C 

sa 

DO  1000  N0PQ=1,NR£D 

1 

S3 

CALL  DZERO ( PSTAR, NN) 

1 

84 

CALL  DZERO (QSTAR, NN) 

1 

85 

RED=REDV(NOPQ) 

1 

as 

WRITERS,*)’  ’ 

1 

87 

WRITECS,*)’  » 

i 

j. 

88 

IF (INNER. EQ. 1 ) THEN 

^ 1 

as 

WRITE(6,  ■«■)  ’ ’ 

1 

30 

WRITE(S,  *)»  . » 

1 

31 

WRITE  (6,  *)  ’ » 

1 

3£ 

WRITE(6,  »)>  ’ 

1 

33 

WRITERS,*)’  ’ 

1 

34 

WRITE  (6,*)  ’ ’ 

i 

35 

ENDIF 

. 1 

36 

WRITE (6, 606) RED 

1 

37 

WRITE<*,606)RED 

1 

98 

606 

FORMAT  (IX,  ’ ***■»**#*#»#•»*»•»***  REDLINE=’ , F8.  4,  ’ **********»***»  ) 

1 

33 

1 

100 

DO  9999  MM=KK+1, NN 

a 

101 

WRITE  (*,  ’ (A)  ’ ) ’ **■*■•»•»♦•*•♦•»■»***********•»•«•**•»■■*•******«*•»#*»***•«•*»*» 

; C. 

loa 

IF  ( INNER.  EQ.  1 ) WRITE  (6,  ’ (A)  ’ ) ’ #*****»*****#«-**»**»»*»*******«-*> 

. a 

103 

BM=DEL(MM, MM) 

i a 

104 

DO  9991  1=1, KK 

3 

105 

TM(I)=DEL(MM, MM-I) 

' 

106 

DO  9991  J=1,KK 

4 

107 

3331 

SM ( I , J ) =DEL ( MM- I , MM- J ) 

' a 

108 

CALL  EQUALS (SM, SMI , KK) 

1 o 

tm. 

103 

CALL  DMINV(SMI, KK, DTM, L3, M3) 

1 c. 

110 

CALL  DSMPRD (TM, SMI , TP, 1 , KK, KK) 

a 

111 

VZM=BM-TP ( 1 ) *TM ( 1 ) -TP ( £ ) *TM ( £ ) -TP ( 3 ) *TM ( 3 ) 

a 

iia 

DZM=DSQRT(VZM) 

iZ. 

113 

DEVXl=DSQRt(SM(l,  1) ) 

a 

114 

DEVX£=DSQRT (SM  <£,£)) 

115 

DEVX3=DSQRT  (SM(3,  3) ) 

o 

d 

116 

C0R1£=SM ( 1,  £) / (DEVXi*DEVX£) 

a 

117 

CORl 3=SM ( 1 , 3) / ( DEVX 1 *D£VX3 ) 

o 

u. 

1 18 

C0R£3=SM  <£, 3) / ( DEVX£*DEVX3) 

D 

Line# 

1 

a 

119 

£ 

lao 

£ 

lai 

£ 

laa 

£ 

1£3 

£ 

1£4 

a 

las 

a 

las 

a 

1£7 

601 

a 

las 

a 

1£9 

3 

130 

4 

131 

991 

3 

* 13£ 

3 

133 

4 

134 

99£ 

3 

135 

3 

136 

3 

137 

3 

138 

3 

139 

3 

140 

3 

141 

3 

14  a 

3 

143 

3 

144 

3 

145 

3 

146 

3 

147 

3 

148 

3 

149 

3 

150 

3 X51 


4 

15£ 

4 

153 

4 

154 

4 

155 

4 

156 

4 

157 

4 

158 

4 

159 

4 

160 

5 

161 

5 

16£ 

5 

163 

5 

164 

5 

165 

5 

166 

5 

167 

5 

168 

6 

169 

6 

170 

6 

171 

6 

17£ 

6 

173 

6 

174 

6 

175 

6 

176 

6 

177 
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C0RR=1-C0R»C0R 
DX 1=XMAX*DEVX 1 /NX 
D X 2=  X Mft  X *DEV  X2 /NX 
DX3«XMftX«-DEVX3/NX 

Q0=DX1*DX2*DX3/  < PI2#*1.  5*DSC3RT  <DABS  <DTM)  ) > 

IF  aNNER%-EQl'l ) WRITE  (S,  601)  DZM,  DEVXl,  DEVX2,  DEVX3,  CORIE,  CDR23, 

1 C0R13 

WRITE (*,  601 ) DZM,  DEVX 1 , DEVX2,  DEVX3,  C0R12, C0R23,  C0R13 
FORMAT (IX, 7 <E9,  4,  IX) ) 

WRITEC*,  ’ <A)» )’  BEFORE  999’ 

DO  999  1I-1,NV0 
DO  991  1*1, NN 

MU ( I ) *RED*AO ( I ) -VO ( 1 1 ) *A 1 ( I ) / VA  R V 
AM*MU(MM) 

DO  992  1*1, KK 
RM<I)*MU(MM-I) 

XMIN1=RM  < 1 ) -XMAX*DEVX1 
XMIN2=RM  (2)  ■*XMAX#DEVX2 
XM I N3»RM ( 3 ) - X MAX  *DE V X 3 
XMAX 1*RM < 1 ) +XMAX*DEVX 1 
XMAX2*RM (2) +XMAX*DEVX£ 

XMAX3*RM (3) +XMAX*DEVX3 
Sl=DMAXl (RED, XMINl) 

S£=DMAX1 (RED,  XMIN2) 

S3=DM A X 1 ( RED , XM I N3 ) 

NS1=(S1~XMIN1) /DXl 
NS2= (S£-XMIN£> /DX2 
NS3* (S3-XMIN3) /DX3 
X1*S1-DX1 
PS*0. DO 
QS*0. DO 

IF(S1.  GT.  XMAXDOOTO  100 
DO  99  JJ«NS1,NX2 
X1*X1+DX1 
Q1*(X1-RM(1) ) 

H1*G!1*Q1*SM  I (1,  1) 

WK2=2*ai*SMI  (1,2) 

WK3=Q1*SMI (1,3) 

AZM1=AM4-TP(1)  *Q1 

IF (S2. GT. XMAX2) GOTO  99 

X2=S£-DX£ 

DO  9 LL=NS2,NX2 
X2*X2+DX2 
Q2=X£-RM(2) 

H2*H1+Q£» (Q£*SMI (2, 2) +WK2) 

WK1=2*  (WK3+Q2*SMI (2,  3)  ) 

AZM£=AZM1+TP (2) *Q£ 

IF (S3. GT. XMAX2) GOTO  9 
X3=S3-DX3 
DO  22  IJK*NS3, NX2 
X3=X3+DX3 
Q3*X3-RM(3) 

H3*H2+  (WK1+Q3*SMI (3,  3) ) *Q3 
AZM*AZM£+TP (3) *03 
GGB=GG1 < (RED-AZM) /DZM, NERR) 
aa=DEXP (-H3/2. DO) 

PS=PS+-aO*GSG 

QS*QS"J"QQ 

Gontinue 
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0 

Line# 

1 

7 * Microsoft 

F0RTRPN77 

V3. 

£0  0£/84 

5 

178 

9 

CONTINUE 

4 

179 

99 

CONTINUE 

3 

ISO 

100 

PS=P9*Q0*PV0(II) 

3 

181 

QS^QS*Q0*PV0(II> 

3 

182 

PSTPR  (MM)  =PSTftR  (MM)  4-PS 

3 

183 

QSTPR(MM) =QSTPR (MM) ^QB 

3 

184 

IF  < INNER. EQ. 1 ) WRITE(6, 60S) 1 1, PVO  <I I ) , PS,  QS,  NSl,  NSS,  NS3, 

185 

1 <RM<I> , I-l, KK) 

3 

186 

WRITE  <*, 60S) 1 1 , POO  < 1 1 ) , PS,  QS,  NSl , NSS,  NS3,  ( RM  < I ) , 1 = 1 , KK) 

3 

187 

60£ 

FORMAT ( IX, 14, 3<1X, ES. 3) , 314, 3 (IX, E8.  3) ) 

3 

188 

999 

CONTINUE 

£ 

189 

PQ=PSTAR(MM) /QSTAR<MM> 

£ 

190 

IF  < MM. EQ. KK+ 1 ) THEN 

£ 

191 

PP1=PSTAR  <MM) 

£ 

19£ 

ELSE 

£ 

193 

ppi=ppi*pa 

£ 

194 

ENDIF 

£ 

195 

IF ( INNER. EQ. 1) THEN 

£ 

196 

WRITE  <6,  ’ <A)  ’ ) » #***  MM  P*  Q#  *■*■** 

P*/Q*  **** 

PP 

2 

197 

ENDIF 

£ 

198 

WRITE(6, 603)MM, PSTAR(MM)  , QSTAR(MM) , PQ,  PPl 

£ 

199 

WRITE (*,’ (A)’)’****  MM  ****  P*  **»*  Q*  **** 

P*/Q*  **** 

PP 

***** 

£ 

£00 

WRITE(*,603)MM, PSTAR (MM) , QSTAR (MM) , PQ, PPl 

£ 

£01 

603 

FORMAT ( 1 X, 15, 4E16. 6) 

o 

c. 

£0£ 

9999 

CONTINUE 

X 

£03 

C 

WRITE(S, *) (PSTAR(I), 1=1, NN) 

1 

£04 

c 

WRITE  (S,*)  (QSTARd) , 1=1,  NN) 

1 

£05 

W1=PSTAR(KK+1> 

1 

£06 

WS=1. 

1 

£07 

DO  1001  I=KK+S, NN 

fIZL 

£08 

W1=W1*PSTAR(I) 

£ 

£09 

1001 

WS=WS*QSTAR(I) 

1 

£10 

IF  (WS.  NE.  O.  DO)  PP=W1/WS 

1 

£11 

c 

WRITE(S,  *)PP 

1 

£1£ 

WRITE(6, » (A) ’ ) ’ FILNAL  PROBABILITY  IS:’ 

1 

£13 

WRITE(6,*) Wl  , WS,  PP 

A 

£14 

WRITE(*, ’ (A) ’ ) ’ FILNAL  PROBABILITY  IS:’ 

1 

£15 

WRITE  (#,  *)  Wl  , WS,  PP 

1 

' £16 

WRITERS,  *) PP 

1 

£17 

1000 

CONTINUE 

£18 

STOP 

£19 

END 

Name  Type  Offset  P Class 


PO 

REPL*8 

£ 

^ Pi 

REPL*S 

40£ 

PM 

REAL *8 

10£4B 

PREP 

REPL*8 

9972 

PZM 

REAL *8 

10496 

1 PZMl 

R£AL*S 

10408 

1 PZM£ 

REAL*8 

10456 

: BM 

REPL*8 

10084 

COR 

REPL*8 

10180 ■ 

CGR1£ 

REAL*a 

10148 

' COR  13 

R£PL*8 

10156 

C0R£3 

REAL*8 

10164 

! CORR 
j DPBS 

R£AL*a 

101 7£ 

INTRINSIC 

1 

II- 

7 


D Line#  1 7 

DATflN 

DEL 

REAL*8 

DEVXl 

REAL*8  . 

DEVX2 

REALMS 

DEVX3 

DEXP 

DMflXl 

DSQRT 

REAL*8 

DTM 

REAL*8 

DVO 

REALMS 

DXl 

REAL*8 

DX2 

REAL*8 

DX3 

REAL*8 

DZM 

REAL*8 

ERR 

REAL*8 

FI 

CHAR* 16 

F2 

CHAR*16 

FOUT  ' 

CHAR* 16 

FRED 

CHAR* 16 

GSG 

REAL*8 

HI 

REAL*8 

H2 

REAL*a 

H3 

REAL*a 

I 

IfiBS 

INTEGER*4 

II 

INTEGER*4 

UK 

INTEGER*4 

INNER 

INTEGER*4 

J 

INTEGER*4 

JJ 

INTEGER*4 

KK 

INTEGER*4 

L3 

INTEGER*4 

LL 

INTEGER*4 

M3 

INTEGER*4 

MM 

INTEGER*4 

MU 

REAL*8 

NERR 

INTEGER*4 

NN 

INTEGER*4 

NOPQ 

INTEGER*4 

NRED 

INTEGER*4 

NSl 

INTEGER*4 

NS2 

INTEGER*4 

NS3 

INTEGER*4 

NVO 

INTEGER*4 

NX 

INTEGER*4 

NX2 

INTEGER*4 

PI 

R£AL*a 

PI2 

REAL*8 

PP 

R£AL*a 

PPl 

REAL*a 

PQ 

R£AL*8 

PS 

REAL*8 

PSTAR 

REAL*8 

PVO 

REAL*8 

QO 

REAL*8 

Q1 

REAL*8 

Q2 

REAL*a 

Q3 

REAL* 8 

QQ 

R£AL*8 

INTRINSIC 

6402 

10124 

10132 

10140 

INTRINSIC 

INTRINSIC 

INTRINSIC 

10100 

9964 

loiaa 

10196 

10204 

10116 

a /COMMQQ/ 
9634 
96S0 
9666 
9610 
10504 
10394 
10440 
i04aa 
9626 

INTRINSIC 

10236 

10472 

9718 

9992 

10368 

9602 

6218 

10424 

6230 

10076 

6242 

9686 

9682 

10000 

9606 

10332 

10336 

10340 

9690 

9702 

9706 

9956 

0 /COMMQQ/ 
10614 
10570 
10562 
10352 
6058 
3202 
10212 
10376 
10432 
10480 
10520 
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1 1 -40 


7 


D Line# 

= 1 

QS 

REAL*a 

QSTftR 

R£AL*6 

RED 

REPL*8 

REDV 

R£ftL*8 

RM 

REftL*8 

SI 

REftL*a 

SE 

REAL*8 

S3 

REAL*8 

SM 

REALMS 

sm 

REAL*8 

TM 

REAL*a 

TP 

REftL*8 

VO 

REALMS 

VOl^ftX 

REftL*8 

VftRV 

R£AL*8 

VZM 

REALMS 

Wi 

REAL*a 

WE 

REAL*a 

WKl 

REAL*a 

UIKS 

REAL*8 

WK3 

R£AL*S 

XI 

REAL* 8 

X£ 

REAL*8 

X3 

REAL*8 

xmx 

R£AL*8 

XMPXi 

REAL*8 

XI^PlX£ 

REAL*8 

XMAX3 

REAL*8 

R£AL*8 

Xr^sINE 

REAL*a 

miH3 

R£AL*8 

1 0360 
560£ 
lOOOS 
5978 
5954 
10308 
10316 
10324 
5762 
5882 
5834 
5858 
802 
9694 
9948 
10108 
10594 
10602 
10448 
10392 
10400 
10344 
10416 
10464 
9710 
10284 
10292 
10300 
10260 
10268 
10276 
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ORIGINAL  PAGE  IS 
OF  POOR  QUALHY 


SUBROUTINE  EQUflL3<fi, B, N> 
REAL*8  ft<3,3),B(3,  3) 

DO  1 1 = 1,  N 
DO  1 J=1,N 
B (I,  J)  =fl(  I,  J) 

RETURN 

END 


Name 

Type 

0 f f se  t P Class 

A 

REAL*8 

0 * 

B 

REAL*a 

4 * 

X 

INTESER*4 

106££ 

J 

INT£SER*4 

10630 

N 

INTESER*4 

8 * 

1 

2 


220 

221 

£22 

£23 

£24 

225 

£26 


££7 

SUBROUTINE  S£RR<£RR 

, NERR) 

££8  C 

AREA  UNDER  Nx<0, 1) 

FRGN  X=0  TO  R 

££9 

REAL*8  £RR(1),PI,C1 

, T,  TPl,  TP£,  DT 

£30 

PI =4- DO*DATAN  < 1.  DO) 

£31 

Cl^DBQRT (£, DO*PI) 

£3£ 

DT=:5-  DO/NERR 

£33 

DT£=:DT/£. 

£34 

T==:0.  DO 

£35 

DO  1 1=^1,  NERR 

11-41 


D;  Line#  1 
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T==T4.DT 

ERR  < I ) =DEXP (-0. 5*T*T) /Cl 
TPl=:£RR<i) 

ERR  < 1 ) « (ERR  <1 ) +1 /Cl ) *DT*0. 5 

DO  2 1=^2,  NERR 

TP£=eRR<I) 

ERR ( I ) =ERR  < I - 1 ) +DT2*  ( TP 1 +ERR ( I ) ) 

TP1=TP2 
CONTINUE 
ERR<0)  =0.  DO 
RETURN 
END 
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X=DABS (XX) 

NN=X*NERR/5,  do 
IF (NN. BE, HERR) THEN 
S61=:0.  5 
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END 
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1 

£79 
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RETURN 
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END 
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£ 
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£ 
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3 

30S 

10 
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3 

307 

15 

BIGR=«(IJ> 

3 

308 

L(K)==:I 

3 

309 

M<K)=J 

3 

310 

£0 

CONTINUE 

1 

311 

J=L (K) 

1 

3 1 £ 

IF(J-K)  35,  35, £5 

4 

X 
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£5 

KI=K-N 

1 

314 

DO  30  1=^1,  N 

£ 

315 

KI=KI-i-N 

£ 
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H0LD=-A(KI) 

£ 
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JI===KI-K4-J 

£ 
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A(KI)==A<JI) 

£ 
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30 
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1 
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35 
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X 
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38 
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3£3 
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3£5 
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£ 
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o 
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45 
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1 

1 

RETURN 
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c! 
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85 
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1 
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£ 

75 

CONTINUE 

1 

354 

D=D*BIBA 

•« 

355 

A<KK)=1- 0/BI BA 

358 

SO 

CONTINUE 

Page  9 
0E-13~8S 
OSsOEs 19 

Microsoft  F0RTRAN77  V3. £0  0£/S4 


15,  £0,  £0 


11-44 


09502s  IS 

Microsoft  F0RTRftN77  V3, £0  0£/S4 


Line# 

1 

7 

357 

K=N 

358 

O 

O 

K=<K-«i) 

359 

IF<K)  150, 150,  105 

360 

105 

I==L(K) 

361 

IF(I-K)  1£0, i£0, 108 

36£ 

108 

JQ=N*  <K-i ) 

363 

JR«N*<I--1) 

364 

DO  110  J==1,  N 

1 

365 

JK*JD+J 

X 

366 

HDLD=fl  < JK) 

il 

367 

JI=JR+J 

% 

368 

F 

369 

110 

A(JI)=HOLD 

1 

370 

leo 

J~M (K) 

371 

IF(J-K)  100,100,  1 £5 

37£ 

i£5 

KI~K--N 

373 

DO  130  1 = 1,  N 

11 

374 

KI=KI+N 

375 

HOLD=A(KI ) 

376 

JI=KI-K*4-J 

>1 

377 

«(KI)=-A CJIX 

378 

130 

A<JI)=HOLD 

i 

379 

00  TO  100 

380 

150 

RETURN 

i 

1 

1 

381 

END 

fNarne 

Type 

Offset  P 

P 

REALMS 

0 * 

BIGft 

R£RL*a 

10778 

D 

R£ftL*8 

8 * 

D«BS 

HOLD 

R£ftL*B 

10830 

I 

INTE0ER*4 

1079B 

IJ 

INTEG£R*4 

10806 

IK 

INTEGER*4 

10858 

IZ 

INTEGER*4 

10794 

J 

INTEGER*4 

10786 

JI 

INT£GER*4 

10838 

JK 

INTEGER*4 

10850 

JP 

INT£GER*4 

1084S 

JQ 

INTEG£R*4 

10878 

JR 

INT£GER-J<‘4 

10882 

K 

INTEGER*4 

10766 

KI 

INT£G£R*4 

10822 

KJ 

INTEGER*4 

10870 

KK 

INTEG£R*4 

10774 

L 

INTEGER*4 

12  ^ 

INT£GER*4 

16  * 

N 

INTEGER*4 

4 ^ 

NK 

INT£e£R*4 

10762 

Name 

Type 

Size 

CGi^.^OQ 

8048 

DSNPRD 

DMINV 

Class 


INTRINSIC 


Class 

CONNON 

SUBROUTINEE 

SUBROUTINE 


H-45/46 


ORIGINAL  PAGE  IS 
OP  POOR  QUALITY 


APPENDIX  A 

WYLE  LABORATORIES  - RESEARCH  STAFF 
TECHNICAL  MEMORANDUM  TM  80-8 


EVALUATION  OF  SPACE  SHUTTLE 
MAIN  ENGINE  FLUID  DYNAMIC 
FREQUENCY  RESPONSE  CHARACTERISTICS 


APPENDIX  A 


WYLE  LABORATORIES  - RESEARCH  STAFF 
TECHNICAL  MEMORANDUM  TM  80-8 


EVALUATION  OF  SPACE  SHUTTLE 
MAIN  ENGINE  FLUID  DYNAMIC 
FREQUENCY  RESPONSE  CHARACTERISTICS 


by 

T.  G,  GARDNER 


for 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 
GEORGE  C.  MARSHALL  SPACE  FLIGHT  CENTER 
MARSHALL  SPACE  FLIGHT  CENTER,  ALABAMA  35812 


Work  Performed  Under  Contract  Number  NAS8-33508 


FOREWORD 


This  report  was  prepared  by  Wyle  Laboratories,  Research  and 
Engineering  Division,  Huntsville,  Alabama,  under  NASA  Con- 
tract No.  NAS8-33508,  for  the  National  Aeronautics  and  Space 
Administration  (NASA),  George  C.  Marshall  Space  Flight 
Center.  The  contract  was  administered  under  the  technical 
direction  of  the  Systems  Dynamics  Laboratory,  with  Mr.  Harry 
Bandgren  acting  as  the  Technical  Contracting  Officer's 
Representative.  Mr.  Duron  Cryder  was  the  contract  admini- 
strator for  NASA. 


TABLE  OF  CONTENTS 


Page 

FOREWORD  ........  . . .....  . . . . . . ..........  fi 

ABSTRACT  ........  . . . . .............  . v 

SPACE  SHUTTLE  MAIN  ENGINE  SYSTEM  DESCRIPTION  ....  . . . . . . . . . 1 

SPACE  SHUTTLE  MAIN  ENGINE  POGO  TEST  PROGRAM.  ...  .....  5 

POGO  DATA  ANALYSIS  SOFTWARE.  ..  ...........  . 7 

APPENDIX  A - POGO  SOFTWARE  FLOWCHART  ......  .....  17 

APPENDIX  B - POGO  SOFTWARE  LISTING  ...  1*1 

APPENDIX  C - POGO  SOFTWARE  TEXT  RECORDS  (FILE  4a) . ....  . . ....  57 

APPENDIX  D - POGO  SOFTWARE  TEXT  BUFFERS  (FILE  4b).  ..........  6l 

APPENDIX  E - POGO  SOFTWARE  SAMPLE  OUTPUT  67 


EVALUATION  OF  SPACE  SHUTTLE  MAIN  ENGINE 
FLUID  DYNAMIC  FREaUENCY  RESPONSE  CHARACTERISTICS 


by 

T.  G.  Gardner 
ABSTRACT 

In  order  to  determine  the  POGO  stability  characteristics  of  the  Space 
Shuttle  main  engine  (SSME)  liquid  oxygen  (LOX)  system,  an  evaluation 
of  the  fluid  dynamic  frequency  response  functions  between  elements  in 
the  SSME  LOX  system  was  performed,  both  analytically  and  experiment- 
ally. This  report  acquaints  the  reader  briefly  with  the  POGO  testing 
program  and,  more  specifically,  documents,  as  a user  note,  the  POGO 
data  analysis  software.  POGO  refers  to  the  effect  the  dynamic  inter- 
action between  devices  in  the  LOX  system  has  on  fluid/mechanical 
vibration  of  the  system.  For  the  experimental  data  evaluation,  a 
software  package  was  written  for  the  Hewlett-Packard  5^51 C Fourier 
analyzer.  The  POGO  analysis  software  consists  of  five  separate  seg- 
ments. Each  segment  is  stored  on  the  5451C  disc  as  an  individual 
program  and  performs  its  own  unique  function.  The  POGO  analysis 
software  includes  two  separate  data  reduction  methods,  a signal 
calibration,  coherence  or  pulser  signal  based  frequency  response 
function  blanking,  and  automatic  plotting  features.  The  5451C  allows 
variable  parameter  transfer  from  program  to  program.  This  feature  is 
used  to  an  advantage  and  requires  only  minimal  user  interface  during 
the  data  reduction  process.  Experimental  results  are  also  included. 
Comparison  of  experimental  results  with  the  analytical  predictions 
permits  adjustments  to  the  general  model  in  order  to  arrive  at  a 
realistic  simulation  of  the  POGO  characteristics. 
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SPACE  SHUTTLE  MAIN  ENGINE  SYSTEM  DESCRIPTION 


The  Space  Shuttle  propulsion  system  consists  of  two  solid  rocket  boosters 
(SRBs)  and  three  main  engines  (SSMEs).  The  SSMEs  are  high-performance, 
liquid-propellant,  variable  thrust  rocket  engines,  operating  at  high 
temperatures,  high  pressures,  and  high  rotational  speeds.  Each  SSME 
operates  at  a chamber  pressure  of  approximately  3000  psia  to  produce  a 
sea  level  thrust  of  375,000  pounds  and  vacuum  thrust  of  470,000  pounds 
and  operates  over  a range  from  50  to  109  percent  of  the  rated  power. 

Each  engine  consists  of  ten  major  components:  two  preburners,  four 

turbopumps,  the  hot-gas  manifold,  the  main  injector,  a heat  exchanger, 
and  the  main  combustion  chamber.  All  other  components  are  attached 
structurally  to  the  hot-gas  manifold,  and  the  entire  arrangement  is 
called  the  SSME  power head. 

A key  to  the  cost  effectiveness  of  the  Space  Shuttle  concept  is  hardware 
reusability.  As  a result,  system  reliability  is  of  paramount  importance. 
The  SSMEs  have  been  subjected,  to  extensive  hot  firing  and  flow  tests. 
Under  these  extreme  operating  conditions,  system  failures  and  malfunc- 
tions have  occurred  due  to  the  self- induced  dynamic  environment.  These 
failures  have  ranged  from  subcritical  wear  of  component  bearings  to 
catastrophic  explosion  and  fire  resulting  from  the  intense  pressure 
oscillations  and  dynamic  stresses  occurring  in  pumps,  valves,  and/or 
propellent  lines.  The  components  in  the  liquid  oxygen  (LOX)  system 
include  the  low  pressure  oxidizer  pump  (LPOP),  POGO  suppression  system 
(or  accumulator),  the  high  pressure  oxidizer  pump  (HPOP) , and  the  main 
combustion  chamber  (MCC).  The  LOX  flows  from  the  external  tank  (ET) 
through  the  above  components  and  is  combined  with  the  liquid  hydrogen 
(LH2)  fuel  at  the  inlet  of  the  MCC  by  the  main  injector. 

The  LPOP  is  an  axial-flow  pump  driven  by  a six-stage  turbine  and  powered 
by  the  LOX.  During  engine  startup  and  mainstage,  the  LPOP  maintains 
sufficient  pressure  in  the  LOX  line  to  permit  the  HPOP  to  operate  at 
high  speeds  without  an  Inducer  and  without  cavitation,  even  under 
worst-case  conditions. 
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The  POGO  suppressor  is  a gas-filled  accumulator,  which  serves  as  a 
capacitance  in  the  LOX  flow  circuit.  The  unit  is  designed  to  prevent 
low-frequency  oscillations,  transmitted  from  the  Space  Shuttle  vehicle, 
from  being  transmitted  into  the  HPOP  and  ultimately  to  prevent  MCC 
pressure  oscillations.  To  provide  suppression  protection  through  the 
startup  and  shutdown  transients,  the  accumulator  is  filled  with  gaseous 
helium.  During  normal  engine  operation,  the  heat  exchanger  provides 
gaseous  oxygen  (GO2),  as  the  compliant  medium,  to  the  accumulator.  The 
system  consists  of  a 0.6  cubic-foot  accumulator,  which  is  attached  to 
the  HPOP  inlet  duct,  an  internal  stand  pipe,  helium  precharge  valve 
package,  gaseous  oxygen  supply  valve  package,  and  recirculation  isola- 
tion valves.  The  liquid  level  in  the  accumulator  is  controlled  by  the 
stand  pipe,  which  is  orificed  to  regulate  the  GO2  overflow  over  the 
engine  operating  power  level  range.  Excess  gaseous  and  liquid  oxygen 
are  recirculated  back  to  the  LPOP  inlet  through  the  engine  oxidizer 
bleed  duct. 

The  HPOP  consists  of  two  single-stage  centrifugal  pumps  on  a common 
shaft,  which  is  directly  driven  by  a two-s.tage  hot-gas  turbine.  The 
main  pump  receives  oxidizer  from  the  LPOP  discharge  and  supplies  LOX  at 
increased  pressure  to  the  LPOP  turbine.  The  HPOP  turbine  is  powered  by 
hot-gas  (hydrogen-rich  steam).  LOX  enters  the  HPOP  main  pump  through 
the  main  pump  housing  and  flows  through  an  inlet  with  a 50-50  flow  split 
into  a double-entry,  common  outlet  impeller.  Several  sets  of  guide 
vanes  direct  the  flow  to  the  impeller  inlets.  The  impeller  has  four 
full  and  four  partial  blades  in  each  half.  After  passing  through  the 
impeller,  the  flow  is  redirected  into  the  HPOP  discharge  by  diffuser 
vanes.  The  HPOP  is  attached  to  the  hot-gas  manifold  by  a flange,  which 
is  canted  ten  degrees  out  from  the  engine  centerline.  The  oxidizer  then 
flows  into  the  MCC. 

The  MCC  is  a cylindrical  structural  chamber,  which  contains  the  burning 
propellant  gases  and  is  flange-attached  to  the  hot-gas  manifold.  The 
MCC  consists  of  a coolant  liner,  coolant  Inlet  and  outlet  manifolds,  and 
a structural  jacket.  The  coolant  liner,  which  provides  the  coolant 
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flow  path  of  the  MCC,  contains  390  milled  axial  coolant  channels,  which 
are  ported  to  the  coolant  inlet  and  outlet  manifolds.  This  network  pro- 
vides an  up-pass  circuit  for  the  liquid  hydrogen  coolant.  The  chamber 
jacket  provides  the  structural  strength  for  the  MCC  and  is  approximately 
20  inches  long.  The  jacket  is  formed  in  two  matching  halves,  which  are 
welded  together  and  to  the  coolant  inlet  and  outlet  manifolds.  A throat 
ring  is  also  welded  to  the  jacket  at  the  MCC  throat  to  provide  added 
strength.  The  MCC  throat  follows  a contraction  from  the  hot-gas  mani- 
fold of  2.96:1  and  is  expanded  to  a ratio  of  5:1  before  attaching  to  the 
engine  nozzle. 

In  order  to  validate  system  performance  and  ensure  equipment  reliability, 
the  SSME  and  components  have  been  and  are  presently  undergoing  extensive 
development  and  qualification  tests.  Testing  of  the  engine  and  components 
is  conducted  at  several  NASA  and  contractor  locations.  Full-scale 
engine  test  firings  for  development  and  flight  acceptance  are  performed 
on  two  single-engine  test  stands  at  the  National  Space  Technology  Labora- 
tories (NSTL),  Bay  St.  Louis,  Mississippi,  and  one  stand  operated  by 
Rockwell  International  near  Santa  Susana,  California,  in  addition,  main 
propulsion  testing  is  performed  at  NSTL  on  a stand  designed  to  accommo- 
date the  Shuttle  main  propulsion  system  element — the  three-engine 
cluster,  external  tank,  and  orbiter  systems. 

Testing  is  being  performed  on  a continuing  basis.  The  length  of  a given 
test  is  dependent  on  specific  test  objectives  and  may  run  from  several 
seconds  to  over  '800  seconds.  During  each  test.  Comprehensive  measure- 
ments are  acquired  to  monitor  system  performance,  including  vibration, 
dynamic  pressure  and  strain  at  critical  engine  locations.  Several  of 
these  latter  measurements  are  utilized  on-line  as  emergency  cutoff 
indicators,  and  all  are  recorded  on  magnetic  tape  for  subsequent  analy- 
sis and  evaluation. 
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SPACE  SHUTTLE  MAIN  ENGINE  POGO  TEST  PROGRAM 


As  part  of  a program  to  determine  the  POGO  stability  characteristics  of 
the  Space  Shuttle  Main  Engine  (SSME)  liquid  oxygen  (LOX)  system,  POGO 
tests  were  conducted  by  Rocketdyne  engineers,  in  California,  on  the  A3 
SSME  test  stand.  In  order  to  define  the  system  fluid  dynamic  characteris- 
tics, frequency  response  functions  [H(f)]  between  various  components  are 
calculated  from  pressure  measurements  taken  at  various  locations  in  the 
system.  To  calculate  the  H(f)s,  a known  pressure  signal  must  be  applied 
to  the  system  and  the  system  response  to  this  signal  measured.  A 
pulser  was  used  to  generate  this  dynamic  pressure  signal,  with  the  tests 
usually  including  both  sine  sweeps  and  sine  dwells.  The  first  device 
located  downstream  of  the  pulser  is  the  low  pressure  oxidizer  pump 
(LPOP) , followed  by  the  accumulator  (or  suppressor),  the  high  pressure 
oxidizer  pump  (HPOP) , and  the  main  combustion  chamber  (MCC).  With  the 
accumulator  active  in  the  system,  the  pulser  signal  is  suppressed  such 
that  it  does  not  reach  the  MCC  with  sufficient  strength  to  be  measured 
adequately.  Consequently,  tests  were  conducted  with  the  accumulartor 
active  and  inactive.  The  general  system  diagram  below  shows  the  rela- 
tive location  of  these  devices.  Pressure  measurements  were  made  at  the 
pulser,  the  LPOP  inlet,  the  HPOP  inlet,  and  the  HPOP  discharge. 


The  pulser  sine  sweeps  ranged  In  frequency  from  2 to  kO  Hz.  At  the 
higher  frequencies  in  the  pulser  sweep,  the  displacement  of  the  pulser 
piston  and  the  corresponding  amplitude  of  the  pressure  pulse  generated 
were  very  small.  The  amplitudes  of  the  pressure  signals  further  down 
the  line  were  essentially  buried  in  the  noise  level  and  very  difficult 
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to  measure.  For  the  sine  dwells,  the  magnitude  of  the  pulsed  signal  can 
be  increased,  providing  greater  dynamic  range  in  the  data.  The  problems 
associated  with  the  dwell  testing  are  the  increased  time  required  to 
conduct  the  test  and  the  fact  that  the  frequency  response  functions  can 
be  defined  only  at  the  frequencies  where  the  dwells  are  located. 

Because  of  this,  a routine  that  will  blank  out  invalid  data  in  frequency 
response  functions  was  included  in  the  POGO  analysis  software.  This 
blanking  can  be  based  on  either  the  pulser  signal  or  the  coherence 
function,  leaving  only  the  frequency  response  function  data  at  the 
frequeneies  where  the  pulser  was  operated.  In  the  initial  tests,  the 
dwell  frequencies  were  chosen  at  approximately  5~Hz  intervals.  In 
subsequent  tests,  the  dwells  were  chosen  closer  together  in  order  to 
better  define  the  system  frequency  response. 
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POGO  DATA  ANALYSIS  SOFTWARE 


For  organizational  purposes  and  because  of  total  size,  the  POGO  data 
analysis  software  is  divided  into  five  segments.  Each  segment  is  stored 
separately  on  the  5^510  disc  and  runs  on  the  5^510  as  an  individual 
program.  The  5^51C  Fourier  operating  system  allows  nesting  of  as  many 
as  ten  programs.  This  is  the  feature  that  allows  these  five  programs  to 
operate  together  and  perform  the  total  POGO  data  analysis.  The  software 
Includes  two  frequency  response  function  [H(f)]  calculation  routines. 
These  two  data  reduction  methods  are  referred  to  as  the  RATIO  and  DIRECT 
methods  for  calculating  frequency  response  functions.  In  the  RATIO 
technique,  all  H(f)s  are  Initially  referenced  to  the  pulser.  These 
H(f)s  are  then  divided,  with  the  resulting  quotient  being  the  H(f) 
across  various  combinations  of  devices.  For  example,  the  H(f)  between 
devices  2 and  3 In  the  following  schematic  may  be  calculated  as  the 
RATIO  of  the  H (f ) between  devices  1 and  3 and  the  H(f)  between  devices  1 
and  2. 


H(f) 


H(f) 


3/1 


3/2 


HlfT 


2/1 


1 

where  H(f)  - 7:— and  n = output,  1 - input  [1,  2].  Here 
n/ I ^ 

cross  spectrum  between  the  Input  and  output,  and  G^j  is  the 
of  the  input  (pulser). 


G„,  is  the 
n I 

autospectrum 


NOTE;  The  POGO  software  outputs  the  autof.  (power)  spectra  and 
cross  (power)  spectra  as  functions  normalized  to  the  bandwidth  (Af) 
or  as  power  spectral  densities  (PSDs)  and  cross-power  spectral 
densities  (XPSDs).  However,  the  5^510  calculates  and  uses  the 
functions  in  their  nonnormal ized  form  to  calculate  frequency 
response  functions  [H(f)s]  and  coherence  functions  [Y^(f)s]  [2]. 

In  the  remainder  of  this  document,  all  spectral  functions  will  be 
referred  to  as  PSDs  and  XPSDs  regardless  of  their  actual  state  of 
normalization.  The  reader  should  understand  that  all  H(f)  and 
Y^(f)  calculations  are  made  prior  to  normal izat ion. 
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To  calculate  the  y^(f)  that  is  compatible  with  the  RATIO  H(f),  the 
following  formula  is  used  [3]: 


Y^(f) 


3/2 


Y^(f)3/1  Y^(f)2/1 


- 1 


1 

where  Y^(i^)  n/1  “ G,,g'  and  n and  1 are  output  and  input,  as  above  [1 , 


2]. 


1 VI 

1 1 nn 


To  arrive  at  y/2  direct  method  involves 

calculating  the  frequency  response  function  as  the  XPSD,  between  devices 
3 and  2,  and  divided  by  the  PSD  of  device  2.  The  coherence  is  calcu- 
lated as  the  magnitude  squared  of  the  XPSD  divided  by  both  the  input  and 
output  PSDs.  These  equations  are  as  follows: 


H(f) 


rMf) 


3/2 


3/2 


;'32 

*^22 

G G 
22  3 3 


where  ~ XPSD,  0^2  ~ input  PSD,  and  * output  PSD. 


The  POGO  software  also  calculates  frequency  response  functions,  both 
RATIO  and  DIRECT,  which  have  been  blanked,  based  on  either  the  coherence 
or  the  pulser  signal.  Blanking  means  that  the  frequency  response  func- 
tion magnitude  arid  phase  are  set  to  zero  at  those  frequencres  where  the 
value  of  the  coherence,  or  pulser  signal,  does  not  exceed  the  user  set 
minimum. 


As  previously  stated,  the  POGO  analysis  software  consists  of  five  seg- 
ments, each  of  which  is  stored  on  the  5^51C  disc  in  a separate  j<eyboard 
£rogram  f_ile  (KPF)  record.  They  are  currently  stored  in  successive 
records  (1  - 5)  on  the  Marshall  Space  Flight  Center,  Systems  Dynamics 
Lab's  5^5 1C  Fourier  software  disc. 
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KEYBOARD  PROGRAM  FILE,  RECORD  1 

KPFl  contains  part  1 of  three  parts,  which  perform  the  RATIO  calcula- 
tions. Part  1 sets  up  the  data  block  size,  calibration  factors,  total 
amount  of  data  to  be  taken,  and  delta  f (Af)  [frequency  resolution  or 
bandwidth].  This  information  Is  used  to  calculate  both  the  RATIO  and 
DIRECT  H(f)s.  The  5^51C  has  a block  arithmetic  mode,  which  allows  any 
block  of  data  to  be  multiplied  or  divided  by  any  number.  This  number, 
however,  must  be  an  integer.  If  this  number  happens  to  be  a floating- 
point variable  parameter  (VP),  the  5^51G  simply  truncates  any  part  of 
the  VP  that  lies  to  the  right  of  the  decimal.  In  reference  to  this 
problem,  the  POGO  analysis  software  first  multiplies  each  calibration 
factor  by  100,  then  calibrates  the  appropriate  data  block,  and  finally 
divides  the  data  block  by  100.  This  allows  for  two  decimal  places  In 
each  calibration  factor.  The  same  type  procedure  is  used  with  regard  to 
the  VP  that  holds  the  Af  parameter,  to  allow  for  one  digit  to  the  right 
of  the  decimal  when  calculating  PSDs,  After  the  initial  setup,  KPFl 
reads  the  analog  data  from  the  tape  and  stores  it  on  the  ADC  Throughput 
File,  beginning  at  record  0.  KPFl  then  retrieves  this  time  domain  data, 
two  channels  at  a time,  and  calculates  PSDs  for  each  channel.  KPFl  also 
calculates  H(f),  Y^(f) » and  XPSDs  for  the  three  pairs,  LPOP/pulser, 
HPOP/pulser,  and  MCC/pulser.  This  data  is  stored  in  data  blocks  12 
through  27.  Program  control  is  then  transferred  to  KPF2,  which  is  part 
2 of  the  POGO  analysis  software. 

KEYBOARD  PROGRAM  FILE,  RECORD  2 

KPF2  begins  the  calculations  for  the  RATIO  H{f)s  and  y^(f)s.  The  5^51C 
has  preprogrammed  keys  that  will  perform  the  arithmetic  necessary  to 
make  these  calculations.  The  keys,  however,  perform  what  will  be 
referred  to  as  block  arithmetic.  The  5^51C  stores  data  in  blocks  inter- 
nally, in  a format  Hewlett-Packard  calls  floating-point-by-block.  This 
means  that  the  data  block  is  actually  a group  of  integers,  all  with  a 
common  scale  factor.  Because  of  this  method  of  data  handling,  the 
arithmetic  operation  performed  by  the  block  command  is  not  floating 
point,  but  is  integer  arithmetic.  These  operations  involved  in  the 
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RATIO  calculations  are  limited  to  a 40-dB  dynamic  range.  When  inverting 
a number,  as  that  number  approaches  zero,  the  inverse  approaches  infinity. 
It  becomes  apparent  that  if  there  are  any  numbers  whose  value  is  less 
than  0.01  in  a data  block,  any  number  greater  than  or  equal  to  1 will  be 
more  than  kO  dB  down  when  that  block  is  inverted.  The  5^51 C autana tic- 
ally  scales  up  the  block  scale  factor  to  eliminate  overflowing,  which 
causes  any  number  more  than  kO  dB  down  to  be  set  to  zero.  This  is  a 
serious  problem  when  doing  the  type  of  calculation  involved  in  arriving 
at  the  RATIO  H(f)s  and  y^(f)s.  in  both  of  these  calculations,  espe- 
cially the  Y^(f)  calculation,  the  quotient  in  the  formulations  is  domi- 
nated by  the  small  numbers  in  the  denominator.  These  small  numbers  are 
not  reliable  data  points,  and  the  valid  data  Is  lost  below  the  40-dB 
dynamic  range. 

The  problem  was  solved  in  this  program  by  using  a floating-point  arith- 
metic scheme  In  the  y^Cf)  calculations  and  a combination  of  floating- 
point arithmetic  and  a filtering  technique  in  the  H(f)  calculations. 

The  floating-point  scheme  Involves  extracting  each  channel  (or  data 
point)  polnt-by-polnt  from  a data  block,  putting  the  value  into  a 
floating-point  format  (via  a floating-point  variable  parameter),  exe- 
cuting the  arithrretic  operation,  and  finally  putting  the  datum  back  Into 
the  data  block.  A magnitude  filtering  technique  Is  used  In  the  H(f) 
calculation  to  eliminate  a loss  of  data  after  the  values  are  put  back 
into  the  data  block.  The  real  and  Imaginary  parts  of  the  H(f)  are 
filtered  so  that  any  values  outside  the  range  of  +7  are  set  to  zero. 

This  allows  a polar  magnitude  maximum  of  10.  The  40-dB  dynamic  range 
then  will  be  10  ^ to  10^.  This  filtering  technique  Is  not  necessary  in 
the  Y^(f)  calculations  because  of  the  final  quantity  inversion  In  the 
formulation.  All  the  values  to  be  inverted,  however,  are  scanned  for 
values  of  identically  zero  before  inversion.  This  is  to  eliminate  any 
numbers  actually  going  to  the  upper  limit  (Infinity). 

The  program  loops  through  the  above  operation  enough  times  to  calculate 
the  RATIO  H(f)  and  y^{f)  for  each  channel  in  the  data  block.  This  loop 
is  nested  inside  a loop,  which  allows  for  all  three  sets  of  ratio 
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calculations  to  be  computed.  Part  2 of  the  P060  analysis  software  then 
transfers  control  to  part  3. 

KEYBOARD  PROGRAM  FILE,  RECORD  3 

KPF3  defines  the  variable  parameter,  corresponding  to  the  ratio  calcula- 
tions, necessary  to  operate  the  plot  routine,  KPF3  also  performs  the 
calculations  involved  in  the  H(f)  blanking  routine.  KPF3  has  several 
entry  points,  which  allot  various  operations  to  be  performed  outside  the 
normal  program  flow.  KPF2  flows  into  the  beginning  of  KPF3,  which  jumps 
to  the  blanking  routine.  The  blanking  routine  begins  by  prompting  the 
user  to  establish  whether  the  blanking  will  be  based  upon  the  y^(f)  or 
the  pulser  PSD.  The  user  will  enter  either  a zero,  to  establish  coher- 
ence blanking,  or  a two,  to  define  pulser  blanking.  This  number  (O  or 
2)  will  also  be  used  to  establ i sh  which  labels  will  be  given  to  the 
plots  of  the  blanked  H(f)s  during  the  plot  routine.  If  pulser  blanking 
is  chosen,  KPF3  will  display  the  pulser  PSD  to  allow  the  user  an  oppor- 
tunity to  view  the  signal  and  decide  on  the  minimum  value  used  for  the 
blanking  routine.  If  Y^(^)  blanking  is  chosen,  KPF3  skips  directly  to 
the  next  step,  which  is  a prompt,  asking  for  the  minimum  value  in  the 
blanking  routine.  Next,  KPF3  sets  up  the  block  that  will  be  used  for 
generating  the  blanked  H(f).  KPF3  then  enters  the  blanking  loop.  This 
loop  operates  individually  on  each  channel  in  the  frequency  domain  data 
block.  In  this  loop,  KPF3  goes  to  either  the  pulser  or  Y^(f)  block  and 
gets  one  channel  at  a time,  beginning  with  channel  one  (channel  zero  is 
the  first  channel)  and  ending  with  the  next  to  last  channel  in  the  data 
block.  The  value  of  this  channel  is  compared  with  the  values  of  the 
channel  immediately  preceding  and  the  channel  immediately  following.  If 
the  value  of  the  channel  in  question  is  greatest,  it  is  determined  to  be 
a local  maximum.  The  value  is  then  compared  with  the  minimum  value 
specified  earlier  by  the  user.  If  the  value  of  this  channel  is  greater 
in  all  three  tests,  it  is  chosen  as  a channel  of  interest.  The  corre- 
sponding channel  in  the  H(f)  data  block  is  gathered  into  a complex  vari- 
able parameter  and  stored  in  the  data  block  that  was  prepared  previously 
for  the  blanked  H(f).  When  the  loop  is  complete,  the  blanked  H(f)  will 
consist  of  calculated  H(f)  data  only  at  those  frequencies  where  the 


11 


pulser  or  y^(^)  passed  all  the  necessary  tests.  The  remainder  of  the 
blanked  H(f)  block  will  be  filled  with  zeros.  Because  the  HP  5^51 C 
stores  frequency  domain  data  with  both  real  and  imaginary  information  in 
a combined  complex  channel,  the  blanked  H(f)  contains  blanked  phase 
information  as  well  as  blanked  magnitude.  KPF3  then  loops  three  times 
to  generate  a blanked  version  of  each  H(f),  HPOP/LPOP,  MCC/HPOP,  and 
MCC/LPOP.  KPF3  then  jumps  to  the  section  that  defines  the  variable 
parameters  necessary  to  run  the  plot  program.  Finally,  KPF3  transfers 
control  to  KPF5,  the  POGO  plot  routine. 

KEYBOARD  PROGRAM  FILE,  RECORD  4 

KPFk  is  a single-part  program,  which  calculates  the  H(f)s,  HPOP/LPOP, 
MCC/HPOP,  and  MCC/LPOP  by  matching  the  outputs  and  inputs  "directly” 
rather  than  ratioing  H(f)s,  which  have  all  been  referenced  to  the 
pulser.  This  program  (KPF4)  uses  several  variable  parameters,  which 
must  have  been  previously  defined  by  running  the  first  three  files 
containing  the  RATIO  software.  These  variable  parameters  are  VPl  = 
block  size;  VP2  = number  of  ADC  Throughput  records;  VP2000  through 
VP2003  * calibration  factors  for  channels  A through  D,  respectively; 
VP2004  = frequency  resolution  (Af);  and  VP2030  = minimum  value  for  the 
blanking  routine.  KPF4  also  retrieves  time  domain  data  from  the  ADC 
Throughput  File.  These  data  must  have  been  written  previously  on  the 
disc  (starting  with  record  zero)  via  KPFl  (RATIO  software,  part  1).  The 
data  may  be  written  to  the  disc  and  the  necessary  variable  parameters 
defined  independently  of  the  RATIO  software.  It  is  recommended,  how- 
ever, that  all  programs  be  operated  in  sequence.  The  direct  POGO 
analysis  software  (KPF4)  writes  over  the  data  blocks  generated  via  the 
RATIO  software.  This  program  therefore  should  not  be  initiated  until 
all  the  required  ratio  data  blocks  have  been  reproduced.  These  data 
blocks  include  all  PSDs  and  XPSDs  (ref.  pulser),  which  are  not  recal- 
culated by  the  DIRECT  software.  Immediately  after  clearing  the  data 
blocks,  KPF4  begins  a loop  that  operates  three  times.  The  loop  reads 
the  correct  data  channels  from  the  disc,  Fourier  transforms  the  data, 
calibrates  the  data,  and  calculates  averaged  PSDs  via  the  special  sub- 
routine of  the  5451c.  KPF4  then  calculates  the  H(f)  and  y^(f)  from  the 
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averaged  PSD  data.  The  H(f),  y^(f),  and  XPSDs  are  then  stored  In  the 
appropriate  data  blocks.  (See  "DIRECT  POGO  Program  Menu,"  appendix  C.) 

KPF4  now  transfers  control  to  KPF3  at  the  point  where  the  H(f)  blanking 
loop  begins.  Before  jumping  to  KPF3,  KPF4  defines  the  variable  param- 
eters that  describe  the  location  of  the  DIRECT  H(F)s.  Once  In  KPF3  the 
blanked  H(f)s  are  calculated  using  the  same  minimum  value  specified  by 
the  user  during  the  RATIO  calculations.  If  the  user  has  entered  the 
POGO  software  in  the  DIRECT  calculations  and  wishes  to  change  the  mini- 
mum blanking  value,  the  user  simply  jumps  to  label  43  of  KPF4. 

Upon  returning  from  the  blanking  routine,  KPF4  defines  the  variable 
parameters  necessary  to  operate  the  POGO  plot  program  (KPF5).  KPF4  then 
transfers  control  to  the  plot  routine. 

KEYBOARD  PROGRAM  FILE,  RECORD  5 

KPF5  is  entitled  "POGO  Plot  Software."  This  software  operates  with  both 
the  RATIO  and  the  DIRECT  software.  KPF5  uses  several  variable  param- 
eters that  must  be  defined  by  the  respective  software  whose  data  is  to 
be  plotted.  The  program  is  organized  into  several  separate  routines, 
including  a plot  routine,  a hard  copy  routine,  a delay  routine,  and  a 
routine  that  sets  up  and  plots  all  the  H(f)s  in  three  different  formats. 
The  program  sets  up  and  automatically  plots  all  the  data  generated  by 
the  respective  RATIO  and  DIRECT  programs.  The  KPF5  prints  a data  block 
menu  and  allows  the  user  to  replot  any  of  the  data  blocks.  The  plot 
program  is  geared  to  run  with  a Tektronics  4052  graphics  terminal  and  a 
Tektronics  hard  copy  unit.  Because  of  various  Incompatibilities  between 
the  5451c  and  the  Tektronics  terminal  with  respect  to  the  time  required 
to  transmit  plot  conmands  and  the  time  required  to  execute  them,  a delay 
routine  is  incorporated  In  the  plot  program.  This  subroutine  simply 
counts  from  one  to  50  and  returns.  The  routine  is  called  In  several 
locations  where  a time  delay  between  plot  commands  is  necessary. 

The  plot  program  uses  several  variable  parameters  that  establish  which 
program  (RATIO  or  DIRECT)  has  generated  the  data  being  plotted.  The 
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plot  program  also  labels  each  plot  with  test  ID  information  in  the  upper 
right-hand  corner  of  the  plot.  This  information  is  stored  as  message 
number  one  in  each  of  the  text  buffers  one  and  two.  Text  buffer  #1  is 
used  for  the  RATIO  data,  while  text  buffer  #2  is  used  with  the  DIRECT 
data.  Message  #1  should  be  changed  for  each  new  set  of  test  data.  The 
message  may  be  changed  by  executing  the  following  set  of  operations  from 
the  5^51C  terminal; 


Step  # Command  5^51 C Status 

1 Y 5403  (nl)  busy 

This  command  calls  text  buffer  nl  (1  for  RATIO,  2 for  DIRECT) 
into  the  computer  core  and  activates  the  editing  commands. 

2 /ROl  busy 

This  command  tells  the  5451 C to  replace  message  #1. 

3 01  busy 

4 (Test  ID)  busy 

5 /*  busy 

This  sequence  of  commands  results  in  the  test  ID  being  recognized 
as  message  01.  The  /*  defines  the  end  of  the  message. 

6 / busy 

This  command  terminates  the  text  buffer  edit  mode. 


•To  allow  the  user  to  enter  the  program  flow  for  special  purposes,  the 
P060  software  has  many  labels  throughout  the  five  KPFs.  Listed  below 
are  the  various  jump  commands  and' the! r uses. 

J 0 1 

This  command  starts  the  POGO  software,  initially  setting  up  the 
calibration  factors,  etc,  and  reading  the  analog  data. 

J 10  1 

This  command  enters  POGO  Software  Part  1 after  all  the  initial 
values  have  been  entered  and  the  analog  data  has  been  stored  on 
disc  in  the  ADC  Throughput  File.  If  this  command  is  to  be  executed, 
the  user  must  first  be  sure  that  the  variable  parameters  that  con- 
tain the  values  for  the  blockslze,  calibration  factors,  number  of 
records  of  data,  and  frequency  resolution  (Af)  are  defined  cor- 
rectly and  that  the  correct  time  domain  data  is  stored  on  the  ADC 
Throughput  File. 
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J 20  2 


This  command  begins  the  RATIO  calculations  for  H(f)s  and  Y^(f)s. 
This  command  may  be  executed  if  the  proper  data  blocks  reside  in 
core  In  the  correct  block  numbers.  The  data  blocks  necessary  are 


H(f) 

LPOP/Pulser 

Block  12 

Y^(f) 

II 

Block  1 3 

H(f) 

HPOP/Pulser 

Block  14 

y2(f) 

II 

Block  15 

H(f) 

MCC/Pul ser 

Block  16 

Y^(f) 

il 

Block  17 

H 30  3 

This  command  may  be  used  when  the  user  wishes  to  recalculate  the 
blanked  RATIO  frequency  response  functions.  The  user  will  have  the 
opportunity  to  re-establish  the  type  of  blanking  [pulser  or  y^(f)] 
and/or  the  minimum  value  for  that  blanking.  The  RATIO  autoplot 
parameters  are  defined  and  the  program  control  is  transferred  to 
the  plot  routine. 

J 31  3 

This  command  is  preliminary  to  the  POGO  Plot  Routine.  This  command 
enters’ KPF3  and  defines  the  variable  parameters  necessary  to  run 
the  automatic  plot  portion  of  the  plot  routine.  After  these  param- 
eters are  defined,  program  control  Is  automatically  transferred 
to  the  plot  routine.  This  command  should  be  used  only  when  the 
user  wishes  to  enter  the  autoplot  routine  with  RATIO  data  blocks. 
This  command  might  be  useful  if,  for  example,  the  user  wants  to 
change  the  test  JD  message  in  the  RATIO  text  buffer  after  the 
POGO  software  has  completed  the  data  calculations.  The  user  can 
stop  program  execution  and  edit  the  text  buffer,  then  return  to 
the  autoplot  sequence. 

J kO  k 

This  command  begins  the  DIRECT  POGO  Software.  This  should  be  used 
only  after  the  RATIO  program  has  been  run,  defining  the  appro- 
priate variable  parameters,  and  writing  the  time  domain  data  on 
the  ADC  Throughput  File. 

J 43  4 

This  command  enters  KPF4  where  necessary  variable  parameters  are 
defined  just  prior  to  transferring  control  to  the  H(f)  blanking 
routine  in  KPF3.  This  command  will  allow  the  user  to  calculate 
DIRECT  blanked  H(f)s. 
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J 44  4 


This  command  causes  the  computer  to  begin  execution  In  KPF4  at 
the  point  where  the  variable  parameters  that  control  the  auto- 
plot routine  for  DIRECT  data  blocks  are  defined.  Program  control 
is  automatically  transferred  to  the  plot  routine. 

J 53  5 

This  command  will  cause  the  plot  program  to  display  on  the  termi- 
nal the  appropriate  menu  corresponding  to  the  data  stored  in  the 
data  blocks.  The  plot  program  prompts  the  user  to  set  the  display 
scale  as  desired,  then  press  the  continue  button  to  allow  the  pro- 
gram to  plot  the  data  block  on  the  graphics  terminal.  This  com- 
mand is  useful  when  the  user  wishes  to  return  to  the  plot  program 
after  the  machine  has  been  Idled  and  the  plot  routine  pointers 
have  been  changed. 


Appendix  A describes  the  flow  of  the  POGO  software.  Appendix  B contains 
a complete  list  of  all  the  program  steps  along  with  a description  of 
each.  Appendix  C contains  a listing  of  the  text  record  messages  used 
with  the  POGO  software.  Appendix  D contains  a listing  of  both  text 
buffers  used  with  the  POGO  software.  Appendix  E contains  sample  data 
output  from  the  POGO  software. 
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APPENDIX  A 


A1 


POGO  SOFTWARE  FLOWCHART 


RATIO  POGO  SOFTWARE  - PART  1 FLOWCHART 
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RATIO  POGO  SOFTWARE  - PART  1 FLOWCHART  (Continued) 


Steps  120  - 125 

Display  the  prompt  - "Start  the 
data  tape  and  press  CONTINUE" 


^^ep  129 
'vv^Pause 


Steps  132  - 136 

Initialize  the  ADC  Throughput  and 
acquire  VP2  records  of  four  channels 
each 


Step  142 
Label  10 


Steps  146  - 176 

Set  up  the  variable  parameters 
for  the  data  handling  loops 


Step  182 
Label  11 


Steps  186  - 213 
Set  up  the  H(f)  loop 


Step  219 

Calculate  the  number  of  the  VP  which 
contains  the  appropriate  cal  factor 
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RATIO  POGO  SOFTWARE  - PART  1 FLOWCHART  (Continued) 

(p 

Step  226 

Initialize  the  off-line  BSFA 


Steps  237  - 266 

Get  the  analog  data  from  the  disc, 
convert  it  to  frequency  domain, 
calibrate  it,  and  leave  the  input 
in  Block  1 and  output  in  Block  2. 
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A4 


RATIO  POGO  SOFTWARE  - PART  1 FLOWCHART  (Concluded) 


Step  377 

Loop  for  all  3 sets  of  POGO  data 

1 

r 

Step  383 

Jump  to  Label  20  Keyboard  Program  2 

T 
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A5 


RATIO  POGO  SOFTWARE  - PART  2 FLOWCHART 


Step  146 


Initialize  the  channel  counter 


0 


Steps  156  - 177 

Get  the  H(f)  and  Y^(f)  data 
for  the  ratio  calcuTation, 
one  channel  at  a time. 

I 
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RATIO  POGO  SOFTWARE  - PART  2 FLOWCHART  (Continued) 


A7 
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A8 


RATIO  POGO  SOFTWARE  - PART  2 FLOWCHART  (Concluded) 
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POGO  SOFTWARE  BLANKING  ROUTINE 

X^tep  1 Ny 
( Label  30  ) 


Steps  5-23 

Set  up  variable  parameters  for  ratio 
H(f)  blanking;  jump  to  label  32 


Label  31 


Steps  31  - 91 

Set  up  variable  parameters  to  plot 
ratio  data  with  POGO  plot  routine 


Step  97 

Jump  to  label  50  in  KPF5  (plot  routine) 


.x^tep  102 
Label  32 
Start  Blanking 
V Routine 


Steps  106  - 115 

Display  the  prompt  and  allow  the  user 
to  determine  the  type  of  blanking 
routine  (coherence  or  pulser  based) 


POGO  SOFTWARE  BLANKING  ROUTINE  (Continued) 


Step  128  - 

Define  the  VP31  for  pulser  blanking 


Steps  135  - 143 

Display  the  correct  data  block  and  ask 
the  user  to  set  the  minimum  blanking 
value 


Steps  148  - 155 

Define  VP39  for  blanking  routine 


Step  164 


Convert  the  block  to  be  blanked  to 
rectangular  coordinates 


Steps  168  - 179 

Clear  a block  and  establish  qualifiers 
for  the  blanked  H(f) 


Default  VP31  to 
coherence  blanking 


n 
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All 


P060  SOFTWARE  BLANKING  ROUTINE  (Continued) 


27 


No 


0 


DIRECT  POGO  SOFTWARE  FLOWCHART 


AI3 


29 


DIRECT  POGO  SOFTWARE  FLOWCHART  (Continued) 


POGO  SOFTWARE  PLOT  ROUTINE 

^ Step  1 

( Start  - Label  50  | 


32 


33 


POGO  SOFTWARE  PLOT  ROUTINE  (Continued) 


( 130  > 
1 Return , 


Steps  120  - 167 

Display  plot  menu,  determine  label 
numbers  for  blanked  plots,  and  allow 
use-to  select  a block  number  to  plot 


^ Step  171  ^ 

Pause  to  set  the 
scale  factor V 


Steps  175  - 179 

Set  up  the  plot  and  jump 
to  the  plot  routine 


^00^ 


II 


Step  183 

Jump  to  the  hard  copy  routine 


^120^ 

^turi^ 


Steps  187  - 191 

Pause  - when  user  continues 
jump  to  Label  53 


POGO  SOFTWARE  PLOT  ROUTINE  (Continued) 


/Step  1 95  \ 
VLabel  100/ 


Steps  199  - 227 

Plot  the  selected  data  block  along 
with  data  and  test  annotation 


Step  232 
^Subreturn, 


Step  235  ^ 
Label  IIQ 


Steps  239  - 255 

Determine  label  numbers  for 
blanked  plots 


Step  262 

Convert  the  H(f)  data  to 
polar  coordinates 


Steps  266  - 273 

Set  up  to  plot  H(f)  phase  and 
jump  to  the  plot  routine 
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APPENDIX  B 

POGO  SOFTWARE  LISTING 
KEYBOARD  FILE  1 - RATIO  POGO  SOFTWARE  PART  1 


1 L 

0 

5 BS 

4096 

9 CL 

12  CL 

1 

16  CL 

2 

20  CL 

3 

24  CL 

4 

28  CL 

5 

32  CL 

6 

36  Y 

5814 

40  MS 

34 

44  MS 

14 

48  Y R 

1 

53  BS 

ID 

57  MS 

14 

61  Y R 

2000 

2003 

67  Y * 

2000 

2000D 

100 

74  Y * 

2001 

2001D 

TOO 

81  Y * 

2002 

2002D 

100 

88  Y * 

2003 

2003D 

TOO 

95  MS 

14 

99  Y R 

2 

104  MS 

14 

108  Y R 

2004 

tThls  means  the  variable  of  parameter  1 


Label  0 

Set  blocksize  to  4096 


Clear  blocks  0 through  6 


Call  user  program  5414  (blank 
terminal  screen) 

Set  text  record  pointer  to  zero 
(default) 

Read  next  text  record  (message  #1) 

Input  from  terminal  VPlt  . 
(blocksize) 

Set  blocksize  to.  VPl 

Read  next  text  record  (message  #2) 

Input  from  terminal  VP2000, 

VP2001,  VP2002,  and  VP2003 
(call brati on  factors ) 

VP2000  - VP2000  * TOO 

VP2001  = VP2001  * TOO 

VP2002  = VP2002  * TOO 

VP2003  = VP2003  * 100 

Read  next  text  record  (message  #3) 

Input  from  terminal  VP2  (number 
of  data  acquisition  loops) 

Read  next  text  record  (message  #4) 

Input  from  terminal  VP2004  (Af) 
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B2 


KEYBOARD  FILE  1 - RATIO  P060  SOFTWARE  PART  1 (Continued) 


113  Y * 

2004 

2004D 

10 

VP2004  = VP2004  * TO 

120  MS 

34 

20 

Set  text  record  pointer  to  20 

125  MS 

14 

Read  next  text  record  (message  47) 

129  D 

Display  block  0 (used  as  program 
pause) 

132  MS 

32 

Set  ADC  Throughput  File  pointer 
to  0 

136  MS 

22 

4 

2D 

Write  analog  data  to  the  ADC 
throughput  file.  No.  of 
records  = 2D  (4  channels  each) 

142  L 

10 

Label  10 

146  Y _ 

3 

18 

Set  VP3  to  18 

152  Y _ 

4 

19 

Set  VP4  to  19 

158  Y _ 

5 

22 

Set  VP4  to  22 

164  Y _ 

6 

23 

Set  VP6  to  23 

170  Y _ 

7 

2 

Set  VP7  to  2 

176  Y _ 

8 

12 

Set  VP8  to  12 

182  L 

11 

Label  11 

186  CL 

189  CL 

1 

193  CL 

2 

197  CL 

3 

Clear  blocks  0 through  6 

201  CL 

4 

205  CL 

5 

209  CL 

6 

213  Y _ 

9 

3 

Set  VP9  to  3 

219  Y A+ 

10 

7D 

1999 

VPIO  = VP7  + 1999 

226  Y 

41 

0 

0 0 

Call  user  program  41;  nl=0= 
center  frequency  of  the  band  of 
interest,  n2=0=width  of  the  band 
of  interest,  n3=0=ADC  throughput 
s ta rt  record.  (Ini ti a 1 i zes  off- 
line BSFA) 

233  L 

12 

Label  12 

237  Y 

45 

2 

0 7D 

Call  user  program  45;  nl=2=block 

nurrtDer  where  the  data  will  be 
stored,  n2=0=number  of  Hannings 
applied  to  the  data,  n3=VP7= 
channel  that  will  be  read  from 
the  disc.  (Off-line  BSFA) 
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B3 

KEYBOARD  FILE  1 - RATIO  POGO  SOFTWARE  PART  1 (Continued) 


244 

Y 

45 

1 

0 1 

Call  user  program  45;  nl=l, 
n2=0,  n3=l 

251 

* 

2 

-lOD 

Multiply  block  2 by  VP(VPIO)^ 

256 

; 

2 

100 

Divide  block  2 by  100 

261 

* 

i 

2000D 

Multiply  block  1 by  VP2000 

266 

: 

1 

100 

Divide  block  1 by  100 

271 

SP 

1 

2 

2 

Compute  the  average  auto-  and 
cross-power  spectrums,  where 
block  nl=input  data  block  and 
block  nl+l=output,  n2=2  speci- 
fies dual  channel,  n3=2  speci- 
fies double  precision 

277 

# 

12 

2D 

0 

Loop  through  label  12  VP2  times 

283 

CH 

1 

2 

2 

Calculate  transfer  function  data 
where  nl -input  data  block,  n2=2 
specifies  dual  channel,  n3=2 
specifies  double  precision 

289 

X< 

1 

Load  block  1 into  block  0 

293 

X> 

8D 

store  block  0 into  block  VPS 
[H(f)] 

297 

Y A+ 

8 

VPS  = VP8  + 1 

303 

X< 

2 

Load  block  2 into  block  0 

307 

X> 

8D 

Store  block  0 into  VP8 

311 

L 

13 

Label  13 

315 

• 

• 

9D 

2004D 

Divide  block  VP9  by  VP2004 

320 

* 

9D 

20 

Multiply  block  VP9  by  20 

325 

X< 

9D 

Load  block  VP9  into  block  0 

329 

x> 

-9D 

Store  block  0 into  block  VP(VP9) 

333 

Y A+ 

9 

VP9  = VP9  + 1 

339 

# 

13 

4 

0 

Loop  through  label  13  four  times 

345 

Y A+ 

4 

VP4  = VP4  + 1 

351 

Y A+ 

5 

5D 

2 

VPS  = VPS  + 2 

358 

Y A+ 

6 

6D 

2 

VP6  = VP6  + 2 

365 

Y A+ 

7 

VP7  = VP7  + 1 

The  minus  sign  in  front  of  the  VP  number  means  to  use  the  value  in  the  VP 
that  is  specified  by  the  value  of  the  VP  given  (i.e.,  the  value  of  the  VP 
specified  by  VPIO). 
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KEYBOARD  FILE  1 - RATIO  PQGO  SOFTWARE  PART  1 (Concluded) 


371  Y A+  8 

377  # 11  3 0 

383  J 20  2 

388  . 


VP8  = VP8  + 1 

Loop  through  label  11  three  times 

Jump  to  label  20  in  keyboard 
stack  2 

End 


kk 
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KEYBOARD  FILE  2 - RATIO 


1 L 

20 

5 TR 

12 

9 TR 

14 

13  TR 

16 

17  Y _ 

3 

3 

23  Y _ 

4 

6 

29  Y _ 

5 

7 

35  Y _ 

6 

14 

41  Y _ 

7 

12 

47  Y _ 

8 

16 

53  Y _ 

9 

14 

59  Y _ 

10 

16 

65  Y _ 

11 

12 

71  Y _ 

12 

13 

77  Y _ 

13 

15 

83  Y : . 

14 

ID 

2 

90  Y A+ 

14 

95  Y _ 

15 

17 

101  Y _ 

17 

13 

107  L 

21 

111  Y TR 

21 

-5D 

0 

118  Y TR 

21 

3D 

1 

125  Y TR 

21 

4D 

1 

132  Y TR 

21 

-12D 

0 

139  Y TR 

21 

5D 

1 

146  Y _ 

0 

0 

152  L 

22 

156  Y X< 

3000 

-4D 

OD 

POGO  SOFTWARE  PART  2 
Label  20 

Convert  blocks  12  \ 

Convert  blocks  14  > ^ rectangular 

convert  diocks  I4  j coordinates 

Convert  blocks  16  / 

Set  VPS  to  3 
Set  VP4  to  6 
Set  VP5  to  7 
Set  VP6  to  14 
Set  VP7  to  12 
Set  VPS  to  16 
Set  VP9  to  14 
Set  VPIO  to  16 
Set  VPll  to  12 
Set  VP12  to  13 
Set  VP13  to  15 
VP14  = VPl/2 
VP14  = VP14  + 1 
Set  VP15  to  17 
Set  VP17  to  13 
Label  21 

Get  the  block  qualifiers  from 
block  VP(VP5)  and  put  them  in 
VP21  - VP25 

Put  the  block  qualifiers  stored 
in  VP21  - VP25  into  block  VPS 

Put  the  block  qualifiers  into 
block  VP4 

Get  the  block  qualifiers  from 
block  VP(VP12)  and  put  them  into 
VP21  - VP25 

Put  the  block  qualifiers  into 
block  VP5 

Set  VPp  to  0 

Label  22 

Get  the  complex  data  from  block 
VP(VP4),  channel  VPO,  and  put 
the  values  into  complex  VP3000 


As 


B6 


KEYBOARD 

FILE 

2 - RATIO 

POGO 

SOFTWARE  PART  2 (Continued) 

16  3 Y X< 

3001 

-5D 

OD 

Get  the  complex  data  from  block 

VP(VP5),  channel  VPO  and  put 

them  into  complex  VP3001 

170  Y X< 

2005 

12D 

OD 

Get  the  real  data  from  block 

VPT2,  channel  VPO  and  put  it 

into  VP2005 

177  Y X< 

2006 

-12D 

OD 

Get  the  real  data  from  block 

VP(VP12),  channel  VPO  and  put 

it  into  VP2006 

184  Y IF 

2005 

0 

7 

0 

If  VP2005  = 0,  skip  seven  steps 

192  Y : 

2005 

1 

2005D 

VP2005  = 1/VP2005 

199  Y If 

2006 

0 

5 

0 

If  VP2006  = 0,  skip  five  steps 

207  Y : 

2006 

1 

2006D 

VP2006  - 1/VP2006 

214  Y A+ 

2007 

2005D 

2006D 

VP2007  = VP2005  + VP2006 

221  Y A- 

2007 

VP2007  = VP2007  - 1 

226  Y : 

2007 

1 

2007D 

VP2007  = 1/VP2007 

233  Y IF 

2007 

0 

1 

2 

If  VP2007  > 0,  skip  one  step 

241  Y _ 

2007 

0 

Set  VP2007  to  0 

247  Y X> 

2007 

5D 

00 

. Put  VP2007  into  the  real 

part 

of  block  VP5j  channel  VPO 

254  Y * 

2009 

2900D 

2902D 

VP2009  = VP2900  * VP2902  (ac) 

261  Y * 

2010 

2901 D 

2903D 

VP2010  = VP2901  * VP2903  (bd) 

268  Y * 

2011 

2901 D 

2902D 

VP2011  = VP2901  * VP2902  (be) 

275  Y * 

2012 

2900D 

29030 

VP2012  = VP2900  * VP2903  (ad) 

282  Y * 

2013 

2902D 

2902D 

VP2013  - (VP2902)2  (c^) 

289  Y * 

2014 

2903D 

2903D 

VP2014  - (VP2903)2  (d2) 

296  Y A+ 

2015 

201 30 

201 4D 

VP2015  = VP2013  + VP2014  (c2+d2) 

303  Y A+ 

2016 

2009D 

201  OD 

VP2016  = VP2009  + VP2010  (ac+bd) 

310  Y A- 

2017 

201  ID 

20120 

VP2017  = VP2009  - VP201; 

2 (bc-ad) 

317  Y : 

2904 

201 6D 

201 5D 

VP2904  = VP2016/VP2015 

(ac+bd) 

_(c2+d2)_ 

324  Y : 

2905 

201 7D 

201 5D 

VP2905  = VP2017/VP2015 

i 

(bc-ad) 

.(c2+d2)J 

331  Y IF 

2904 

7 

1 

2 

If  VP2094  > 7,  skip  one  step 

339  Y IF 

2904 

-7 

1 

1 

If  VP2904  > -7,  skip  one  step 

347  Y 

2904 

0 

Set  VP2904  to  0 

46 


KEYBOARD  FILE  2 - RATIO  POGO  SOFTWARE  PART  2 (Concluded) 


B7 


S5S 

Y 

IF 

2905 

7 

1 

S61 

Y 

IF 

2905 

-7 

1 

S69 

Y 

2905 

0 

S75 

Y 

X> 

3002 

4D 

OD 

S82 

Y 

A+ 

0 

S87 

# 

22 

140 

0 

S9S 

Y 

A+ 

3 

398 

Y 

A+ 

4 

4D 

2 

405 

Y 

A+ 

5 

5D 

2 

412 

Y 

A+ 

12 

12D 

2 

419 

# 

21 

3 

0 

425 

J 

30 

3 

430 

« 

If  VP2905  > 7,  skip  one  step 

If  VP2905  > -7,  skip  one  step 

Set  VP2905  to  0 

Put  complex  VP3002  into  block 
VP4,  channel  VPO 

VPO  = VPO  + 1 

Loop  through  label  22  VP14  times 

VPS  « VPS  + 1 

VP4  - VP4  + 2 

VP5  - VP5  + 2 

VP12  = VP12  + 2 

Loop  through  label  21  three  times 
Jump  to  label  SO,  keyboard  stack  S 
End 
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KEYBOARD  FILE  3 - RATIO  POGO  SOFTWARE  PART  3 


1 L 

30 

Label  30 

5 Y _ 

31 

7 

Set  VP31  to  7 

n Y _ 

35 

6 

Set  VP35  to  6 

17  Y _ 

36 

3 

Set  VP36  to  3 

23  J 

32 

Jump  to  label  32 

27  L 

31 

Label  31 

31  Y _ 

3 

1 

Set  VP3  to  1 

37  Y _ 

4 

27 

Set  VP4  to  27 

43  Y _ 

5 

16 

Set  VP5  to  16 

49  Y _ 

6 

14 

Set  VP6  to  14 

55  Y _ 

7 

12 

Set  VP7  to  12 

61  Y _ 

8 

10 

Set  VP8  to  10 

67  Y _ 

9 

8 

Set  VP9  to  8 

73  Y _ 

10 

6 

Set  VPIO  to  6 

79  Y _ 

11 

25 

Set  VPll  to  25 

85  Y _ 

12 

24 

Set  VP12  to  24 

91  Y _ 

13 

125 

Set  VP 13  to  125 

97  0 

50 

5 

Jump  to  label  50  in  Keyboard 
File  5 

102  L 

32 

Label  32 

106  MS 

34 

82 

Set  the  text  record  pointer 

111  MS 

14 

Read  next  text  record  (type 
blanking  routine 

115  Y R 

30 

Input  from  terminal  VP30 

120  Y IF 

30 

0 

1 

2 

IF  VP30>0,  skip  1 Step 

128  Y A+ 

31 

31D 

11 

VP31  = VP31  + 11 

135  X< 

31 D 

Load  blanking  block  (VP31) 

139  MS 

14 

Read  the  next  text  record 
(minimum  blanking  value) 

143  Y R 

2030 

Input  from  terminal  VP2030 

148  Y : 

39 

lOD 

2 

VP39  = VPl/2 

155  Y A- 

39 

VP39  = VP39  - 1 

160  L 

33 

Label  33 

164  TR 

35D 

48 


KEYBOARD 

FILE  3 

- RATIO  POGO 

168  CL 

36D 

172  Y TR 

21 

350 

0 

179  Y TR 

21 

36D 

1 

186  Y _ 

32 

0 

192  Y _ 

33 

1 

198  Y _ 

34 

2 

ZQk  L 

34 

208  Y _ 

3000 

0 

214  Y X< 

2032 

3ID 

320 

221  Y X< 

2033 

3ID 

330 

228  Y X< 

2034 

31 D 

340 

235  Y IF 

2033 

203OD 

3 -2 

243  Y IF 

2033 

2032D 

2 -2 

251  Y IF 

2033 

2034D 

1 -2 

259  Y X< 

3000 

35D 

330 

266  Y X> 

3000 

360 

330 

273  Y A+ 

32 

278  Y A+ 

33 

283  Y A+ 

34 

288  # 

34 

390 

294  Y A+ 

35 

350 

2 

301  Y A+ 

36 

306  Y A+ 

31 

310 

300 

313  # 

33 

3 

319  < 

322. 

B9 

SOFTWARE  PART  3 (Concluded) 

Clear  block  BP 36 

Get  the  block  qualifiers  from 
block  VP35  and  put  them  Into 
VP21  - VP25 

Put  the  block  qualifiers  stored 
In  VP21-VP25  into  block  VP36 

Set  VP32  to  0 

Set  VP 33  to  1 

Set  VP34  to  2 

Label  3k 

Set  VP3000  to  0 

Get  the  real  data  from  block 
VP3I , channel  VP32  and  store  It 
In  VP2032 

Get  the  real  data  from  block 
BP31 , channel  VP33  and  store  It 
In  VP2033 

Get  the  real  data  from  block 
VP3I,  channel  VP3^  and  store  It 
in  VP2034 

If  VP2033<VP2030,  skip  3 steps 

If  VP2033<VP2032,  skip  2 steps 

If  VP2033<VP2034,  skip  1 step 

Get  the  complex  data  from  block 
VP35>  channel  VP 3 3 and  store  It 
in  VP3000 

Put  the  complex  data  stored  1n 
VP3000  Into  block  VP36,  channel 
VP33 

VP32  = VP32  + 1 
VP33  = VP33  + 1 
VP34  = VP34  + 1 

Loop  through  label  3^  VP39  times 

VP35  = VP35  + 2 

VP36  » VP 36  + 1 

VP31  = VP 31  + VP30 

Loop  through  label  33  three  times 

Subreturn 

End 
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BIO 

KEYBOARD  FILE  4 - 

1 L 

40 

5 X< 

18 

9 X> 

20 

13  Y _ 

3 

8 

19  Y _ 

4 

14 

25  Y _ 

8 

3 

31  Y _ 

10 

4 

37  Y _ 

12 

4 

43  Y _ 

14 

2 

49  Y _ 

16 

3 

55  Y _ 

18 

2 

61  L 

41 

65  CL 

68  CL 

1 

72  CL 

2 

76  CL 

3 

80  CL 

4 

84  CL 

5 

88  CL 

6 

92  Y 

41 

0 

0 0 

99  Y A+ 

5 

-4D 

1999 

106  Y A+ 

7 

-3D 

1999 

113  L 

42 

117  Y 

45 

2 

0 -3D 

124  Y 

45 

1 

0 -4D 

131  * 

2 

-7D 

136  : 

2 

100 

141  * 

1 

-5D 

146  : 

1 

100 

151  SP 

1 

2 

2 

DIRECT  POGO  SOFTWARE 
Label  40 

Load  block  18  into  block  0 

Store  block  0 into  block  20 

Set  VP3  to  8 

Set  VP4  to  14 

Set  VP8  to  3 

Set  VPIO  to  4 

Set  VP12  to  4 

Set  VP14  to  2 

Set  VP16  to  3 

Set  VP18  to  2 

Label  41 


Clear  blocks.  0 through  6 


Call  User  Program  41,  where 
nl-n2=n3=0  (Initialize  off- 
line BSFA) 

VP5  = 1999  + VP(VP4) 

VP7  = 1999  + VP(VP3) 

Label  42 

Call  User  Program  45;  nl=2,  n2=0, 
n3=VP(VP3)  (off-line  BSFA) 

Call  User  Program  45;  nl=l, 
n2=0,  n3=VP(VP4) 

Multiply  block  2 by  VP(VP7) 

Divide  block  2 by  100 

Multiply  block  1 by  VP(VP5) 

Divide  block  1 by  TOO 

Calculate  averaged  2-channel- 
double  precision  power  spectrum 
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KEYBOARD  FILE  4 - DIRECT 


157  # 

42 

2D 

163  CH 

1 

2 

169  X< 

1 

173  X> 

3D 

177  Y A+ 

3 

183  X< 

2 

187  X> 

3D 

191  : 

5 

2004D 

196  * 

5 

20 

201  : 

6 

2004D 

206  * 

6 

20 

211  X< 

5 

215  X> 

4D 

219  Y A+ 

4 

225  X< 

6 

229  X> 

4D 

233  Y A+ 

3 

239  Y A+  . 

4 

245  # 

41 

3 

251  L 

43 

255  Y _ 

31 

9 

261  Y _ 

35 

8 

267  Y _ 

36 

5 

273  J 

32 

3 

278  L 

44 

282  Y _ 

3 

2 

288  Y _ 

4 

19 

294  Y _ 

5 

12 

300  Y _ 

6 

10 

306  Y _ 

7 

8 

312  Y _ 

8 

8 

318  Y 

9 

8 

Bll 

POGO  SOFTWARE  (Continued) 

Loop  through  label  42  VP2  times 

Calculate  H(f)  and  y^{f)  from 
the  averaged  power  spectrum  data 

Load  block  1 into  block  0 CH(f)] 

Store  block  0 into  block  VP3 

VP3  = VP3  + 1 

Load  block  2 into  block  0 [Y^(f)] 
Store  block  0 in  block  VP3 
Divide  block  5 by  VP2004 
Multiply  block  5 by  20 
Divide  block  6 by  VP2004 
Multiply  block  6 by  20 
Load  block  5 into  block  0 
Store  block  0 in  block  VP4 
VP4  = VP4  + 1 
Load  block  6 into  block  0 
Store  block  0 into  block  VP4 
VP3  = VPS  + 1 
VP4  = VP4  + 1 

Loop  through  label  41  three  times 
. Label  43 
Set  VP31  to  9 
Set  VP35  to  8 
Set  VP36  to  5 

Jump  to  label  33  in  Keyboard  File  3 

Label  44 

Set  VP3  to  2 

Set  VP4  to  19 

Set  VPS  to  12 

Set  VPS  to  10 

Set  VP7  to  8 

Set  VP8  to  8 

Set  VP9  to  8 
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B12 


KEYBOARD  FILE  4 - DIRECT  POGO  SOFTWARE  (Concluded) 


324  Y _ 

10 

8 

330  Y _ 

11 

15 

336  Y _ 

12 

58 

342  Y _ 

13 

150 

348  J 

50 

5 

353  . 

Set  VPIO  to  8 

Set  VPll  to  15 

Set  VP12  to  58 

Set  VP13  to  150 

Jump  to  lable  50  in  Keyboard 
File  5 (Plot  Routine) 

End 
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KEYBOARD  FILE  5 - POGO  PLOT  ROUTINE 


B13 


1 L 50 


5 

Y 

5838 

3D 

10 

Y 

5821 

6 

15 

Y 

5865 

3 

20 

L 

51 

24 

Y 

26 

4D 

30 

Y 

IF 

4 

5D 

5 

0 

38 

Y 

IF 

4 

6D 

4 

0 

46 

Y 

IF 

4 

7D 

3 

0 

54 

Y 

IF 

4 

8D 

2 

0 

62 

Y 

IF 

4 

9D 

1 

0 

70 

Y 

IF 

4 

TOD 

1 

2 

78 

J 

no 

82 

Y 

5809 

0 

0 

0 

89 

J 

100 

93 

J 

120 

97 

L 

52 

101 

Y A- 

4 

106 

# 

51  IID 

0 

112  L 53 
116  J 130 
120  Y 5814 


124  MS  34  12D 

129  MS  14 

133  Y R 4 


Label  50 

Call  User  Program  5838;  nl=VP3 
(Read  text  buffer  VP3  into  core) 

Call  User  Program  5821;  nl=6 
(Set  up  the  terminal  as  the 
plot  device) 

Call  User  Program  5865;  nl=3 
(Set  log  plots  to  plot  3 decades , 
vertical  scale) 

Label  51 

Set  VP26  to  the  value  of  VP4 

If  VP4  = VPS,  skip  five  steps 

If  VP4  = VPS,  skip  four  steps 

If  VP4  = VP7,  skip  three  steps 

If  VP4  = VP8,  skip  two  steps 

If  VP4  = VP9,  skip  one  step 

If  VP4  > VPIO,  skip  one  step 

Jump  to  label  110 

Call  User  Program  5809;  nl=0, 
n2=0,  n3=0  (Set  up  to  plot  real/ 
magnitude,  linear  horizontal 
scale,  scale  swi tch  at  12 
o'clock) 

Jump  to  label  100 
Jump  to  label  120 
Label  52 
VP4  = VP4  - 1 

Loop  through  label  51  VPll  times 
Label  53 

Jump  to  label  130 

Call  User  Program  5814  (Blank 
the  terminal  screen) 

Set  the  text  record  pointer  to 
VP12 

Read  the  next  message  (plot 
menu) 

Input  VP4  (data  block  number  to 
be  plotted 
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B14 


KEYBOARD  FILE  5 - POGO  PLOT  ROUTINE  (Continued) 


138  Y 

26 

4D 

Set  VP26  to  the  value  of  VP4 

144  Y 

IF 

4 

lOD 

2 

1 

If  VP4  > VPIO,  skip  2 steps 

152  Y 

IF 

30 

0 

1 

0 

If  VP30  = 0,  skip  1 step 

160  Y 

A+ 

26 

260 

30 

VP26  = VP26  + 30 

167  MS 

14 

Read  the  next  message  (Set  the 
scale  and  press  "CONTINUE") 

171  D 

4D 

Display  block  VP4 

175  Y 

5809 

Call  User  Program  5809;  nl=n2= 
n3=default  (Set  up  the  plot 
according  to  the  display 
switches) 

179  J 

100 

Jump  to  label  100 

183  0 

120 

Jump  to  label  120 

187  D 

40 

Display  block  VP4 

191  J 

53 

Jump  to  label  53 

195  L 

100 

Label  100 

199  Y 

5814 

Call  User  Program  5814 

203  Y 

5829 

Call  User  Program  5829  (Plot 
the  full  data  block) 

207  Y 

5815 

40 

Call  User  Program  5815;  nl=VP4 
(Plot  the  data  from  block  VP4) 

212-  Y 

5816 

Call  User  Program  5816  (Plot 
the  defaulted  axis) 

216  Y 

5817 

260 

Call  User  Program  5817;  nl=VP26 
(Label  the  vertical  axis  with 
text  buffer  message  VP4  and  the 
horizontal  axis  with  the 
default  label) 

221  Y 

5808 

975  600 

Call  User  Program  5808;  nl=970, 
n2=600  (Position  the  cursor  to 
y=97.0%  of  the  vertical  plot 
window  and  x=60.0%  of  the 
horizontal  plot  window) 

227  Y 

5819 

1 

Write  text  buffer  message  #1 

232  < 

Subreturn 

235  L 

no 

Label  110 

239  Y 

IF 

4 

100 

2 

1 

If  VP4  > VPIO,  skip  2 steps 

247  Y 

IF 

30 

0 

1 

0 

If  VP30  = 0,  skip  1 step 
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KEYBOARD  FILE  5 - POGO  PLOT  ROUTINE  (Concluded)  BI5 

VP26  = VP26  + 30 


255 

262 

Y A+ 
TP 

26 

4D 

26D 

30 

266 

Y 

5809 

1 

0 

273  J 

100 

277  Y 

5808 

970  1 3D 

283  Y 

5819 

26D 

288  J 

120 

292  Y 

5809 

0 0 

299  J 

130 

303  J 

100 

307  J 

120 

311  TL 

4D 

315  Y 

5809 

0 0 

322  J 

130 

326  J 

100 

330  J 

120 

334  J 

52 

338  L 

120 

342  J 

130 

346  Y 

5820 

350  < 

353  L 

130 

357  Y _ 

19 

1 

363  L 

131 

367  Y A+ 

19 

372  # 

131 

50  0 

378  < 

381  . 

Convert  block  VP4  to  polar 
coordinates 

Call  User  Program  5809;  nl=l, 
n2=0,  n3=0  (Set  up  to  plot 
imaginary/phase,  linear  hori- 
zontal, scale  switch  to  12 
o'clock) 

Jump  to  label  100 

Call  User  Program  5808;  nl=970, 
n2=VP13 

Call  User  Program  5819;  nl=VP26 

Jump  to  label  120 

Call  User  Program  5809;  nl-0, 
n2=0,  n3=0 

Jump  to  label  130 

Jump  to  label  100 

Jump  to  label  120 

Convert  block  VP4  to  log 
polar  coordinates 

Call  User  Program  5809;  nl=0, 
n2=0,  n3=0 

Jump  to  label  130 

Jump  to  label  100 

Jump  to  label  120 

Jump  to  label  52 

Label  120 

Jump  to  label  130 

Call  User  Program  5820  (Hard 
copy  the  terminal) 

Subreturn 

Label  130 

Set  VP19  to  1 

Label  131 

VP19  = VP19  + 1 

Loop  through  label  131  fifty  times 

Subreturn 

End 
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Cl 


APPENDIX  C 

POGO  SOFTWARE  TEXT  RECORDS  (FILE  4a) 


Line  # Contents 


Message  #1 

0 

1 This  is  a 5451 C program,  written  for  POGO  data  analysis.  The  prograra 

2 reads  4 analog  channels  simultaneously  and  writes  them  directly  to  the 

3 disc  via  the  ADC  Throughput.  The  time  domain  data  is  then  called  from 

4 the  disc  and  H(f)s  calculated  re.  Pulser.  H(f)s  across  devices  are 

5 calculated  by  ratios.  28  JANUARY  80 

6 

7 Setup  the  ADC  for  the  desired  frequency  resolution  (Delta  f),  and 

8 input  signal -to-noise  ratio.  Then,  enter  Blocksize.  (INTEGER  LE  1024) 

9 /* 

Message  #2 

10  Enter  a calibration  factor  for  each  channel.  (4  - FLOATING  POINT) 

11  /* 

Message  #3 

12  Enter  the  number  of  data  acquisition  loops.  (INTEGER) 

13  /* 

Message  #4 

14  Enter  Delta  f as  specified  by  the  ADC.  (FLOATING  POINT) 

15  /* 

Message  #7 

20  Start  the  data  tape.  Press  "CONTINUE." 

21  /* 
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Contents 


C2 

Line  # 

Message  #8 

RATIO  POGO  DATA  PLOT  MENU 

24  The  program  output  is  located  in  core  memory  as  follows 

25 

26  BLOCK#  DATA 

27 


28 

3 

Weighted  H(f)  HPOP/LPOP 

29 

4 

Weighted  H(f)  MCC/HPOP 

30 

5 

Weighted  H(f)  MCC/LPOP 

31 

6 

H(f) 

HPOP/LPOP 

32 

7 

Coh(f) 

II 

33 

8 

H(f) 

MCC/HPOP 

34 

9 

Coh(f) 

II 

35 

10 

H(f) 

MCC/LPOP 

36 

11 

Coh(f) 

II 

37 

12 

H(f) 

LPOP/Pulser 

38 

13 

Coh(f) 

11 

39 

14 

H(f) 

HPOP/Pulser 

40 

15 

Coh(f) 

II 

41 

16 

H(f) 

MCC/Pulser 

42 

17 

Coh(f) 

II 

43 

18 

PSD 

Pul ser 

44 

19 

PSD 

LPOP 

45 

20 

PSD 

HPOP 

46 

21 

PSD 

MCC 

47 

22 

Re  XPSD  LPOP/Pulser 

48 

23 

Im  XPSD 

II 

49 

24 

Re  XPSD  HPOP/Pulser 

50 

25 

Im  XPSD 

II 

51 

26 

Re  XPSD  MCC/Pulser 

52 

27 

Im  XPSD 

II 

53 

54  Please  enter  the  Block#  you  wish  to  plot. 

55  /* 

Message  #9 

56  Set  the  scale  and  press  "CONTINUE." 

57  /* 
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C3 


Line  # Contents 

Message  #10 

DIRECT  POGO  DATA  PLOT  MENU 

58  The  program  output  is  located  in  core  as  follows: 

59 


60 

BLOCK# 

DATA 

61 

62 

5 

Weighted  H(f)  HPOP/LPOP 

63 

6 

Weighted  H(f)  MCC/HPOP 

64 

7 

Weighted  H(f)  MCC/LPOP 

65 

8 

H(f)  HPOP/LPOP 

66 

9 

Coh(f)  " 

67 

10 

H(f)  MCC/HPOP 

68 

11 

Coh(f) 

69 

12 

H(f)  MCC/LPOP 

70 

13 

Coh(f) 

71 

14 

Re  XPSD  HPOP/LPOP 

72 

15 

Im  XPSD 

73 

16 

Re  XPSD  MCC/HPOP 

74 

17 

Im  XPSD 

75 

18 

Re  XPSD  MCC/LPOP 

76 

19 

Im  XPSD 

77 

78  Please  enter  the  block  # you  wish  to  plot. 

79  /* 


Message  #11 

80  Set  the  scale  and  press  "CONTINUE". 

81  /* 


Message  #12 

82  Do  you  want:  Coherence  blanking  (enter  2) 

83  Pulser  blanking  (enter  0)? 

84  /* 

Message  #13 

85  Enter  the  minimum  blanking  value  (floating  point). 

86  /* 
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D1 


APPENDIX  D 

P060  SOFTWARE  TEXT  BUFFERS  (FILE  4b) 


Text  Buffer  #1 
RATIO  P060  PLOT  LABELS 

01 

(TEST  ID) 

/* 

13 

HPOP/LPOP  Ratio  H(f) 

Pulser  Blanked 

/* 

04 

MCC/HPOP  Ratio  H(f) 

Pulser  Blanked 

/* 

05 

MCC/LPOP  Ratio  H(f) ' 

Pulser  Blanked 

/* 

06 

HPOP/LPOP  Ratio  H(f)  - Polar 
/*  ‘ 

07 

HPOP/LPOP  Ratio  Coh(f) 

/* 

08 

MCC/HPOP  Ratio  H(f)  - Polar 

/* 

09 

MCC/HPOP  Ratio  Coh(f) 

/* 

10 

MCC/LPOP  Ratio  H(f)  - Polar 

/* 

11 

MCC/LPOP  Ratio  Coh(f) 

/* 
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D2 

Text  Buffer  #1  (Continued) 

12 

LPOP/Pulser  H(f)  - Polar 

/* 

13 

LPOP/Pulser  Coh{f) 

/* 

14 

HPOP/Pulser  H(f)  - Polar 

/* 

15 

HPOP/Pulser  Coh(f) 

/* 

16 

MCC/Pulser  H(f)  - Polar 

/* 

17 

MCC/Pulser  Coh(f) 

/* 

18 

Pulser  PSD 

/* 

19 

LPOP  PSD 

/* 

20 

HPOP  PSD 

/* 

21 

MCC  PSD 

/* 

22 

LPOP/Pulser  Real  XPSD 

/* 

23 

LPOP/Pulser  Imag  XPSD 

/* 

24 

HPOP/Pulser  Real  XPSD 

/* 
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D3 

Text  Buffer  #1  (Concluded) 

25 

HPOP/Pulser  Imag  XPSD 

/* 

26 

MCC/Pulser  Real  XPSD 

/* 

Z1 

MCC/Pulser  Imag  XPSD 

/* 

33 

HPOP/LPOP  Ratio  H(f) 

Coherence  Blanked 

/* 

34 

MCC/HPOP  Ratio  H(f) 

Coherence  Blanked 

/* 

35 

MCC/LPOP  Ratio  H(f) 

Coherence  Blanked 

/* 


Text  Buffer  #2 
DIRECT  POGO  PLOT  LABELS 

01 

(TEST  ID) 

/* 

05 

HPOP/LPOP  Direct  H(f) 

Pulser  Blanked 

/* 

06 

MCC/HPOP  Direct  H(f) 

Pulser  Blanked 

/* 

07 

MCC/LPOP  Direct  H(f) 

Pulser  Blanked 

/* 
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D4 

Text  Buffer  #2  (Continued) 

08 

HPOP/LPOP  Direct  H(f)  - Polar 

/* 

09 

HPOP/LPOP  Direct  Coh(f) 

/* 

10 

MCC/HPOP  Direct  H(f)  - Polar 
/* 

n 

MCC/HPOP  Direct  Coh(f) 

/* 

12 

MCC/LPOP  Direct  H(f)  - Polar 
/* 

13 

MCC/LPOP  Direct  Coh(f) 

/* 

14 

HPOP/LPOP  Direct  Real  XPSD 

/* 

15 

HPOP/LPOP  Direct  Imag  XPSD 

/* 

16 

MCC/HPOP  Direct  Real  XPSD 

/* 

17 

MCC/HPOP  Direct  Imag  XPSD 

/* 

18 

MCC/LPOP  Direct  Real  XPSD 
/* 

19 

MCC/LPOP  Direct  Imag  XPSD 

/* 

35 

HPOP/LPOP  Direct  H(f) 

Coherence  Blanked 

/* 
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Text  Buffer  #2  (Concluded) 


D5 


36 

MCC/HPOP  Direct  H(f) 
Coherence  Blanked 

/* 

37 

MCC/LPOP  Direct  H(f) 
Coherence  Blanked 

/* 
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INTRODUCTORY  SUMMARY 


This  report  consists  of  a series  of  charts,  equations,  and  the  computer  prc^am  listing 
for  determining  the  frequencies  associated  with  the  roUer  bearings  of  the  SSME  high 
pressure  fuel  turbopump  (HPFTP)  and  high  pressure  oxygen  turbopump  (HPOTP).  The 
frequencies  are  in  terms  of  speed  ratio  normalized  by  the  rotational  speed  of  the 
shaft.  Speed  ratio  was  calculated  for  the  cage,  ball  spin,  outer  race,  and  inner  race 
for  the  different  bearing  parameters  of  the  SSME  HPFTP  and  HPOTP. 

Several  methods  of  data  presentation  are  provided  for  ease  in  interpretation.  This 
includes  charts  for  the  charge  in  contact  angle  for  both  the  inner  and  outer  race, 
constant  inner  race  contact  angle  versus  outer  contact  ai^le,  and  charts  of  inner 
contact  angle  equal  to  outer  contact  angle,  which  is  only  applicable  for  the  ease  of 
slow  rotational  speed  and/or  applied  radial  load  of  large  magnitude. 

The  data  from  two  static  firing  tests  are  also  included  where  the  cage  and  outer  race 
frequencies  were  noted  on  internal  strain  gages  mounted  directly  on  the  bearing 
holder.  Several  attempts  were  made  usii^  different  analysis  techniques  (including 
adaptive  filtering,  cross  PSDs,  etc.)  to  isolate  these  frequencies  in  the  external 
accelerometer  data.  To  date  these  efforts  have  not  been  successfuL 
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TECHNICAL  DISCUSSION 


The  speed  ratio  for  the  cage,  baE  spin,  and  outer  and  inner  race  are  defined  as  foUows: 


1 


Cf^e  Speed  Ratio 


1 - y cos 
1 + cos 


Outer  Race  Speed  Ratio 


Njj  (1-r  cos  a.) 

1 + cos  (a.- a) 
1 o 


Inner  Race  Speed  Ratio 


N|^  cos  y*  cos  Q. 

''  1 + cos  (a.-a J 

I o 


where 


and 


Ball  Spin  Ratio  = < I (cos  + tan  yS  sin  (1  + y*  cos 


-1 


+ (cos  + tan/S  sin  cn.)  (1  - y*  cos 


-1 


.-1 


y cos  jSi 


sin 

tan^  - cos  + y' 

= outer  raceway  contact  ai^le 
or.  = inner  raceway  contact  ar^le 

r=  D/d^ 

D = ball  diameter 

= pitch  diameter 
m 

Nj^  = number  of  balls 


1 


Tedric  A.  Harris.  Rolling  Bearing  Analysis,  John  Wiley  & Sons,  New  York,  1966. 
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The  bearing  parameters  used  for  this  study  are  as  follows: 

High  Pressure  Fuel  Turbopump  (HPFTP) 

Pump  and  turbine  bearii^s: 

Ball  diameter  = 0.343750  inch 
Pitch  diameter  = 2.34  inches 
Number  of  balls  = 14 

Thrust  bearings: 

Ball  diameter  = 0.656250  inch 
Pitch  diameter  = 3.84  inches 
Number  of  balls  = 13 

High  Pressure  Oxygen  Turbopump  (HPOTP) 

Turbine  end  bearings: 

Ban  diameter  = 0.50  inch 

Pitch  diameters  = 3.196,  3.17,  3.20,  and  3.19  inches 
Number  of  balls  = 13 

Pump  end  bearings: 

Ball  diameter  = 0.4375  inch 
Pitch  diameter  = 2.56  inches 
Number  of  balls  = 13 
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SPEED  RATIO  OF 

HPFTP  PUMP  AND  TURBINE  BEARINGS 

Ball  Diameter  = 0.343750  inch 
Pitch  Diameter  = 2.34  inches 
Number  of  Balls  = 14 
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HPFTP  PUMP  A TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  D1AMETER*0. 34375  PITCH  DI AMETER*2 . 3400  NUMBER  BALLS*  14 
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SPEED  RAHO  OF 
HPFTP  THRDST  BEARINGS 

Ball  Diameter  = 0.656250  inch 
Pitch  Diameter  = 3.84  inches 
Number  of  Balls  = 13 
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HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERS0. 65625  PITCH  DlAMETER^S . 8400  NUMBER  BALLS*  13 
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OUTER  RACE  CONTACT  ANCLE  r DEGREES I 


SPEED  RATIO  OF 
HPOTP  TURBINE  BEARINGS 

Ball  Diameter  ==  0.50  inch 
Pitch  Diameter  s 3.196  inches 
Number  of  BaUs  - 13 
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SPEED  RATIO  OF 
HPOTP  PUMP  END  BEARINGS 

Ball  Diameter  = 0.4375  inch 
Pitch  Diameter  = 2.56  inches 
Number  of  Balls  = 13 
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OUTER  RACE 


HPOTP  PUMP  END  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAnET£R»0.43750  PITCH  DIAMETERa2 . 5600  NUMBER  BALLS=  IS 


— Z2UJOC  (T-<UUi  UOZI-<Ol-  <ZU^UJ 


to 


24 


OUTER  RACE  CONTACT  ANCLE  (DEGREES) 


HPOTP  punp  END  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETER*0.43750  PITCH  01 AnETERs2 .5600  NUMBER  BALLS*  13 


>-<ZZUiQ£  a<UUJ  UOZ^-<U»-  <ZU-JUi 


25 


OUTER  RACE  CONTACT  ANGLE  (DEGREES) 


HPOTP  PUMP  END  BEARINGS  PARAMETRIC  STUDY 
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OUTER  RACE  CONTACT  ANGLE  f DEGREES) 


SPEED  RATIO  OF 

HPFTP  PUMP  AND  TURBINE  BEARINGS 
3-,  10-,  20-,  AND  30'DEGREE  a. 
(Inner  Race  Contact  Angle) 

Ball  Diameter  = 0.343750  inch 
Pitch  Diameter  = 2.34  inches 
Number  of  Balls  = 14 
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SPEED  RATIO  OF 
HPFTP  THRUST  BEARINGS 
3-,  10-,  20-,  AND  30-DEGREE  «. 
(Inner  Race  CkHitaet  Angle) 

Ball  Diameter  = 0.656250  inch 
Pitch  Diameter  = 3.84  inches 
Number  of  Balls'  = 13 


44 


HPFTP  THRUST  BEARINGS  PARANETRIC  STUDY 
BALL  DlAMETER:t0.65625  PITCH  DIANETERa5.8400  NUhBER  BALLS«  13 

INNER  RACE  CONTACT  ANGLE  * 3 


U < U LU  </)  a.  UJ  UJ  Q OC  < I- 


OUTER  RACE  CONTACT  ANCLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAMETERa0.65625  PITCH  OIAriETER*3.8400  NUMBER  BALLS®  13 

INNER  RACE  CONTACT  ANCLE  ® 10 


U < O Ui  to  0.  LU  UJ  Q OK  >—  O 


46 


OUTER  RACE  CONTACT  ANGLE  1 DEGREES J 


HPFTP  THRUST  BEARINGS  PARAHETRIC  STUDY 
BALL  DIAriETER=0.65625  PITCH  0IAMETER*3 .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANCLE  * 20 


OUTER  RACE  CONTACT  ANGLE  fOEGREES) 


HPFTP  THRUST  BEARINGS  PARAflETRIC  STUDY 
BALL  DIAMETER*0.65625  PITCH  DI AnET£R*3 .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANGLE  « 30 


48 


OUTER  RACE  CONTACT  ANGLE  f DEGREES) 


HPFTP  thrust  bearings  PARAMETRIC  STUDY 
BALL  01 AMETERsB. 65625  PITCH  01 AMETERaS .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANGLE  = 3 


OUTER  RACE  CONTACT  ANGLE  (DEGREES 1 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERs0.65625  PITCH  DIAMETERsS .8400  NUMBER  BALLSs  13 

INNER  RACE  CONTACT  ANCLE  « 10 


50 


OUTER  RAGE  CONTACT  ANGLE  I DECREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DlAMETERsB. 65625  PITCH  0IAMETER«3 .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANGLE  * 20 


51 


OUTER  RACE  CONTACT  ANGLE  (DEGREES J 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAMETER«0.65825  PITCH  DIAMETER*3 .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANGLE  * 30 


CD  < a.  *-  2 q:  < o 


52 


OUTER  RACE  CONTACT  ANGLE  ( DEGREES J 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERs0.65625  PITCH  0 I AMETERs? . 8400  NUMBER  BALLSs  IS 

INNER  RACE  CONTACT  ANGLE  » S 


00 

w 

CM 

GO 

CO 

CM 

ca 

63 

C9 

CD 

0> 

03 

CD 

03 

CM 

s 

a 

N 

CO 

in 

tr 

to 

to 

w 

CO 

to 

w 

in 

in 

ID 

ID 

ID 

ID 

o ID  I-  ui  q: 

a < 

U LU 

(/)Ql 

LU  UJ  o 

o:  < 

s 

ID 

ID 


ID 

to 


lO 


ID 


M 

ID 


G2 


ID 


53 


OUTER  RACE  CONTACT  ANGLE  IDEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAf1ET£Rs0.65625  PITCH  DIAf1ETER«3.8400  NUMBER  BALLS*  15 

INNER  RACE  CONTACT  ANGLE  * t0 


54 


OUTER  RACE  CONTACT  ANCLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERs0.65625  PITCH  OIAnET£Rs3.8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANGLE  * 20 


uioc  a<uuj  (/)Q.uiujQ  o:<H 


OUTER  RAGE  CONTACT  ANGLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETER=0.65625  PITCH  DIAMETER*3 .8400  NUMBER  BALLS*  13 

INNER  RACE  CONTACT  ANCLE  « 30 


56 


OUTER  RACE  CONTACT  ANGLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAnETRlC  STUDY 
BALL  OIAriETERs0.e5625  PITCH  01 AMETERaS .8400  NUMBER  BALLS®  13 

INNER  RACE  CONTACT  ANGLE  * 3 


57 


OUTER  RACE  CONTACT  ANCLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  01 AMETER«0. 65625  PITCH  DIAMETER*5.8400  NUMBER  BALLS#  13 

INNER  RACE  CONTACT  ANCLE  « 10 


58 


OUTER  RACE  CONTACT  ANCLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAttETRIC  STUDY 
BALL  DIAMETER«0.65625  PITCH  01 ANETERsS .8400  NUMBER  BALLS«  tS 

INNER  RACE  CONTACT  ANGLE  * 20 


zzuicc  ce<uui  (noLUiuia 


OUTER  RACE  CONTACT  ANGLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DlAMETERsB. 65625  PITCH  OIAMCTER»3.8400  NUMBER  BALLSs  13 

INNER  RACE  CONTACT  ANGLE  = 30 


60 


OUTER  RACE  CONTACT  ANGLE  (DEGREES) 


SPEED  RATIO  OF 
HPOTP  TURBINE  BEARINGS 
3-,  10-,  20-,  AND  30-DEGREE  «. 
(Inner  Race  Ck>ntact  Angle) 

Ball  Diameter  = 0.50  inch 
Pitch  Diameter  = 3.196  inches 
Number  of  Balls  - 13 
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HPOTP  PUMP  END  BEARINGS 
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(Emer  Race  Contact  Angle) 

Ball  Diameter  = 0.4375  inch 
Pitch  Diameter  = 2.56  inches 
Number  of  Balls  = 13 
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OUTER  RACE 


SPEED  RATIO  OF 
DIFFERENT  PITCH  DIAMETERS 
OF  HPOTP  TURBINE  BEARINGS 

Ball  Diameter  - 0.50  inch 
Pitch  Diameter  = 3.17,  3.19,  and  3.20  inches 
Number  of  Balls  = 13 

Inner  Contact  Angle  = Outer  Contact  Angle 
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HPOTP  TURBINE  BEARINGS  PARAflETRIC  STUDY 
BALL  OIANETERS0.50000  PITCH  DIAnETER=:3 . 1700  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 


CONTACT  ANGLE  (DECREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETER=0.50000  PITCH  DIAMETER*? . 1 900  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANGLE  (DEGREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETER =0.50000  PITCH  DIAMETER«5.2000  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANGLE  (DEGREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DlAMeTERaB. 50000  PITCH  DIAMETERaS. 1700  NUMBER  BALLSa  13 
INNER  CONTACT  ANCLE  a OUTER  CONTACT  ANCLE 
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CONTACT  ANGLE  IDECREESl 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAMETER«0.50000  PITCH  D I AMETER*S . 1900  NUMBER  BALLS*  \% 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 


CONTACT  ANCLE  (DECREES) 


HPOTP  TURBINE  BEARINGS  PARANETRIC  STUDY 
BALL  DlANETERaB. 50000  PITCH  0IAMETER«3 .2000  NUMBER  BALLS*  15 
INNER  CONTACT  ANCLE  * OUTER  CONTACT  ANCLE 
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CONTACT  ANCLE  ( DEGREES 1 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BAU.  DIAMETERS0.50000  PITCH  DI AMETERsS. 1 700  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DJAMETER=0. 50000  PITCH  DIAMETER«3. 1900  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANCLE  C DECREES) 


HPOTP  TURBINE  BEARINGS  PARAOETRIC  STUDY 
BALL  DIAf1ETERs0.50000  PITCH  OIAMETERaJ . 2000  NUMBER  BALLSa  13 
INNER  CONTACT  ANGLE  « OUTER  CONTACT  ANGLE 
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CONTACT  ANCLE  (DEGREES J 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETER*0.50000  PITCH  OIAMETERaS . 1700  NUMBER  BALLS»  13 
INNER  CONTACT  ANCLE  = OUTER  CONTACT  ANGLE 


CONTACT  ANGLE  1 DEGREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAMETER*0.50000  PITCH  DIAMETER«3. 1900  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  « OUTER  CONTACT  ANCLE 


CONTACT  ANGLE  (DECREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAf1ETERs0.50000  PITCH  DIAMETERsJ .2000  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  » OUTER  CONTACT  ANGLE 


107 


CONTACT  ANCLE  (DEGREES) 


SPEED  RATIO  WITH  DYNAMIC  EFFECTS  NEGLECTED 


Charts  of  this  section  are  for  simple  rolling  motion,  which  is  only  applicable  for  the 
case  of  slow  rotational  speed  and/or  applied  radial  load  of  large  magnitude.  The 
dynamic  effects  were  neglected:  i.e.,  outer  contact  angle  equal  to  inner  contact  angle. 
The  equations  for  this  case  are  as  follows: 

C£^e  speed  ratio  = 1/2  (1  - D/d^^  cos  0) 

Outer  race  speed  ratio  = Nj^/2  (1  - D/d^^  cos  <!>) 

Inner  race  speed  ratio  = Nj^/2  (1  + D/d^  cos  <ih 
Ban  spin  ratio  = 1/2  d^^/D  j"!  - (D/dj^)^  cos^ 

D = BaU  diameter 

d_  = Pitch  diameter 
m 

"b=  Number  of  balls 
0 = Contact  angle 
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HPFTP  punp  & TURBINE  BEARINGS  PARANETRIC  STUDY 
BALL  OIAMETERsB. 34375  PITCH  DIAnETER*2.3400  NUMBER  BALLS*  H 
INNER  CONTACT  ANGLE  = OUTER  CONTACT  ANGLE 
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contact  angle  (DEGREES 1 


HPF7P  PUMP  & TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERsB. 34575  PITCH  DIAMETER«2 . 3400  NUfIBER  BALLS»  14 
INNER  CONTACT  ANGLE  « OUTER  CONTACT  ANGLE 
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CONTACT  ANGLE  I DECREES) 


HPFTP  PUMP  & TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  OIAMETERS0.34375  PITCH  OIAMETER»2.3400  NUMBER  BALLS#  14 
INNER  CONTACT  ANGLE  « OUTER  CONTACT  ANCLE 
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CONTACT  ANGLE  (DECREES) 


HPFTP  PUMP  & TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DlAMETERaB. 34375  PITCH  DIAMETER*2 * 3400  NUMBER  BALLS®  14 
INNER  CONTACT  ANGLE  » OUTER  CONTACT  ANCLE 
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CONTACT  ANGLE  I DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAf1ETER*0.65e25  PITCH  DI AMETER*3 .8400  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  = OUTER  CONTACT  ANCLE 
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CONTACT  ANCLE  (DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERsB. 65625  PITCH  DIAMETER*S .8400  NUMBER  BALLS*  IS 
INNER  CONTACT  ANGLE  » OUTER  CONTACT  ANCLE 


CQ  •<  j -a  U)  Ql  Z Qi  < i-  *-•  O 


114 


CONTACT  ANGLE  f DEGREES) 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAnETER«0.65625  PITCH  DI AMETERsJ .8400  NUMBER  BALLS«  13 
INNER  CONTACT  ANCLE  « OUTER  CONTACT  ANGLE 
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CONTACT  ANGLE  I DECREES 1 


HPFTP  THRUST  BEARINGS  PARAMETRIC  STUDY 
BALL  01 AMETERS0. 65625  PITCH  OIAMETER-3.8A00  NUMBER  BALLS*  IS 
INNER  CONTACT  ANGLE  « OUTER  CONTACT  ANCLE 
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CONTACT  ANGLE  (DEGREES) 


hpotp  puhp  end  beapings  parametric  study 

BALL  DIAflETERsB. 43750  PITCH  DIAMETER*2.56B0  NUMBER  BALLS*  13 
INNER  CONTACT  ANCLE  » OUTER  CONTACT  ANCLE 
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CONTACT  ANCLE  I DECREES) 


HPOTP  PUMP  END  BEARINGS  PARAMETRIC  STUDY 
BALL  01 AMETERsB. 43750  PITCH  DIAMETER«2 .5600  NUMBER  BALLSa  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANCLE 


CONTACT  ANGLE  (DECREES) 


HPOTP  PUMP  END  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERaB.43750  PITCH  DIAMETERs2.5600  NUMBER  BALLS*  13 
INNER  CONTACT  ANCLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANGLE  (DECREES) 


HPOTP  PUMP  END  BEARINGS  PARAnETRIC  STUDY 
BALL  DIAhETERs0.43750  PITCH  DI AriETER»2 .5600  NUMBER  BALLS*  13 
INNER  CONTACT  ANCLE  * OUTER  CONTACT  ANGLE 


CONTACT  ANCLE  « DEGREES! 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAf1ETER*0.50000  PITCH  DiAMETERsB. I960  NUMBER  BALLSs  13 
INNER  CONTACT  ANCLE  * OUTER  CONTACT  ANCLE 
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contact  angle  lOECREES) 


HPOTP  TURBINE  BEARINGS  PARANETBIC  STUDY 
BALL  DIAnET£Rs0.50000  PITCH  DIAriETER«3. 1960  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANCLE  fOECREES) 


HPOTP  TURBINE  BEARINGS  PARAMETRIC  STUDY 
BALL  DIAMETERS0.50000  PITCH  DIAMETER*3. 1960  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANGLE 
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CONTACT  ANCLE  (DECREES) 


HPOTP  TURBINE  BEARINGS  PARANETRIC  STUDY 
BALL  OIAHETERa0.50000  PITCH  DIAMETER»3. 1060  NUMBER  BALLS*  13 
INNER  CONTACT  ANGLE  * OUTER  CONTACT  ANCLE 
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CONTACT  ANGLE  f DEGREES) 


COMPUTER  PROGRAM  LISTING 
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420  A«="HP0TP  OBBOBD  BEARINGS  PARAhETRIC  STUDY" 

420  A^=’*HP0TP  TURBINE  BEARINGS  PARAMETRIC  STUDY 

420  A«  = *‘HP0TP  TURBINE  BEARINGS  PARAMETRIC  STUDY" 
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HPOTP  TURBINE  BEARINGS  PARAflETRIC  STUDY 
BALL  OIANETER»0.50000  PITCH  DIAHETER=S . I960  NUMBER  BALLS*  IS 

INNER  RACE  CONTACT  ANGLE  * 3 
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1 15  20  25  30  35 

OUTER  RACE  CONTACT  ANGLE  f DEGREES I 
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INTERNALLY  MEASURED  BEARING  FREQUENCIES 
OF  TWO  STATIC  TEST  FIRINGS 


This  section  includes  selected  data  in  the  form  of  working  papers  for  two  static  firing 
tests,  750-215  and  750-234.  The  two  major  points  of  interest  are  the  cage  frequency  on 
test  750-215  and  the  outer  race  frequency  on  test  750-234.  Of  special  interest  is  the 
beating  of  the  bearing  frequencies  with  the  synchronous  frequency  of  the  pump.  The 
data  in  this  section  is  only  included  to  illustrate  or  demonstrate  the  applicability  of  the 
charts  of  the  previous  sections.  A comparison  of  the  measured  frequency  on  test  750- 
234  and  the  applicable  chart  of  the  HPOTP  outer  race  speed  is  also  shown.  Several 
points  of  possible  contact  angles  for  the  measured  speed  ratio  are  shown  on  the  chart. 
Also  included  is  a plot  of  the  outer  race  frequency  for  different  power  levels,  including 
a slow  ramp  down.  For  this  test,  the  contact  ai^le  does  not  appear  to  be  a function  of 
power  level.  When  additional  data  from  normal  bearing  operation,  and  especially  from 
bearings  that  have  experienced  damage,  becomes  available,  it  is  recommended  that 
detailed  analysis  be  performed  and  included  in  a report. 
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1.0  INTRODUCTORY  SUMMARY 


This  report  documents  the  results  of  a statistical  evaluation  of  RMS  (root-mean- 
square)  vibration  levels  measured  during  flight  of  the  Space  Shuttle.  The  data  was 
recorded  on-board  the  Orbiter,  from  accelerometers  located  on  the  Space  Shuttle  Main 
Engine  (SSME)  high  pressure  oxidizer  turbopumps  (HPOTP)  and  high  pressure  fuel 
turbopumps  (HPFTP).  It  should  be  noted  that  the  data  recorded  on-board  the  Shuttle 
vehicle  has  a 50-800  Hz  bandpass  filter  and  was  analyzed  using  a 50-1000  Hz  band- 
width. Ground  test  data  are  acquired  over  a wider  frequency  range.  Therefore,  any 
direct  comparison  of  the  composite  vibration  levels  of  the  flight  data  to  measure- 
ments recorded  at  the  ground  test  stands  without  the  bandpass  filter  will  require 
the  application  of  a small  correction  factor.  However,  the  synchronous  vibration 
levels  will  compare  directly. 

The  evaluation  included  the  composite  and  synchronous  vibration  levels  at  65%,  100% 
and  104%  power  levels.  Results  are  presented  in  the  form  of  summary  sheets  listing 
the  mean  and  standard  deviation  for  each  individual  measurement;  plots  of  the 
cumulative  distribution  with  a Gamma  distribution  overlay,  and  a plot  of  the  density 
function.  Also  presented  is  the  ratio  of  the  synchronous  to  composite  vibration 
levels  for  the  high  pressure  fuel  and  oxidizer  turbopumps  during  100%  and  104%  power 
levels. 

Plots  are  presented  of  the  Normal,  Rayleigh,  and  Gamma  Probability  Distribution 
functions  for  comparison  with  over  a thousand  measured  ground  test  data  points.  The 
data  represents  the  oxidizer  composite  and  synchronous  vibration  levels  from  ground 
static  firings  conducted  on  the  Al,  A2,  and  A3  test  stands.  An  examination  of  the 
plots  demonstrates  the  Gamma  distribution  function  should  provide  reasonable 
approximations  to  the  flight  data.  The  application  of  a classical  distribution  is 
desirable  for  data  characterization  since  this  permits  continuous  statistical 
definition  and  manipulation  from  discrete  flight  measurement  observations.  Sketches 
are  included  of  the  accelerometer  measurement  locations  and  the  power  profile  plots 
for  each  flight. 

All  of  the  plots  contained  in  this  report  are  the  automatic  output  of  the  menu 
driven  MSEC  Diagnostic  Data  Base  Program.  Details  of  the  program  are  documented  in 
Reference  1. 
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2.0  TECHNICAL  DISCUSSION 

2.1  Space  Shuttle  Flights 

Table  I summarizes  the  19  flights  of  the  Space  Shuttle  as  of  August  1985.^  The  data 
analysis  is  based  upon  18  flights,  since  only  the  mean  square  vibration  levels 
sampled  at  100  msec  intervals  were  recorded  on  STS-9.  This  type  of  data  is  not 
compatible  with  the  analysis  methods  used  for  this  report.  Included  in  Table  I are 
the  launch  date.  Space  Shuttle  Vehicle  identification,  maximum  power  level  during 
flight,  engine  serial  number,  HPOTP  serial  number,  and  HPFTP  serial  numbers. 

2.2  Measurement  Location 

The  measurement  locations  on  the  HPOTP  are  shown  in  Figure  1 with  the  block  arrange- 
ments for  the  135®-1  and  135®-2  locations  shown  in  Figure  2.  A few  of  the  early 
flights  had  measurements  located  at  the  180®  position,  while  the  present  flight 
measurement  program  is  standardized  to  3 measurements  at  the  45®,  135®-1  and  135®-2 
positions.  No  measurements  during  flight  are  located  on  the  High  Pressure  Oxidizer 
Turbine  (HPOT)  or  the  low  pressure  oxidizer  turbopumps  (LPOTP). 

Figure  3 shows  the  measurement  locations  for  the  HPFTP.  The  measurement  location  of 
early  flights  was  also  somewhat  different  for  the  fuel  turbopump  than  the  present 
standardized  plan  of  3 measurements  at  the  0®,  174®  and  186®  positions.  Data 
recorded  at  180®  position  on  the  early  flights  is  listed  at  the  186®  location  for 
the  purpose  of  storage  in  the  computer  data  bank.  In  addition,  a few  flights 
included  data  at  the  90®  position  rather  than  the  present  0®  location.  No  measure- 
ments are  located  on  the  high  pressure  fuel  turbine  (HPFT)  end  or  the  low  pressure 
fuel  turbopump  (LPFTP)  during  flight.  All  measurements  are  in  the  radial  direction 
for  both  the  HPOTP  and  HPFTP.  The  relationship  of  the  accelerometer  measurement 
location  plane  to  the  SSME  powerhead  components  is  shown  in  Figure  4. 


^ Appendix  A includes  results  from  flights  27  (51-1)  and  28  (51-J)  which  were 
not  available  at  the  time  the  statistical  analysis  was  performed.  The  vibration 
levels  of  both  flights  were  nominal  and  will  not  provide  a significant  change  to  the 
statistical  values  of  this  report.  As  additional  flight  data  becomes  available, 
this  report  will  be  periodically  updated. 
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Figure  1.  High  Pressure  Oxidizer  Turbopump  Accelerometer  Locations,  Flight 
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Figure  2.  HPOTP  Accelerometer  Block  135®-!  and  135®-2,  Flight 
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Figure  4.  Accelerometer  Location  Plane,  Flight 


2.3  Power  Profiles 

The  power  profiles  for  each  of  the  18  flights  are  shown  in  Figures  5 through  8.  The 
profiles  for  all  ground  tests  and  flights  are  stored  in  the  MSFC  Diagnostic  Data 
Bank  for  retrieval  and  plotting.  It  should  be  noted  that  two  profiles  are  shown  for 
Mission  STS  51-F  indicating  the  early  cutoff  of  engine  number  1 during  this  flight. 

2.4  Mean  and  Standard  Deviation  of  Vibration  Measurements 
Given  a data  sample  of  size  N the  mean  and  variance  are  defined  by 


X =»  1/N  ^ Xj 

N 


(Sample  Mean) 


(2.4-1) 


“ 1/(N  - 1)  XI  (^i  - x)^  (Sample  Variance) 

N 


(2.4-2) 


The  standard  deviation  is  the  positive  square  root  of  the  variance  and  defined  by 


SD 


X (Xi  - x)2 
N 


(Sample  Standard  Deviation) 


(2.4-3) 


2.4.1  Combined  Data  Groups 

Each  of  the  above  statistics  was  first  calculated  for  each  measurement  location  and 
then  combined  as  necessary  for  printout  Or  plotting.  The  following  is  a discussion 
of  procedures  for  combining  data  depending  on  whether  the  sums  and  the  sum  of 
squared  terms  or  only  the  mean  and  standard  deviation  are  stored. 


Assume  we  wish  to  combine  the  vibration  levels  measured  on  the  LOX  PBP  135®-!  with 
the  levels  measured  on  LOX  PBP  135®-2,  where  the  mean  and  standard  deviation  were 
calculated  for  each  measurement.  If  the  sums  and  sum  of  squared  terms  are  stored 
the  combined  mean  and  standard  deviation  are  as  follows. 


H (Xi)l  + X(Xi>2 

Nj  Nj 

Xi2  = — (2.4-4) 

Nj  + Nj 
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Figure  6.  Flight  Power  Profile,  STS-07  to  STS-41D 
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where  and  (Xj)j  are  the  number  of  measurements  and  vibration  value  respectively 
of  the  LOX  PBP  135°-1  location,  Nj,  (Xj)2  the  values  of  the  LOX  PBP  135°-2  location 
and  Xj^2  is  the  combined  mean  vibration  value. 


Equation  (2.4-3)  for  calculation  of  the  standard  deviation  is  modified  to  a format 
where  the  only  requirement  is  storage  of  the  sums.  Using  Equations  (2.4-1)  and 
(2.4-3) 


SD 


(2.4-5) 


(SD)2  (N-1) 


E 

N 


(Xi  - X)2 


= (xj®  - 2X|  X + x^) 

N 


= E (Xi)^  - E 2 Xj  X + X (x)^ 
•N 

Since 


E 2 Xj  = 2 N X 

N ■ 


and 


E Cx)2  = N(x)2 

N 


(SD)2  (N-1)  = 23  (Xi)2  - 2N(x)2  + N(x)2 

N 


= E (Xi)*  - N(x)2 

N 


and 


19 


SD 


(2.4-6) 


(Xi)2-  N(x)2 

N-1 


or 


SD 


(2.4-7) 


As  can  be  easily  seen  from  examination  of  Equation  (2.4-7),  only  the  number  of  data 
points  and  the  sum  of  Xj  and  (xj)^  values  are  required  for  calculation  of  the 
standard  deviation.  Most  hand  calculators  and  computer  programs  use  this  method 
rather  than  Equations  (2.4-1)  and  (2.4-3).  Equations  (2.4-1)  and  (2.4-3)  require 
the  calculation  of  the  mean  for  the  complete  data  set  before  the  sum  of  the  differ- 
ence squared  can  be  calculated. 


The  standard  deviation  for  the  combined  LOX  PBP  135°-1  and  LOX  PBP  135°-2,  assuming 
the  sum  of  the  Xj  and  (xj)*  terms  are  available  is 


(2.4-8) 


However,  if  the  sums  are  not  available  the  combined  mean  and  standard  deviation  can 
be  calculated  from  the  individual  number  of  data  points,  means  and  standard  devia- 
tions as  follows. 

From  Equation  (2.4-1) 

2 X}  = N X (2.4-9) 

N 

and  with  Equation  (2.4-4)  the  combined  mean  of  the  LOX  PBP  135°-1  and  LOX  PBP  135°-2 
is 
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(2.4-10) 


^ Ni  Xj  + Nj  X2 

12  “ — 

Ni  + Nj 

where  the  subscripts  1 and  2 of  N and  x indicate  the  number  of  data  points  and  mean 
of  the  individual  data  sets,  respectively. 

For  the  calculation  of  the  standard  deviation  using  Equations  (2.4-1),  (-6),  and 
(-8) 

2 Xj  = N X 
N 

2 (x.)2  - (SD)2  (N-1)  + N(x)2 

N 

and  substitution  into  Equation  (2.4-8) 


(SDi2)2  (Ni  + N2  - 1)  = (SD)i2  (Ni  - 1)  + Ni  (x>i2  + (SD)22  (N2-I)  + N2(x)22 

(Nj  Xi  + N2  Xj)® 

N1  + N2  . 


and  the  combined  standard  deviation  of  the  LOX  PBP  135°-1  and  LOX  PBP  135°-2  is 


(2.4-11) 


The  procedure  can  be  repeated  if  an  additional  data  set,  for  example  the  LOX  PBP  RAD 
45°,  is  to  be  combined  with  the  LOX  PBP  RAD  135°- 1 and  LOX  PBP  RAD  135°-2  data  set, 
A more  elegant  derivation  of  the  equations  for  combining  multiple  data  sets, 
augmenting  a data  set  and  deleting  a data  point  from  unpublished  work  of  Tom  Coffin 
is  included  in  the  next  section.  A procedure  for  combining  data  based  on  the  small 
sample  size  theorem  for  the  high  order  moments  (skewness  and  kurtosis)  will  be 
forthcoming  in  future  documentation. 
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2.5  Formulae  for  Estimating  Basic  Engine  Parameter  Statistics 


2.5.1  Background 


Given  a sample  of  size  N from  a population  of  size  P,  the  mean  and  variance  are 
defined  for  N and  P by 


m * (1/N)  ^ ^ X. 

N 

(Sample  mean) 

(2.5-1) 

m = (1/P)  ^ ^ X. 
P 

(Population  mean) 

(2.5-2) 

= (1/N-l)  ^ (x.-in)^ 
N 

(Sample  variance) 

(2.5-3) 

= (1/P)  ^ (x.-m)^ 
P 

(Population  variance) 

(2.5-4) 

Note  that  for  large  N, 

(N/N-D-*!, 

so  formula  3 approaches 

- (1/N)  (x.-m)^ 

(for  large  N) 

(2.5-3A) 

N 


The  asymptotic  convergence  between  formulae  2.5-3  and  -3 A is  more  strongly  emphasized  by 
observing  that  the  standard  deviation  is  the  primary  statistic  of  interest,  so  the 
practical  deviation  between  formulae  2.5-3  and  -3 A is  as 

Since  the  population  of  possible  engine  parameter  outcomes  is  infinite,  our  statis- 

'a 

tical  estimates  are  based  on  formulae  2.5-1  and  2.5-3  or  -3 A. 
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2.5.2  Combined  Multiple  Data  Groups 

Next,  assume  a sample,  S,  consisting  of  three  subgroups,  |x},  {y\ , ^z|.  Technically 

S - XUYUZ  (2.5-5) 

A practical  example  would  be  to  let  X,  Y,  Z represent  a parameter  measured  on  test 
stands  Aj,  Ajj,  Ajjj,  respectively.  The  set  S then  represents  the  parameter  variation 
givii^  equal  weight  to  each  test  from  all  test  stands. 


Assume  also  that  we  have  computed  in  and  a for  each  of  these  groups  and  wish  to 
compute  the  same  statistics  for  the  complete  set  (2.5-5).  It  is  desirable  to  use  the 
group  statistics  already  calculated  to  obtain  these  results,  as  opposed  to  applying 
formulae  2.5-1  and  2.5-3  on  the  complete  data  set.  More  specifically,  given 

m^,  mY)  ^2’  ^Z 


compute  mg,ag,  where  S is  defined  by  formula  2.5-5. 


First,  consider  the  mean  iiig.  By  formula  2.5-1, 


mg  = (1/Ng) 


E V E V E 

«X  «Y  **2 


where  * ^X  ^Y  ^ ^z’ 


(2.5-6) 


(2.5-7) 


But  also  by  formula  2.5-1, 

^X^X 


X.  = N.„m, 


N 


X 


E Yi  = N 

E ^ 


y”"y 


z”^z 


N. 


SO 


mg  = (1/Ng)  [Nj^m^  + NyinY  + N^m 


z] 


(2.5-8) 


(2.5-9) 
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Or,  more  compactly, 


where 


(2.5-10) 


Next  consider  the  variance,  CTg.  By  formula  2.5-3, 


?!  = (1/N, 


1 2]  * 2^  <yr“s>^  * 2^  I 


(2.5-11) 


Expanding  the  first  sum  in  brackets, 


(x.-mg)^  = ^ X?  - 2N5^m2jmg  + N^mg, 


(2.5-12) 


and  noting  that  by  formula  2.5-3, 


the  first  sum  in  formula  2.5-11  may  be  written 


(x.-mg)^  = (Nj^-1)  + Njr(m j._mg)'‘ 


(2.5-14) 


and  an  identical  form  follows  for  the  remaining  two  sums. 


The  variance  of  the  combined  data  set  may  therefore  be  conveniently  computed  by  the 
formula 


4 = ^l/Ng-l)  + 2 Ni(nii-mg)2 


(2.5-15) 
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Again,  if  the  data  groups  are  large  (say,  N.  is  greater  than  30), 


(2.5-16) 


2.5.3  Augmenting  a Data  Set 

Often,  means  and  variances  may  be  computed  for  a set  of  data  and  a new  data  point 
then  obtained.  It  is  desirable  to  compute  these  statistics  for  the  newly  augmented  data 
set  using  those  already  obtained  (as  opposed  to  starting  anew). 


Consider  a group  |x|,  of  size  N^,  for  which  we  have  already  computed  m^,  by 
formuale  2.5-1  and  2.5-3,  respectively,  and  let 

S = |xp  X2,  . . .,  Xj^,  z|  (2.5-17) 

then,  by  formula  2.5-1, 


ms  = (1/Nx+l) 


(2.5-18) 


or  simply 


ins  ""  (1/N^+l)  (N^m^^+Z) 


(2.5-19) 


Finally,  by  formulae  2.5-3  or  2.5-11, 


N, 


and  using  2.5-17, 


4 


(l/Nj^)  |(Nx-l)  ^ (Z-mg)^ 


2.5-19  and  2.5-21  provide  the  desired  formulae. 


(2.5-20) 


(2.5-21) 
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2.5.4  Deleting  a Data  Point 

A third  situation  is  encountered  when  the  mean  and  variance  for  a data  set  have  been 
computed,  after  which  it  is  found  that  an  included  test  result  is  invalid.  It  is  then 
convenient  to  have  formulae  available  to  obtain  the  statistics  for  the  reduced  set,  for 
example,  neglecting  the  invalid  data  point. 

Let  S,  S’  represent  the  complete  and  reduced  data  set,  respectively.  By  simply 
rearranging  2.5-19,  the  mean  for  the  reduced  set  is 

- _ Ngmg  - Z (2.5-22) 

Ng-1 

where  Z represents  the  invalid  data  point  to  be  deleted.  Similarly,  the  desired 
variance  is  obtained  by  manipulating  formula  2.5-21: 

(Ng-l)[?|  - - (Z-5is)4  <2'5-23) 

The  above  formulae  are  easily  programmed  on  a small  calculator  for  ready  reference, 
and  have  been  implemented  in  the  diagnostic  data  base  software  for  data  base  editing- 
and  update. 


<7g,  = (l/Ng-2) 
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2.6  Classical  Distribution  Functions 


Three  classical  continuous  probability  distributions  were  selected  for  comparison 
with  the  SSME  vibration  data.  They  were  the  Normal,  Rayleigh  and  Gamma  functions 
and  are  listed  as  programmed  in  the  Diagnostic  Data  Bank  in  the  next  paragraph. 
Since  the  choice  of  a classical  probability  function  for  a data  set  is  arbitrary, 
additional  investigations  are  planned  to  evaluate  other  classical  distributions  and 
functions  of  higher  order.  Best-Fit  (Chi-Square  Goodness-of-Fit,  Bonferroni-Type 
Inequalities,  etc.)  will  also  be  investigated.  However,  a visual  inspection  of 
Figures  9 to  32  indicates  the  Gamma  function  is  reasonable  for  the  SSME  high 
pressure  turbopump  vibration  data  from  static  firing  tests  and  should  be  applicable 
to  the  flight  data.  Plots  are  shown  for  both  the  composite  and  synchronous  vibra- 
tion levels.  Static  firing  tests  were  utilized  for  the  comparison,  since  a larger 
data  base  was  available  consisting  of  over  1,000  data  samples. 

The  application  of  a classical  distribution  is  desirable  for  data  characterization 
since  this  permits  continuous  statistical  definition  and  manipulation  from  discrete 
flight  measurement  observations.  The  cumulative  distribution  plots  provide  a quick- 
look  assessment  of  flight  results  with  the  historical  distribution  of  measurements 
from  previous  flights.  Density  plots  are  useful  for  an  assessment  of  the  historical 
data  scatter  or  dispersion  around  the  mid-point  or  mean  value.  The  density  func- 
tions are  shown  in  Figures  15  to  20  for  the  composite  values  and  Figures  27  to  32 
for  the  synchronous  values  from  which  the  cumulative  distribution  was  calculated 
from 

N 

F(NAx)  = ^ f(Xi)  (2.6-1) 

i=l 

where  N Ax  is  the  Grms  level  used  in  plotting. 

Normal  Distribution  _ 


f(x) 

= 

r. 

l/syiF  e L J Ax 

(2.6-2) 

M 

= 

Mean  Grms  Level 

S 

Standard  Deviation 

AX 

= 

Step  Size  = 0.5  Grms 
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Rayleigh  Distribution 


fW  - 4x  (2.6-3) 

a* 

AX  = Step  Size  = 0,5  Grms 

The  parameter  a of  the  Rayleigh  density  can  be  the  mean,  standard  deviation,  and 
the  first  moment  or  second  moment  of  the  Rayleigh  function  (Reference  2).  A more 
detailed  study  and  discussion  of  the  Rayleigh  density  and  other  classical  functions 
will  be  included  in  future  documentation. 

Rayleigh  (M) 
a 

Rayleigh  (SD) 
a 

Rayleigh  (MR) 
a 


Rayleigh  (TC) 

a 


Mean 


Standard  Deviation 


M 


y nl2 


M “ Mean 


y/l-nlZ 


S = Standard  Deviation 


Gamma  Distribution 


The  Gamma  distribution  will  be  presented  in  more  detail,  since  it  appears  to  provide 
the  most  reasonable  fit  to  the  vibration  data.  It  should  be  noted  that  the  Gamma 
distribution  contains  two  parameters  related  to  both  the  mean  and  standard  devia- 
tion. The  four  Rayleigh  functions  are  one  parameter  systems  and  utilize  the 
mean  or  standard  deviation,  but  not  both.  The  Gamma  probability  distribution  is 
defined  as 

f (x)  = — ( X x)^*^  e ■ ^ * X > 0 (2.6-4) 

T(r) 
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where  Fir)  is  the  Gamma,  or  generalized  factorial  function. 


The  parameter  X and  r in  terms  of  the  mean  and  variance  from  Reference  (3)  are: 


A.’  (T 


2 = r/X^ 


Solving  each  equation  for  r 


X A.  , r = ^2 


.22  a X X 

X a - AX 


A = x/a^ 

From  the  first  of  Equation  (2.6-5)  and  (2.6-8) 


x2/<7 


Using  Equations  (2.6-4),  (-8)  and  (-10) 


a2rtxV<T*] 


Some  important  relationships  of  the  Gamma  function  are; 


r(n+l)  = 


x"  e dx,  which  leads  to 


r(n)  = (n-1)  r(n-l) 

r(n)  » (n-1)  ! 

r(n+l)  = n ! 


Thus  the  definition  as  a generalized  factorial  function. 
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Therefore  in  terms 

of  the  mean  Grms 

and  standard 

deviation 

the  Gamma  density 

function  is 

m mx 

[(mVs^)  -1] 

(mjc) 

" c2 

f(x) 

, , ■ ■ 

s^[mVs^-l]!  s* 

♦ 

s 

e 

AX  (2.6-13) 

m = 

Mean  Grms  Level 

s = 

Standard  Deviation 

A X 

Step  Size  = 0.5 

Although  the  Gamma  distribution  appears  to  provide  a reasonable  fit  to  the  data, 
some  additional  improvements  should  be  investigated.  These  include  a change  in  the 
class  interval  width,  unequal  class  intervals  and  maybe  some  type  of  truncated 
Gamma  or  other  classical  function  to  account  for  the  low  synchronous  values  ap- 
proaching the  noise  floor  of  the  instrumentation.  It  may  be  noted  F(x)  and  f(x) 
are  not  the  continuous  form  of  the  cumulative  probability  function  and  probability 
density  function,  but  represent  percentiles  within  discrete  class  intervals.  This 
formulation  permits  direct  comparison  between  classical  function  approximations  and 
discrete  frequency  distributions  representing  empirical  measurement  characteristics. 
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Figure  9.  Gumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Composite)  Normal  Overlay 
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Figure  14.  Cumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Composite)  Gamma  Overlay 
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Figure  15.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Composite)  Normal  Overlay 
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Figure  16.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Composite)  Rayleigh  (M)  Overlay 
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Figure  17.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Composite)  Rayle^h  (SD)  Overlay 
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Figure  18.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Composite)  Rayleigh  (MR)  Overlay 
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Figure  19.  Probability  Density  (LOX  PBP  RAD,  Static  Firii^,  Composite)  Rayle^h  (TC)  Overlay 
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Figure  20.  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Composite)  Gamma  Overlay 
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Figure  21.  Cumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Normal  Overlay 
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Figure  22.  Cumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Rayleigh  (M)  Overlay 
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Figure  24.  Cumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Rayleigh  (MR)  Overlay 
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Figure  26.  Cumulative  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Gamma  Overlay 
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Figure  27.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Normal  Overlay 
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Figure  28.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Rayle^h  (M)  Overlay 
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Figure  29.  Probability  Density  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Rayleigh  (SD)  Overlay 
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Figure  30.  ProbabiUty  Density  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Rayleigh  (MR)  Overlay 
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Figure  32.  Probability  Distribution  (LOX  PBP  RAD,  Static  Firing,  Synchronous)  Gamma  Overlay 


2.7  Sort  Routine  for  Probability  Density  Function 


A variation  of  the  following  quick  sort  routine  is  used  in  the  MSFC  Diagnostic  Data 
Base  Program  to  calculate  the  density  functions.  The  procedure  follows  Reference 
4.  This  method  requires  only  one  data  pass  and  therefore  is  basically  the  fastest 
most  simple  routine  available. 

Consider  the  N data  values  for  the  composite,  synchronous  etc.  for  each  measurement. 
The  probability  density  function  can  be  estimated  by 


f(x) 


Nx 
N W 


where  W is  a narrow  interval  centered  at  x and  is  the  number  of  data  values 

A 

which  fall  within  the  range  x + W/2.  A general  procedure  to  estimate  f(x)  can  be 
obtained  digitally  by  dividing  the  full  range  of  x into  an  appropriate  number  of 
equal  width  class  intervals.  Tabulate  the  number  of  data  values  in  each  class 
interval,  and  divide  by  the  product  of  the  class  interval  width  W and  the  sample 

A * 

size  N.  The  estimate  f(x)  is  not  unique  since  it  is  dependent  upon  the  number  of 
class  intervals  and  their  width  selected  for  the  analysis.  Equal  class  intervals  of 
0.5  Grms  were  used  in  the  MSEC  Diagnostic  Data  Bank  Program.  The  number  of  class 
intervals  is  equal  to  twice  the  maximum  value  of  the  plot  with  the  data  plotted  at  x 
rather  than  x + W/2.  This  is  not  in  strict  accordance  with  the  above  definition, 
but  for  data  interpretation  the  plots  show  the  percentage  of  data  values  less  than 
or  equal  to  x,  where  x is  the  Grms  value. 

A general  procedure  for  evaluation  is  as  follows.  Let  K define  the  number  of  class 
intervals  selected  to  cover  the  entire  range  of  the  data  values  from  a^,  to  bj^. 
Then  the  width  of  each  interval  is  given  by 


and  the  end  point  of  the  ith  interval  is 

nj  = a„  + i W i = 0,  1,  2,  . . . K 
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where  = a,,  and  = b^.  The  bin  numbers  will  now  be  defined  as  a sequence  of 
K+2  numbers  [NJ,  i « 0,  1,  2,  . . . K+1,  by  the  conditions 

Nj,  = [number  of  x such  that  x < n J 

Nj  =>  [number  of  x such  that  n^  < x < n^] 

Nj  = [number  of  x such  that  ni_i  < x < Uj] 


Njj  = [number  of  x such  that  n^.^  < x < nJ 

= [number  of  x such  that  x > nJ 

This  procedure  will  sort  out  the  N data  values  of  x so  that  the  bin  number  sequence 
[Nj]  satisfies 

K+1 

N - 2 N, 

i=G 

A quick  one  pass  method  of  sorting  on  a digital  computer  is  to  evaluate  each  value 
Xjjj  n=l,  2, . , . , N in  turn  as  follows. 

a.  if  Xj,  < a„,  add  the  integer  one  to  Np. 


b. 


If  a^  < < bj^,  compute  I = 


then  select  i as  the  largest  integer  less  than  or  equal  to  I,  and  add 
the  integer  one  to  Nj. 

c.  If  x^  > b„j,  add  the  integer  one  to 

The  four  output  forms  are  histogram,  percentage  of  data  in  each  class  interval, 
probability  density  estimate  and  probability  distribution  estimate. 
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The  histogram  is  simply  Nj  without  changes 


Hj  = Nj  Histogram 

where  Hj  is  the  number  of  data  values  in  each  class  interval. 

The  second  output  is  the  sample  percentage  of  data  in  each  class  interval  defined 
for  i = 0,  1,  2,  . . . K+1  by 

^ Nj 

Fj  % = — — X 100  Percentage  of  Data 

N 


The  third  output  is  the  probability  density  estimate  defined  at  the  midpoints  of  the 
K class  intervals  in  [a„,  b^]  by 


i=l,  2,  3 . 


. . K 


Probability  Density 


and  the  probability  distribution  estimate  defined  at  the  class  interval  end  points 
where  i = 0,  1,  2,  . . . , K+1. 


j=0  j=0 


Ratio  of  Synchronous  to  Composite 


The  analysis  of  the  ratio  of  synchronous  to  composite  vibration  levels  measured  in 
flight  is  documented  for  future  reference.  This  data  will  be  utilized  in  studies 
to  relate  flight  to  static  firing  data,  the  effect  of  . vibration  level  amplitude  on 
the  ratio  of  synchronous  to  composite,  evaluation  of  the  FASCOS  system  and  justifi- 
cation of  vibration  redlines. 
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The  ratio  of  synchronous  to  composite  vibration  levels  of  Figures  33,  34,  35,  and 
36  for  100%  and  104%  power  levels  were  calculated  using  the  data  stored  in  the  MSFC 
Diagnostic  Data  Base.  Both  the  composite  (50-1000  Hz)  and  synchronous  vibration 
levels  in  the  data  base  are  derived  from  a sliding  11  point  average  (0.4  second 
PSDs)  which  is  equal  to  4.4  seconds  of  data  during  each  power  level.  The  highest* 
or  maximum  value  of  this  average  is  then  stored  in  the  data  base**  to  represent  the 
characteristic  Grms  vibration  level  at  that  power  level.  Since  the  Diagnostic  Data 
Base  Program  does  not  record  or  track  the  time  of  occurrence  during  the  flight  the 
ratio  may  or  may  not  represent  the  same  time  period  of  the  flight.  However,  past 
history  has  shown  the  composite  and  synchronous  track  in  a reasonable  manner  during 
most  static  firing  tests  and  flights.  This  pseudo  ratio  shown  in  the  figures  should 
therefore  be  fairly  representative  of  a ratio  calculated  from  data  recorded  at  the 
same  time  period.  Additional  investigations  are  planned  to  evaluate  the  skewed 
distribution  on  the  high  pressure  fuel  turbopump  and  the  appearance  of  a bi-modal 
type  of  distribution  on  the  high  pressure  oxidizer  turbopump.  One  possibility  is 
the  ratio  of  synchronous  to  composite  is  a function  of  the  vibration  level.  Also 
the  investigation  will  include  a comparison  with  the  measured  static  ground  test 
data. 


*The  vibration  level  stored  in  the  Diagnostic  Data  Base  and  reported  at  data  review 
meeting  could  be  called  the  "highness"  vibration  level.  Highness  as  defined  in  the 
Webster  Dictionary  "the  highest  of  the  nobility"  with  noble  defined  as  "having 
excellent  qualities;  superior."-  The  MSFC  Technical  Monitor  for  this  report  has 
suggested  contiguous  average  maximum  vibration  level  as  a descriptive  term,  while 
the  Wyle  Technical  Reviewer  has  suggested  noble  vibration  level.  Only  historical 
usage  will  finally  define  the  terminology, 

**The  computer  also  skips  random  clipped  signals  and  the  output  data  is  analyzed  to 
purge  the  data  base  of  invalid  and  extreme  noisy  data. 
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Figure  34.  Probability  Density  of  Synehronous/Composite  Ratio 
for  the  HPFTP  at  104%  Power  Level  During  Flight 
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F^ure  36.  Probability  Density  of  Synchronous/Gomposite  Ratio 
for  the  HPOTP  at  104%  Power  Level  During  Flight 


2.9  TurboBumn  Rotational  Speed  History 


The  rotational  speed  of  the  high  pressure  oxidizer  turbopump  (HPOTP)  and  high 
pressure  fuel  turbopump  (HPFTP)  are  shown  in  Figures  37  and  3$.  The  speed  is 
calculated  from  the  synchronous  frequency  of  the  measured  vibration  data  and 
represents  the  maximum  speed  of  the  pump  during  the  flight  at  100%  and  104%  power 
levels.  The  notation  on  the  plot  lists  the  flight  number  STS  xx,  engine  position 
Ex  and  turbopump  serial  number. 
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Figure  37.  HPOTP  Maximum  Rotational  Speed  During  Fl^ht 
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Figure  38.  HPFTP  Maximum  Rotational  Speed  During  Flight 
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4.0  RESULTS 


A summary  of  the  mean  vibration  levels,  the  standard  deviation,  and  maximum  level 
measured  on  the  SSME  high  pressure  oxidizer  and  fuel  turbopump  during  flight  is 
shown  in  Figures  39  through  50.*  The  results  are  for  the  SSME  operating  at  104%, 
100%  and  65%  power  levels.  Statistics  are  for  each  measurement  location  and  the 
average  of  all  locations  which  represents  a spatial  average  around  the  turbopump. 
Both  the  composite  (50-1000  Hz)  and  synchronous  vibration  levels  were  considered  in 
this  investigation. 

Cumulative  probability  distribution  plots  are  Shown  in  Figures  51  through  62  with 
the  classical  Gamma  distribution  plotted  as  an  overlay  dashed  line.  These  plots  as 
discussed  in  Section  2.6  can  provide  a quick-look  assessment  of  flight  results  as 
compared  to  this  historical  distribution  of  the  previous  flight  data.  The  probabil- 
ity density  for  each  power  level  for  the  oxidizer  and  fuel  turbopumps  is  shown  in 
Figures  63  through  74.  These  plots  provide  a quick  assessment  of  the  historical 
flight  data  scatter  or  dispersion  from  the  mid-point  or  mean  value  whichever  may  be 
of  interest.  The  tabulated  data  on  each  plot  is  the  number  of  data  points  in  each 
0.5  Grms  interval.  Also  the  number  of  tests  is  equal  to  the  number  of  test  samples 
rather  than  the  number  of  flights. 

For  future  reference,  flight  numbers  keyed  into  the  computer  are  listed  for  the 
flights  used  in  this  analysis.  The  data  sheets  for  each  measurement  are  also 
included  for  reference. 


The  data  of  STS  Launch  27  (51-1)  and  28  (51-J)  are  included  in  Appendix  A. 
Both  flights  were  nominal  and  will  not  have  a significant  effect  on  the 
statistical  data  of  Figures  39  through  74. 
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FLIGHT  summary  SHEETS 
HPOTP  AND  HPFTP 
104%,  100%  AND  65%  POWER  LEVEL 
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F^ure  39.  Summary  of  Composite  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  104%  Power  Level  During  Flight 


Synchronous  9I04X  Power  Level 

FUEL  PUriP  RAD  0 FUEL  PUMP  RAD  90  FUEL  PUMP  RAD  174 
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F^re  40.  Summary  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fu^  Turbopump  at  104%  Power  Level  During  Flight 


Composite  6104X  Power  Level 

LOX  PBP  RAD  45  LOX  P8P  RAD  135-1  LOX  PBP  RAD  135-2 
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Figure  41.  Summary  of  Composite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  104%  Power  Level  During  Flight 


Synchronous  9104X  Power  Level 

LOX  PBP  RAD  45  LOX  PBP  RAD  135-t  LOX  PBP  RAD  135-2 
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High  Pressure  Oxidizer  iSirbopump  at  104%  Power  Level  During  Flight 


Compos! te  8i00X  Power  Level 

FUEL  PUMP  RAD  0 FUEL  PUMP  RAD  90  FUEL  PUMP  RAD  174 
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Synchronous  dl00X  Power  Level 

FUEL  PUMP  RAD  0 FUEL  PUMP  RAD  90  FUEL  PUMP  RAO  174 
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Compos? 76  6^00%  Power  Level 

LOX  P8P  RAD  45  LOX  PBP  RAD  135-1  LOX  PBP  RAD  135- 
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Figure  45.  Summary  of  Ck>mposite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  100%  Power  Level  During  Flight 


Synchronous  8100%  Power  Level 
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Compoeite  9 65X  Power  Level 

FUEL  PUMP  RAD  0 FUEL  PUMP  RAD  90  FUEL  PUMP  RAD  174 
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F^ure  47.  Summary  of  Composite  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  65%  Power  Level  During  Flight 


Synchronous  6 65X  Power  Level 

FUEL  PUMP  RAD  0 FUEL  PUMP  RAD  90  FUEL  PUMP  RAD  174 
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Figure  48,  Summary  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  65%  Power  Level  Durii^  Flight 
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Figure  49.  Summary  of  Ck>mposite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  65%  Power  Level  During  Flight 


Synchronous  8 65X  Power  Level 
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FLIGHT  CUMULATIVE  PROBABILITY  DISTRIBUTIONS 
HPOTP  AND  HPFTP 
104%,  100%  AND  65%  POWER  LEVEL 
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Figure  52.  Cumulative  Distribution  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  104%  Power  Level  During  Flight 
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Figure  53.  Cumulative  Distribution  of  Composite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  104%  Power  Level  During  Flight 


Figure  54.  Cumulative  Distribution  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  104%  Power  Level  During  Flight 
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Figure  57.  Cumulative  Distribution  of  Composite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  100%  Power  Level  During  Flight 
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Figure  58.  Cumulative  Distribution  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  100%  Power  Level  During  Flight 
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Figure  59.  Cumulative  Distribution  of  Composite  Vibration  Leve^  on  the 
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Figure  60.  Cumulative  Distribution  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  65%  Power  Level  During  Flight 


Figure  61.  Cumulative  Distribution  of  Composite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  65%  Power  Level  During  Flight 


Figure  62.  Cumulative  Distribution  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  65%  Power  Level  During  Flight 
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Figure  63.  Probability  Density  of  Composite  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  104%  Power  Level  During  Flight 
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F^ure  64.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  104%  Power  Level  During  Flight 
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Figure  66.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  104%  Power  Level  Durir^  Flight 
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Figure  67.  Probability  Density  of  Composite  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  100%  Power  Level  During  Fl^ht 
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Figure  68.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  100%  Power  Level  During  Flight 


Figure  70.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  100%  Power  Level  During  Flight 


Figure  71.  Probability  Density  of  Composite  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  65%  Power  Level  Durii^  Flight 
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Figure  72.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Oxidizer  Turbopump  at  65%  Power  Level  Durir^  Flight 


Figure  73.  Probability  Density  of  Composite  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  65%  Power  Level  During  Flight 
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Figure  74.  Probability  Density  of  Synchronous  Vibration  Levels  on  the 
High  Pressure  Fuel  Turbopump  at  65%  Power  Level  During  Plight 
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1.0  INTRODUCTION  AND  CONCLUSIONS 


The  purpose  of  this  report  is  threefold:  1)  to  document  an  investigation  of  higher 

than  expected  (number  of  occurrences  and  amplitude)  synchronous  vibration  levels  on 
the  high  pressure  oxidizer  turbopump  (HPOTP),  2)  to  illustrate  the  value  of  basic 
fundamental  statistical  methods  of  data  analysis,  incorporated  in  the  NdFC  Diag- 
nostic Data  Base  Program,  and  3)  to  document  a data  list  of  tests  where  the  HPOTP 
vibration  data  are  not  contaminated  by  high  running  main  impellers,  whirl,  instru- 
mentation malf  unctions,  etc.  Both  the  HPOTP  Preburner  Pump  (PBP)  and  Turbine  (TURB) 
measurements  were  considered  in  the  analysis,  operating  at  100%,  104%  and  109%  power 
levels. 

The  terminology  "High  Running  Main  Impellers"  was  applied,  after  the  study  was 
completed,  to  an  investigation  in  the  last  quarter  of  1982  and  early  1983  of  higher 
than  expected  (number  and  amplitude)  synchronous  vibration  levels  on  the  HPOTP. 
When  the  original  study  was  performed  the  vibration  data  was  hand  plotted  in  a 
format  similar  to  the  plotting  routines  now  available  using  the  MSFC  Diagnostic  Data 
Base  Program.  Data  in  this  reports  documents  a re-evaluation  and  update  of  the 
original  effort  with  maximum  utilization  of  the  MSFC  computer  program  for  data 
sorting,  analysis  and  plotting.  Analysis  included  calculation  of  the  mean,  standard 
deviation,  spatial  average,  cumulative  distribution  and  probability  density  with  the 
classical  Gamma  distribution  plotted  as  an  overlay.  Elegant  and/or  sophisticated 
statistical  techniques  were  not  used  or  required  in  the  analysis,  since  visual 
methods  were  adequate  to  show  a distinct  difference  between  data  groups. 

This  report  does  not  change  the  conclusions  and  results  of  the  original  study  which 
indicated  a strong  correlation  between  six  main  impellers  and  the  high  synchronous 
vibration  levels.  The  analysis  does,  however,  illustrate  the  application  of 

probability  density  estimates  for  detecting  embedded  measurement  groups  representing 
differing  statistical  (and  likely  physical)  behavior.  In  addition,  the  statistical 
characterization  of  the  vibration  data  provided  herein  represents  nominal  HPOTP 
operation  for  future  reference. 
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2.0  TECHNICAL  DISCUSSION 


2.1  Diagnostic  Data  Base 

Table  I is  a listing  of  the  relevant  vibration  measurements  stored  in  the  MSFC 
Diagnostic  Data  Base  Program.  The  listing  includes  the  test  number,  date  of  test, 
engine  serial  number,  HPFTP  serial  number,  HPOTP  serial  number  and  the  power  level 
of  which  the  vibration  data  was  acquired.  For  this  analysis  the  number  of  tests  at 
each  power  level  was: 

407  tests  @ 100% 

155  tests  @ 104% 

234  tests  @ 109%. 

These  tests  resulted  in  the  following  synchronous  vibration  data  sample  for  analysis 
at  the  HPOTP  preburner  pump  and  turbine  measurement  locations.  Each  measurement 
was  treated  as  an  individual  data  point  and  therefore  the  results  represent  the 
spatial  average  over  the  transducer  locations. 


Number  of 

Power 

Data  Samoles 

Level 

Measurement  Location 

1299 

100% 

LOX  PBP  RAD  45,  135-1,  135-2, 

542 

104% 

180,  225 

712 

109% 

971 

100% 

348 

104% 

LOX  TURB  RAD  45,  90,  135 

546 

109% 

Serial  numbers  of  HPOTP’s  with  the  suspect  main  impellers  discussed  in  Section  2.4 
are  highlighted  and  the  tests  marked  with  an  "H".  Also  the  tests  that  experienced 
whirl  are  marked  with  an  "SS"  and  the  power  level  noted  in  the  remarks  column. 
These  tests  were  also  deleted  before  calculation  of  the  statistics  for  normal 
turbopump  operation. 

2.2  Vibration  Measurement  Location 

The  plane  in  which  the  preburner  pump  and  turbine  accelerometers  are  located  is 
shown  in  Figure  1.  All  accelerometers  were  oriented  to  measure  the  vibration  in 
the  radial  direction.  In  Figure  2 the  major  components  of  the  HPOTP  pump  section 
are  identified  including  the  main  impeller. 
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TABLE  I.  SSHE  HIGH  PRESSURE  QXjDIZER  AND  FUEL  TURBOPUMP 
VIBRATION  DATA  IN  MSFC  DIAGNOSTIC  DATA  BASE 
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Figure  1.  Accelerometer  Location  Piane  HPOTP 


2.3  Total  Data  Base  HPOTP  Synchronous  Vibration 

All  tests  listed  in  Table  I were  used  to  generate  the  plots  shown  in  Figures  3 to 
14.  Separate  plots  are  shown  for  100%,  104%,  and  109%  power  levels.  These  plots 
are  the  cumulative  and  density  histograms  of  all  the  relevant  synchronous  Grms 
vibration  data  presently  in  the  MSFC  Diagnostic  Data  Base.  Immediately  obvious, 
especially  at  the  104%  power  level,  is  the  bi-modal  type  of  distribution  which 
indicates  the  possibility  of  two  different  data  groups.  Also  of  note  is  the  fact 
that  a bi-modal  type  of  distribution  does  not  appear  on  the  turbine  end  vibration 
measurements. 

In  the  last  quarter  of  1982  a major  effort  was  instigated  to  evaluate  the  second 
group  of  tests  i.e.  tests  with  synchronous  levels  greater  than  approximately  4 
Grms.  Numerous  MSFC  and  contractor  personnel  were  - involved  in  the  effort  and 
investigated  different  facets  of  the  problem  such  as  imbalance,  tear  down  inspection 
records,  common  components,  etc.  for  each  vibration  data  group.  The  conclusion  from 
the  original  study  was  that  the  increased  synchronous  vibration  levels  were  associ- 
ated with  the  HPOP  Main  Impeller, 

2.4  Listing  of  Suspect  Main  Impeller 

The  results  of  the  MSFC/NASA  investigation  specifically  identified  six  main  impel- 
lers common  to  the  second  data  group  of  tests,  i.e.  tests  with  synchronous  vibration 
levels  greater  than  4 Grms.  A total  of  33  main  impellers  were  investigated.  The 
serial  numbers,  tests,  power  levels  and  the  serial  number  of  the  HPOTP  where 
the  main  impeller  was  installed  are  shown  in  Tables  II  through  VII.  These  tables 
are  a modified  version  of  output  of  a computer  routine  available  within  NASA/MSFC 
for  tracking  components  of  the  SSME. 
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Figure  3.  Cumulative  Distribution  HPOTP-PBP  Total  Data  Base  104% 
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Figure  4.  Probability  Density  HPOTP-PBP  Total  Data  Base  104% 

30 


______  Synchronous  1 04’:  PMR  iUL  12-M0U-85 

TEST  »'S  ftl2Sl-495,A2201,A2202-384,A3l25-26l, 

Hunber  of  tests  =348 


— — Synchronous  104H  PHR  LML  12-H0U-85 

TEST  #'S  A1291-495,A2201iA22O2-334,A3l25-26l, 

Nunlosr  of  tests  >348 


Figure  6.  Probability  Density  HPOTP-TURB  Total  Data  Base  104% 
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Figure  7.  Cumulative  Distribution  HPOTP-PBP  Total  Data  Base  100% 
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Figure  8.  Probability  Density  HPOTP-PBP  Total  Data  Base  100% 
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Figure  9.  Cumulative  Distribution  HPOTP-TURB  Total  Data  Base  100% 
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Figure  12.  Probability  Density  HPOTP-PBP  Total  Data  Base  109% 
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Figure  13.  Cumulative  Distribution  HPOTP-TURB  Total  Data  Base  109% 
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TABLE  n,  HPOP  MAIN  IMPELLER  RI3TORY  - S/N  2427800 
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Total  S/N  2427800  5992.68  581.14  2212.32  2288.70  19.80 


TABLE  ra.  HPOP  MAIN  IMPELLER  HISTORY  - S/N  3134124 
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TABLE  IV.  HPOP  MAIN  IMPELLER  HISTORY  - S/N  3134446 
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TABliE  VI.  HPQP  MAIN  IMPELLER  HISTORY  - S/N  UUm 
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TABLE  Vn.  HPOP  MAIN  IMPELLER  HISTORY  - S/N  7363066 
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2.5  Synchronous  Vibration  Data  of  Six  Main  Impellers 

The  next  step  in  the  investigation  was  to  repeat  the  analysis  with  the  six  main 
impeller  vibration  test  data  separated  from  the  data  sample.  The  results  are  shown 
in  Figures  15  to  22  at  the  104%  power  level.  These  plots  present  a very  strong 
implication  of  two  distinct  data  groups  and  indicate  the  listed  six  main  impellers 
as  very  suspect.  This  conclusion  is  based  upon  the  obvious  visual  difference 
between  the  density  histograms  of  Figures  4,  15  and  16.  No  elegant  statistical 
mathematical  methods  are  required.  Also  it  should  be  noted  that  the  cumulative 
distribution  after  removal  of  the  suspect  main  impeller  data  provides  a much  better 
fit  to  the  classical  Gamma  distribution  discussed  in  Reference  1.^  Figures  23  to  38 
are  the  plots  at  100%  and  109%  power  level  and  also  indicate  a difference  between 
groups,  while  not  as  dramatic  as  the  104%  power  level.  Very  little  change  appears 
in  the  LOX  TURB  RAD  vibration  cumulative  distribution  and  probability  density  plots 
which  is  further  evidence  to  suspect  components  on  the  preburner  end  of  the  turbo- 
pump. Since  the  purpose  of  this  document  is  to  record  the  results  of  the  study, 
rather  than  a detailed  discussion  of  statistical  techniques  and  physical  effects,  a 
detailed  analysis  of  each  plot  will  not  be  addressed.  The  difference  in  the 
calculated  mean  Grms  and  standard  deviation  for  the  three  data  groups  are  listed  in 
Table  VIII. 


^ Reference  I.  Swanson,  W.  "Statistical  Analysis  of  the  Vibration  Data 
for  the  SSME  High  Pressure  Turbopump  During  Flight,"  Wyle  TM  64058-01  October  1985. 
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TABLE  VIIL  LOX  PBP  RAD  SYNCHRONOUS  VIBRATION  LEVELS 


Total  Data  Base 

100%  Power  Level 
Six  Main  ImDellers 

Normal  Ooeration 

Mean 

2.35  Grms 

4.61  Grms 

1.89  Grms 

Std  Dev 

1.69 

2.16 

1.12 

Data  Sample 

1299 

219 

1080 

104%  Power  Level 

Mean 

3.06  Grms 

6.57  Grms 

2.09  Grms 

Std  Dev 

2.09 

1.86 

1.05 

Data  Sample 

542 

111 

412 

109%  Power  Level 

Mean 

2.72  Grms 

5.42  Grms 

2.15  Grms 

Std  Dev 

1.72 

2.14 

1.03 

Data  Sample 

712 

114 

576 

Spatial  Average  LOX  TURB  RAD  Synchronous 

100%  Power  Level 

Total  Data  Base 

Six  Main  Impellers 

Normal  Ooeration 

Mean 

1.62  Grms 

2.07  Grms 

1.55  Grms 

Std  Dev 

0.84 

1.08 

0.76 

Data  Sample 

971 

138 

833 

104%  Power  Level 

Mean 

1.64  Grms 

1.80  Grms 

1.59  Grms 

Std  Dev 

0.75 

0.76 

0.73 

Data  Sample 

348 

61 

276 

109%  Power  Level 

Mean 

1.55  Grms 

1.70  Grms 

1.52  Grms 

Std  Dev 

0.74 

0.82 

0.73 

Data  Sample 

546 

82 

449 
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Normal  HPOTP  Operation  104% 


TESTS  USED  IH  THIS  ANALYSIS; 

TEST  #'S  A1291-323,A1390-430,A1433-440.A1446-495.A2201-307,A23U  -319, 

A2323-326, A2329-343, A2345-354, A2357-374, A2383-384, A3125- 181 , A321 1-215, A3 
221-261, 

ENTER... 1)  LIST  STANDARD  TABLE 

2)  PLOT  PROBABILITY  DISTRIBUTION 
3>  PLOT  PROBABILITY  DENSITY 
4>  SELECT  A NEW  SET  OF  TESTS 

5)  SELECT  NEH  DATA  TYPE  AND  PMRLML 

6)  RETURN  TO  MAIN 


• Six  Main  Impellers  104% 

TESTS  USED  IN  THIS  ANALYSIS: 

TEST  #/S  A 1389,  A 1432,  A 1442-443,  A2310,  A2375-376,  A3 1 83-209,  A32 16-2 19 
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Synchronous  104r<  PMR  LML  14-H0M-85 


Hunber  oi  tests  =412 


' Figure  15.  Probability  Density  for  HPOTP-PEP  Normal  Operation  it)4%;~ 


— __ — __  Synchronous  104/:  PWR  LUL  14-110U-85 
TEST  * ' S A I 389 , A 1432 , A 1 442-443 , A23 1 0 , A2375-376 , A3 1 83-209 , A32 16-2 19 , 
Hunber  of  tests  =111 


25 


G rns 

Figure  16.  Probability  Density  for  HPOTP-PBP  Six  Main  Impellers  104% 
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— — - Synchronous  104':  PHR  L0L  14-M0M-85 


tluHber  of  tests  =276 
30 


G r«s 

Figure  17.  Probability  Density  HPOTP-TURB  Normal  Operation  104% 


Synchronous  104*:  PHR  LML  14-M0M-85 
TEST  i-^S  01388,01432, 01442-443, A2310,A2375-3?6,A3i83-209, 03216-219, 
Hunher  of  tests  =61 


Figure  18.  Probability  Density  HPOTP-TURB  Six  Main  Impellers  104% 
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Janchronous  104*<  PHR  LUL  14-H0U-85 
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Figure  19.  Cumulative  Distribution  HPOTP-PBP  Normal  Operation' 104% 


Sanchronous  104J;  PMR  LUL  14-H0M-83 

TEST  H'S  A1389,A1432,A1442-443,A2310,A2375-376,A3183-209,A3216-219, 
Hunber  of  tests  =111 
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Figure  20.  Cumulative  Distribution  HPOTP-PBP  Six  Main  Impellers  104% 
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Synchronous  104';  PMR  LML  14-H0U-85 


HuHber  of  tests  -276 


G r»s 


Figure  21.  Cumulative  Distribution  HPOTP-TURB  Normal  Operation  104% 


Synchronous  104?;  PMR  LML  14-H0M-85 

TEST  #'S  fil389,ftl432,Al442-443,H23i0,H23?5-376, 03183-209, 03216-219, 
lluHber  of  tests  -61 


Figure  22.  Cumulative  Distribution  HPOTP-TURB  Six  Main  Impellers  104% 
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Normal  HPOTP  Operation  100% 
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TESTS  USED  IH  THIS  ANftLYSIS: 

TEST  #-'S  A1 194-280. A1284-341 , ftl358-363, A1 3?3, A1384, A1390-391 , ftl394-395, A 
1406-430, A1433-440, A1445-495, A2145-264, A2271-2F4 ,A2289-303j A2307,A231 1-3 
19, A2323-374, A2377-384, A3674-162, A3177-132, A3184- 188, A3210-215, A3220-262 


ENTER... 1)  LIST  STANDARD  TABLE 

2>  PLOT  PROBABILITY  DISTRIBUTION 

3)  PLOT  PROBABILITY  DENSITY 

4)  SELECT  A NEW  SET  OF  TESTS 

5>  SELECT  NEH  DATA  TYPE  AND  PHRLUL 
6)  RETURN  TO  HAIN 


Six  Main  Impellers  100% 


TESTS  USED  IN  THIS  ANALYSIS; 

TEST  » " S A 1 282 , A 1 347 , A 1366 , A 1 375-383 , A 1 389, A 1393, A 1 397-405 , A 1432 , A 1 442-4 
44, A2265-268, A2285, A2306, A2309-310, A2322, A2375-376, A3170-173, A3183, A3190 
-209,A3216-219, 


ENTER... 1> 
2) 
3) 
4> 
5> 
6) 


LIST  STANDARD  TABLE 
PLOT  PROBABILITY  DISTRIBUTION 
PLOT  PROBABILITY  DENSITY 
SELECT  A NEH  SET  OF  TESTS 
SELECT  NEH  DATA  TYPE  AND  PHRLUL 
RETURN  TO  MAIN 
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Figure  23.  Probability  Density  HPOTP-PBP  Normal  Operation  100% 
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Figure  24.  Probability  Density  HPOTP-PBP  Six  Main  Impellers  100% 
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Figure  25.  Probability  Density  HPOTP-TURB  Normal  Operation  100% 
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Figure  26.  Probability  Density  HPOTP-TURB  Six  Main  Impellers  100% 
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Figure  27.  Cumulative  Distribution  HPOTP-PBP  Normal  Operation  100% 
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Figure  28.  Cumulative  Distribution  HPOTP-PBP  Six  Main  Impellers  100% 
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Normal  HPOTP  Operation  109% 


TESTS  USED  IN  THIS  ANALYSIS: 

TEST  # S A1322-343,A1351-355,A1361-362,A1364-372.A1384-395.A1433-436,A14 
45-493, A21S3-261 ,A226r,A22?2-276, A2286-372. A2377-381, A3125-152, A3154-168 
,A3174-181,A3184-188,h3211-262. 


ENTER... 1)  LIST  STANDARD  TABLE 

2>  PLOT  PROBABILITY  DISTRIBUTION 
3)  PLOT  PROBABILITY  DENSITY 
4>  SELECT  A NEW  SET  OF  TESTS 
5>  SELECT  HEW  DATA  TYPE  AND  PWRLUL 
6)  RETURN  TO  MAIN 


Six  Main  Impellers  109% 


TESTS  USED  IN  THIS  ANALYSIS: 


TEST  # S A1 344-349 , A 1358, A 1 374-38 1 , A 1 398-399 , A 1 49 1 -405 , A 1 443, h2268, A2277 
-284, A2375-376, A3171 , A3 1 72 , A31 83, , A3198- 191 , 


ENTER... I)  LIST  STANDARD  TABLE 

2)  PLOT  PROBABILITY  DISTRIBUTION 
3>  PLOT  PROBABILITY  DENSITY 
4)  SELECT  A HEW  SET  OF  TESTS 
5>  SELECT  NEW  DATA  TYPE  AND  PWRLUL 
6>  RETURN  TO  HAIN 
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Figure  31.  Probability  Density  HPOTP-PBP  Normal  Operation  109%  ^ 
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Figure  32.  Probability  Density  HPOTP-PBP  Six  Main  Impellers  109% 
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Figure  33.  Probability  Density  HPOTP-TURB  Normal  Operation  109% 
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Figure  34.  Probability  Density  HPOTP-PBP  Six  Main  Impellers  109% 


a.  Ul  Ul  X H-  <S  Ul 


— ■ — SyrichroTiOus  IflS’l  PWR  LML  14~H0M-35 
Hunber  o#  tests 


G rns 

Figure  35,  Cumulative  Distribution  HPOTP-PBP  Normal  Operation  109% 
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Figure  36.  Cumulative  Distribution  HPOTP-PBP  Six  Main  Impellers  109% 
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Figure  37.  Cumulative  Distribution  HPOTP-TURB  Normal  Operation  109% 
— ...  SuTichronous  109';  PHR  LML  I5-M0M-85  ' " 
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Figure  38.  Cumulative  Distribution  HPOTP-TURB  Six  Main  Impellers  109% 


58 


2.6  Random  Selection  of  Tests 

A random  sample  of  29  tests  was  selected,  which  is  the  number  of  tests  in  the  six 
main  impeller  data  sample,  to  evaluate  the  effect  of  sample  size.  The  same  analysis 
was  used  for  this  new  data  group  and  the  results  are  shown  in  Figures  39  and  40. 
While  it  is  recognized  this  method  of  analysis  is  not  rigorous,  it  seems  reasonable 
to  assume  a random  selection  should  have  followed  a relative  smooth  continuous  type 
of  distribution  unless  the  data  is  from  two  different  data  groups.  The  distribution 
is  again  bi-modal  and  presents  additional  evidence  that  the  total  data  sample 
consists  of  two  different  data  groups. 
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Randoni  Selection  of  29  Tests 


TESTS  USED  IN  THIS  AHALYSIS: 
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3218, A3221-222, A3227, A3257, 

ENTER...!)  LIST  STANDARD  TABLE 

2)  PLOT  PROBABILITY  DISTRIBUTION 

3)  PLOT  PROBABILITY  DENSITY 

4)  SELECT  A HEM  SET  OF  TESTS 

5)  SELECT  NEW  DATA  TYPE  AND  PHRLUL 

6)  RETURN  TO  MAIN 
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Figure  39.  Cumulative  Distribution  HPOTP-PBP  Random  Selection  of  29  Tests  104% 
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Figure  40.  Probability  Density  HPOTP-PBP  Random  Selection  of  29  Tests  104% 


2.7  Vibration  Test  History  of  the  Six  Main  Impellers 

To  complete  this  analysis,  the  synchronous  spatial  average  vibration  test  history  of 
the  six  main  impellers  is  plotted.  In  the  original  study  the  data  was  hand  plotted, 
and  later  incorporated  as  a routine  in  the  MSFC  Diagnostic  Data  Base  Program.  The 
output  is  shown  in  Figures  41  to  46.  Listed  below  are  the  main  impeller  serial 
numbers,  HPOTP  serial  number,  last  test  and  date  of  test. 


Main  Imoeller  S/N 

HPOTP  # 

Last  Test 

Date 

2427800 

2412 

A3-207 

7-21-83 

3134124 

231  IRl 

A 1-405 

2-4-83 

3334446 

2208R1 

Al-432 

1-4-84 

3135444 

0210R2 

A3- 192 

3-2-83 

7326066 

0307 

A2-376 

11-11-83 

7363708 

911  IRl 

A3-191 

2-12-83 

It  is  beyond  the  scope  of  this  study  to  fully  investigate  the  variation  in  synchron- 
ous vibration  levels  from  test  to  test,  pump  build  to  pump  build  or  verify  the  six 
main  impellers  were  the  sole  contributor  to  the  bi-modal  distribution  of  the  HPOTP 
synchronous  vibration  data.  However,  since  the  original  study  only  main  impeller 
S/N  7326708,  HPOTP  #0307,  Test  A2-375  and  A2-376  (Figure  46)  has  deviated  from  the 
pattern  of  all  the  identified  main  impellers  being  associated  with  abnormal  syn- 
chronous vibration  levels.  For  reference,  the  mean,  Ic,, , 2o  , and  3a  vibration 
levels  for  normal  turbopump  operation  are  as  follows. 

HPOTP  PBP  RAD  Spatial  Average  Synchronous 


10Q% 

104% 

109% 

Mean 

1.89  Grms 

2.09  Grms 

2.15  Grms 

Std  Dev 

1.12  Grms 

1.05  Grms 

1.03  Grms 

Mean  + 1 a 

3.01  Grms 

3.14  Grms 

3.18  Grms 

Mean  + 2 a 

4.13  Grms 

4.19  Grms 

4.21  Grms 

Mean  + 3 a 

5.25  Grms 

5.24  Grms 

5.24  Grms 

The  vibration  levels 

for  104%  power 

level  are 

shown  on 

each  plot  to  illustrate  the 

deviation  from  the 

expected  normal 

operation 

value  and 

relationship  to  the  sigma 

levels. 
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Figure  41.  Synchronous  Vibration  Test  History, 
Spatial  Average  PBP,  Impeller  S/N  2427800 
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Figure  44.  Synchronous  Vibration  Test  History, 
Spatial  Average  PBP,  Impeller  S/N  3135444 
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Figure  45.  Synchronous  Vibration  Test  History, 
Spatial  Average  PBP,  Impeller  S/N  7363066 


